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THE SCATTERING OF ELECTROMAGNETIC WAVES 
BY A CORRUGATED SHEET! 


T. B. A. SENIOR 


ABSTRACT 


The physical optics method is used to determine the scattering of a plane 
wave by a perfectly conducting sheet having sinusoidal corrugations. The 
only approximation concerns the current distribution on the sheet and the 
scattered field then appears as a spectrum of plane waves whose amplitudes are 
given by a simple integral expression. 


1. INTRODUCTION 


The problem of scattering by a corrugated sheet is one which is easy to 
formulate in mathematical terms but extremely difficult to solve, and in spite 
of its importance in many applications (for example, the reflection of low- 


frequency radio waves from an idealized sea) an exact solution has yet to 


be found. Indeed, very little progress has been made since the original treat- 
ment by Lord Rayleigh (1878) and the approximations which he proposed 
have been the basis for much of the later work. 

A direct attack upon the problem starts by assuming an expression for 
the scattered field as an angular spectrum of plane waves, the amplitudes of 
which have to be chosen to satisfy the boundary conditions. This leads to a 
single equation, valid at all points of the sheet, from which to determine 
the unknown amplitudes. It is here that approximations come in and some 
of those which have been suggested are referred to in Section 2. 

To avoid the difficulties associated with the series formulation several 
theories have been suggested which do not use the above equation as their 
starting point. Two of these will be mentioned because of their affinity to the 
physical optics approach, and both have been considered by LaCasce and 
Tamarkin (1956) in connection with the reflection of sound waves at an 
air—water interface. 

The first is due to Eckart (1953) and is a scalar theory. The source of 
radiation is at a finite distance with the field confined to a narrow two-dimen- 
sional beam, and the scattered field ¢° is then expressed as a Helmholtz 
integral involving ¢° and its normal derivative on the sheet. Although these 
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are not independent and cannot be specified independently of one another, 
the crux of the method is to replace both ¢° and 0¢°/dn by the values appro- 
priate to the incident field alone. A further assumption is to regard the reflec- 
tion as taking place at an equivalent flat surface and the scattered field then 
appears as a spectrum of waves whose amplitudes have a simple Bessel 
function form. 

The second theory (Brekhovskikh 1952) is also based on the Helmholtz 
integral. Here again the crucial step is the determination of one boundary 
condition having specified the other, and Brekhovskikh estimates 0¢°/dn by 
assuming that the incident field is scattered by each element of the corrugated 
sheet as though undergoing ordinary specular reflection. As such the method 
is similar to physical optics, but the advantages which should accrue are 
sacrificed by an ill-advised use of approximate techniques in evaluating the 
resulting integrals. 

The physical optics method itself is described in Section 3. The essential 
approximation is to replace the current distribution on the sheet by twice 
the tangential component of the incident magnetic field and the scattered 
field is then given by a double integral over the surface of the sheet. By a 
careful use of a steepest descent technique, the integrals can be evaluated 
in such a way as to preserve all the features of an exact solution, and some 
of the consequences are examined in Section 4. 


2. SERIES FORMULATION 


The sheet is assumed to be perfectly conducting and to have the equation 
(1) y = acos Kx 


where a is the amplitude of the corrugations and 27/K is the period. It is 
infinite in both the x and z directions and can be visualized as a tin roof 
whose corrugations are parallel to the z axis of the co-ordinate system. At 
any point (x, y, z) on the sheet the unit vector normal is 


1 
V/{1+(aK sin Kx)*}° 
The simplest case to consider is that in which the plane wave is incident 
normally, corresponding to propagation in the negative y direction. For ease 
of comparison with results obtained for the analogous problem in the theory 
of sound, it is convenient to take the magnetic vector entirely in the z direction 
and hence along the corrugations. The incident field is then* 


(2) n = (aK sin Kx, 1, 0) 


(3) E' = (E;,0,0), H'= (0,0,H:) 
with 
(4) E} = exp(—iky) and H} = Y exp(—iky) 


where Y is the intrinsic admittance of free space and k is the propagation 
constant. 


*Throughout this paper, m.k.s. units are employed and a time factor exp (—iwt) suppressed 
throughout. 
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Concerning the scattered field, the symmetry of the problem requires 
that its magnetic vector also lie in the z direction. If this is written as an 
angular spectrum of plane waves each of which is propagated at some angle 
6, to the positive y axis, we have 


Ht = Y > A, exp{ik(x sin 0,+y cos 0) } 
which can be inserted into the boundary condition 


| 
én on 





at the surface y = acos Kx to give 


(5) 7 A, (aK sin Kx sin 0,+cos 0.) exp[ik{x sin 0,,+ (1+ cos 0,,) a cos Kx}]=1. 


Because of the nature of the sheet, however, equation (5) must be periodic 
in x with period 27/K, and this implies that 


ksin 0, = mK 
where m can take all integer values, positive, negative, or zero. Hence 
k cos Om = ++/{k?—(mK)?} 


and since the positive sign must be chosen in order to comply with the con- 

dition of outgoing radiation at infinity, the equation from which to determine 

the A,, is 

(6) >> A,[amK* sin Kx+~/{k°—(mK)’}] exp(imKx+ ip» cos Kx) = k 
m=—co 

where yp, = alk++/{k?—(mK)?*}]. This is valid for all x, and, in theory at 

least, is sufficient to enable the A,, to be calculated. 

Once the A, have been found, H is given by 


(7) Hi = VY >> An explimKx+in/ {k’— (mK)"}y] 
m=—aco 

and the components of the electric vector can then be obtained from Maxwell's 
equations. For an incident field of the above polarization the entire problem 
can be expressed in terms of the single component H, and this is directly 
analogous to the velocity potential in the theory of sound. If E and H are 
interchanged it is still true that one component suffices (indeed, the equation 
for the A,, is somewhat simpler in this case), but the solution no longer applies 
to the scattering of sound by a rigid surface. 

Although the exact solution of equation (6) is still unknown, a certain 
amount of information about the scattered field can be obtained without a 
knowledge of the A,,. If, in equation (6), m is replaced by —m and x by —x, 
the equation is unaltered except for the fact that A_,, appears instead of 
A, Hence A_m = Am and, from equation (7), the scattered field is seen to 
consist of a fundamental wave propagated in the direction of the y axis 
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(specular reflection by the equivalent flat sheet), together with an infinity 
of pairs of higher-order modes. Each pair is composed of two waves of equal 
amplitude and phase whose directions of propagation are at angles +6, to 
the y axis, where 

cepsieiit i cincrss 
VIE = (mK) ' 
As m increases these directions become more nearly horizontal until, finally, 
the value of m is such that k? < (mK)?. The corresponding mode and all 
those of higher order now travel horizontally and, being attenuated in the 
y direction, are surface waves whose influence is confined to the immediate 
vicinity of the sheet. Similarly, if attention is concentrated on a particular 
mode, the mth say, and if k is slowly decreased, a frequency will be reached 
at which propagation in the y direction ceases. The appropriate frequency 
has k = mK and is often referred to as the ‘cutoff’ for that mode. All except 
the fundamental exhibit this behavior. 

Unfortunately, only a limited amount of information can be obtained in 
this manner and sooner or later it becomes necessary to determine the ampli- 
tudes A,,. The main difficulty with equation (6) is due to the form of the 
functions which multiply the individual terms. If these functions were ortho- 
gonal the equation could be inverted and though the presence of the expo- 
nentials prevents this from being done, the exponentials themselves can be 
expanded in terms of an orthogonal set of functions. The most direct way of 
doing this was proposed by Lord Rayleigh (1878). Since the functions 
exp (ipKx), p = 0, +1, +2,... are orthogonal over any 27/A range of x, 
the expansion 


0. tan: iS Ee Ooi 


oc 


exp (ium cos Kx) = > 1” Jy(um) exp(ipKx) 
P=—a 
can be used to generate an infinite set of equations from which to calculate 
the modal amplitudes. To proceed any further, however, approximations have 
to be introduced and a typical one is to assume A’/k is small, with ak of order 
unity. This leads to an expansion of A, in ascending powers of A/k. Similar 
results have been derived by LaCasce and Tamarkin (1956). 

Once the need to approximate has been admitted, a more direct approach 
is to modify the exponentials in the original equation (6). If, for example, 
it is assumed that K is sufficiently small for (mK)? to be neglected in com- 
parison with &? for all values of m, equation (6) can be inverted to give a 
recurrence relation for the A,,. Whilst it is clear that the assumption on which 
this is based will ultimately fail no matter how small K is, the solution is in 
gocd agreement with the results of other approximate methods. 

All the methods described above impose similar restrictions on the magni- 
tudes of K/k and ak, and in essence the requirements are K Kk, ak = O0(1). 
These conditions are aimed at minimizing the influence of the cutoff modes 
and reducing the effect of any one mode on those of lower order. In some 
cases it is actually assumed that all cutoff modes can be omitted from equation 
(6), and since the sum then contains only a finite number of terms, the equation 
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can be solved ‘exactly’. The assumption, however, is quite invalid. Just 
because a mode does not carry energy away from the sheet its amplitude 
need not be zero, and cannot be neglected in an equation which refers to 
the very place where such a mode may predominate. 


3. PHYSICAL OPTICS 


The only approximation inherent in this method is one which is implied by 
an assumption about the current distribution excited on the surface of the 
sheet. The exact current is determined by the total field, and this is unknown, 
but according to physical optics a reasonable approximation can be found 
by reference to the incident field alone and is 


J = 2(a A BB. 


In other words, the distribution at every point of the surface is the same 
as if the reflection were taking place at an infinite tangent plane. 
For the plane wave whose components are given by equations (3) and (4), 
the postulated current distribution is 
ee 2Y exp(— ak cos Kx) 
J (, ~eKk aekz, 9) —- SS he 
( ) V/{1+(aK sin Kx)*} 
This generates a scattered field 


H’ = -1 fs A vPdS 
4nr 


where P = {exp (1kR)}/R, R being the distance of the receiver from a variable 
point on the surface. If the co-ordinates of the receiver are (xo, yo, 0), then 


R = V{(xo—x)?+(yo—a@ cos Kx)?+27} 


and since 


, _ (xe-* yo—a cos Kx a 
ve ( ke" R 70 oR’ 
it follows that 


H® = (0,0, 772) 


where 
(83) Ht = Y r f ° (a cos Kx— yo) +(x—xo0)aK sin Kx 
—_ oS oe ae R V{1+(aK sin Kx)*} 
oP ; 
x aR exP (—1ak cos Kx) dx dz 
with 


os (a! ) exp(kR) 
nS oe x 


There are two features of this result which should be noticed. If xo is 
increased by 2mx/K, where m is an integer, the integral can be reduced to 
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its original form by writing x’ = x—2m2/K. The expression for H is there- 
fore periodic in x9 with period 27/K, a property which is essential in view 
of the corrugated nature of the sheet. Secondly, when a = 0, 


e el fs oP 
HH, = 3 ari @ 


and by introducing polar co-ordinates the integrals can be evaluated exactly 
to give 

H: = Y exp (ikyo), 
as required for reflection from a flat sheet. 

When a is not zero, however, the integrals can only be evaluated by 
approximate techniques. The one parameter at our disposal is yo and since 
our main interest is in the field at large distances, yp can be assumed large 
compared witlf everything except the lateral dimensions of the sheet. 

Equation (8) can be written as 





Hi= Y ‘ ~ (a cos Te ee mks exp(—iak cos Kx) I(x) dx 
where 
Ler 
I@) = x J. Rar” 
with 


R= V(22+6?) and b = V{(xo—x)?+(yo—a cos Kx)?}. 


The z integration will be carried out first. Since 6 is a constant in this 
integration, the substitution z = b sinha reduces I(x) to the form 


ee exp (7bk cosh a) 
I(x) T J (i 6b cosh -) b cosh a da 


b > |yo—a cos Kx, 


But 


and the right-hand side can be made as large as we please by choosing yo 
appropriately. This makes possible an asymptotic evaluation of I(x) by the 
method of steepest descents. The saddle point is a = 0 and hence 


I(x) ~ - =. exp{i(bk—1/4)}, 


which is accurate to the first order in yo. By insertion into equation (8), we 


now have 


” (a cos Kx—yo)+(x—x0)aK sin Kx 


eee. - 
(9) Hp=—- “FZ. exp(—im/4) /{1+(aK sin Kx)*} 


—oo 


exp{ik(b—a cos Kx)} 
pees ais dx. 


Xx 
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Rather than attempt to carry out the x integration immediately, it is 
more convenient to find the Fourier transform with respect to xo. The trans- 
form is defined as 


A = f exp (ivxo)H? dxo 


and consequently, from equation (9), 


L(x) 


_sV114 GK sin Kx)*) P(— tak 008 Kx) de 


(10) H = 
where 
(141) L(x) = - ys exp(—im/4) : { (a cos Kx—yo) +(x—X0) aK sin Kx} 


exp (ivxo+71bk) 
bb 


This last integration can be performed by making the substitution 


x9 = x+c sinh B 


x dx». 


where ¢ = |yo—a cos Kx}. Since 6 = c cosh 8, the exponential factor in the 
integrand becomes 
exp 7 {yx-+d cosh (8+ y)} 


with dcosh y = kc, dsinh y = vc, and d = cv/(k?—v’*). Here again c, and 
therefore d, can be made as large as we please providing v ¥ k, and if this 
particular value of v is excluded for the time being, the integral can be treated 
by the steepest descents method. The saddle point is 8 = —y and to the 
first order in yo, 

v 


an 
— 7ivaey sin Kx¢ , 


which can be substituted into equation (10) to give 


L(x) ~ exp{i(vx+d)} 1 


= Vexn{; 2) ( Late/ V(k*—»')} aK sin Kx 
(12) A = Yexpfiyo/(k’—»’)} : Vil+ GK sin Kx)"} 


X explivex—ia{k++/(k’—v’)} cos Kx] dx. 


Having reached this stage, it is a simple matter to express H$ in a form 
which brings forth the dependence on the co-ordinates x9 and yo. For this 
purpose equation (12) is written as 


(13) A = Yexp{iyo V(k’—v’)} f F(Kx) exp(ivx) dx 


where 








(14) F(Kx) = A= ear = Ss expl—ia tht V a | 
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Since F(Kx) is a periodic function of x of period 27/K, it can be expanded 
as a Fourier series 


co 


F(Kx) = >> A,,(v) exp(imKx) 


m=—oo 


with amplitudes 4,,(v) given by 
K mn /K 

(15) Ax ly): == J F(Kx) exp(—imKx) dx. 
2a —/K 

If this expansion is inserted into equation (13), we have 


bi Pest i DRE he J expli(o-+mk)x| de 


m=—a@ 


and the integral is simply a delta function having unit area when vy = —mK 
but zero otherwise. Hence 


2 «Tatiana Hat Adi 


m=—oc 


and application of the inverse transform now gives 


(16) Ht = VY >> A,(—mK) explimKxo+i V{k?— (mK)"} yo] 
(cf. equation (7)). This is a representation of the scattered field as an angular 
spectrum of plane waves. The form is precisely that which would be expected 
of an exact solution and, from equations (14) and (15), the amplitude of the 
mth mode is 

L ("1+ [mK/ V{k*—(mK)"}] aK sinx 


(17) An = 5- & Vii+(aK sin x)'} p(—imx —ipm cos x) dx. 


It will be observed that for a = 0 


7" 


An = on _exp(imx) dx, 
which reduces to unity when m = 0 and to zero otherwise. 


4. DISCUSSION 


In evaluating the above integrals the only approximation which has been 
made is to take kyo large compared with unity. The expression for H$ then 
appears as a spectrum of plane waves propagating in the directions 0p, 
lm| < k/K, together with an infinity of waves which are attenuated in the 
direction normal to the sheet. For large yo, however, the amplitudes must 
be the same as would be obtained by rigorous evaluation of the integrals in 
(8) and the fact that the entire dependence on yo is through the complex 
exponentials representing the plane waves now shows that the solution is 
valid for all yo. In other words, since the ‘error’ terms are inverse powers 
of yo which cannot be expressed as exponentials (so as to produce a modification 
to the A,,), the solution must correspond to an exact evaluation of the original 


integrals. 
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Let us consider the behavior of the spectrum as a function of k. From 
equation (16) the number of modes carrying energy from the sheet depends 
upon the values of A and &. If attention is concentrated on the mth mode, 
where m > 0, then for k > mK the mode propagates in the direction specified 
by the angle @,, and as k decreases, 0, increases towards 7/2. The amplitude 
A, also increases due to the presence of the factor /{k?—(mK)?} in the 
denominator and since equation (17) is certainly valid for all k # mK, the 
increase must continue until, ultimately, A, becomes infinite. This occurs 
when k = mK, the cutoff value for the mth mode. The wave is now travelling 
horizontally and any further decrease in k introduces an attenuation in the 
y direction which confines the wave to the immediate vicinity of the corru- 
gated sheet. In effect we have a surface wave whose amplitude decreases with 
k but is not necessarily small. 

In a formal sense the behavior when k = mK is still unknown because of 
the assumption v # k in evaluating the expression for L(x). Continuity, how- 
ever, obviously requires that 4, be infinite there and this can be confirmed 
by putting v = k in equation (11). The substitution x» = x+c¢ sinh 8 then 
gives 
L(x) = 4 exp{i(kx— 71/4)} i (i+eK sin Kx sinh 8 a 3 


and since sinh 8 is an odd function of 8, 


Ete). = 4/e -exp{i(kx— 1/4)} t oy oo 8 
ae 2 
4/ be exp{i(kx— 71/4)} ‘ r(2) \ 


for large yo, where I denotes the gamma function. The rest of the analysis 
is as before with the result that, at cutoff, 


—1mx —tak cos x) 


se exp(— tr 4) e rt j2 exp : Si 
(18) An = “Veyo ir 5 Vi{l+(aK sin x)’} ~ 


The dependence on yo through a positive power is the only means by which 
4,, can become infinite when k = mK, since yo is the one variable (and large) 
parameter appearing in equation (11) and L(x) cannot be infinite per se. As 
cutoff is approached, yo must be taken larger and larger in order to fulfill 
the condition vow (k2—v*) > 1 in the evaluation of L(x) and hence, in (18), yo 
must be treated as infinite. 

These conclusions only apply to the amplitude factor A, and for other 
field components the cutoff behavior may differ from the above. With the 
incident field of Section 2, the non-zero components in the scattered field 


9 


are H3, E8, and E&’. From equation (16), we have 


a Jf a r\2 9 a, @ 
tie ee VIE MKT 4, explimKxo+i V {k= (mK)*| yol, 
ee “ mK ee . 2 r\ 2 
t= > Am explimKxo+i V {k'— (mK) } yo), 


m=—co 
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and whereas E§ is similar to H$ in that its mth mode becomes infinite at 
the cutoff frequency, the corresponding mode in the expansion for E$ remains 
finite when k ee mK and its amplitude has the limiting value 


iaK sin x sin mx a 
(19) 7 de ieee exp(—imaK cos x) dx. 


This difference is a direct consequence of the boundary conditions and pro- 
vides an explanation for the discrepancy between the formulae of Eckart 
and Brekhovskikh near to cutoff (see LaCasce and Tamarkin 1956). In both 
theories conditions are imposed on the function and its derivative at the 
conducting sheet, and since the subsequent analysis is approximate, the 
amplitude factors will depend on the relative weights attached to these 
conditions in the derivation of the solution. 

All the field components are determined by the amplitudes A, and an 
expression for A,, as a finite integral is given in equation (17). Unfortunately, 
no method for evaluating this has been found. The integral is a function of 
the three parameters a, K, and k, and the case which would be the most 
interesting to study is that in which aK is not small compared with unity. 
The theories referred to in Sections 1 and 2 then fail and though it is difficult 
to justify the use of an optics current if aK is as large as k/K, the integral 
has been evaluated numerically for a few values of the parameters satisfying 
this condition. The results are in good agreement with the experimental 
measurements of LaCasce and Tamarkin. 

If aK < 1 it is relatively easy to obtain an approximation to the integral. 
The denominator in the integrand can then be expanded in ascending powers 


of aK to ag 


Re. 
-+f ome me 30K sin x)*40(ak*)! 


X exp(—imx—iptm cos x) dx 


and since 


s.. 
20 e 


f exp(—imx—ipm cos x) dx = (—1)”" Jm(um), 


it follows that 
ae ee (aK) 
(20) An Por ( 1) | sat (mK)? y Tis (tm) — Am 


X {(m— 1) Sess (um) + (m+1) Susi | 


where terms of order (aK)* have been omitted. Of particular interest are the 
amplitudes of the fundamental and first harmonic, and these are 


Ay= Je(2ak) — 9 J,(2ak), 


— ek Ilut 


k 
A, = ~ itz ~K") Ji (41) — 
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which are similar to the results of Lord Rayleigh (1878). However, the case 
aK <1 is not a good one for demonstrating the advantages of the present 
solution. In the absence of an exact evaluation of the integral for A,, it 
would be most desirable to have an expansion corresponding to (20) but 
appropriate to large aK. This would permit a more complete comparison 
with the experimental data, but no way of adequately treating the case 
aK > 1 has yet been found.* 


5. CONCLUSIONS 

Physical optics is known to be an effective method for treating scattering 
by bodies whose dimensions (including radii of curvature) are large compared 
with the wavelength, and this is particularly true when the entire surface is 
illuminated by the incident field. The problem of a corrugated sheet satisfies 
these requirements. 

The optics method as such is, of course, not new, but in the past there 
has been a tendency to squander some of the advantages by a lack of care in 
the evaluation of the resulting integrals. Inaccuracies which then appear may 
arise as a consequence of the evaluation and should not be automatically 
ascribed to a failure in the approximation to the current. 

The main feature of the present paper is the evaluation of the optics 
integral. Although a steepest descents technique is used on two separate 
occasions, these are such that the only condition is upon the distance of the 
receiver from the sheet (kyo > 1), and since the solution is a spectrum of 
plane waves, this must be equivalent to an exact evaluation of the integral. It 
is interesting to note that even the cutoff modes, which would not be expected 
to contribute to the field at infinity, do appear in the final solution. 

The amplitudes of the plane waves are given by a simple integral over a 
finite range, but no method for evaluating this has been found. For aK < 1, 
however, the integrand can be expanded to give a series expression for the 
amplitude, and this is in good agreement with the available experimental 
data. For certain larger values of aK the integral has been evaluated numeri- 
cally for different ak, and the fact that these results also agree with experiment 
suggests that the solution has a wide range of applicability. 
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*Since this was written a series evaluation has been obtained which is convergent (and 
useful) for all values of aA. This will be described in a future paper. 








THERMAL CONDUCTIVITY OF POWDERS! 


M. J. Lausitz 


ABSTRACT 


The effective thermal conductivity of several powders of reasonably uniform 
particle size, and one powder of graded particle size, was measured in the tem- 
perature range of 100° C to 1000° C. Most of the results can be very satisfactorily 
explained by assuming that the effective conductivity of the powder is composed 
of two parts, the first part being a two-phase gas~solid conductivity, and the 
second an equivalent radiation conductivity. 


INTRODUCTION 

There has been in the past a continuous interest in the effective conductivity 
of multiphase media, particularly of the two-phase gas-solid medium. The 
conductivity of such a two-phase medium depends on which of the two is the 
continuous phase, and which is what may conveniently be called an obstacle. 
The case of dense ceramics, where the ceramic may be considered as the 
continuous phase, has been studied in detail by Francl and Kingery (1954). 
They compared their results with several of the available theoretical formulae, 
and found reasonably good agreement. Powdered ceramics, however, where 
air or another gas is the continuous phase, have not yielded particularly con- 
sistent results, and the agreement with the theoretical expressions derived 
from various models has not been spectacular. The chief shortcomings of the 
theoretical expression are due to the necessary simplicity and artificiality of 
the models for which they have been calculated. The difficulty here lies in 
two directions: primarily, the construction of a model that represents reason- 
ably well the two-phase system under consideration, and secondarily, the 
calculation of the effective conductivity of that model. The first part is 
difficult, particularly if irregular particles of varying size are considered. The 
second is impossible except for very simple geometry, or under some simplifying 
though highly unreasonable assumptions. So far, only one accurate expression 
for the effective thermal conductivity has been obtained. This was calculated 
by Lord Rayleigh (1892) for a cubical array of uniform spheres, and his result 
can be put into the following form: 

, K _ (2+v)/(l—v)—2p 
5 Kn (2+»)/(I—») +p 
where K = the effective conductivity of the two-phase medium, 
the conductivity of the continuous phase, 

y SU hh. 
K, = the conductivity of the obstacles, 

p = the fractional volume of obstacles. 


. 
I 


This expression was later generalized by Burgers (1919) to include arrays of 
arbitrary ellipsoids in a cubical arrangement. 
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Fairly recently the cubical array of spheres has been considered again, with 
certain simplifying assumptions. Deissler and Eian (1952) envisaged such an 
array of spheres with the heat flow lines being parallel straight lines. This 
is equivalent to assuming that the conductivity of the components of the 
medium is vanishingly small in directions perpendicular to the heat flow. 
Their calculation is done for a closely packed structure, i.e. p = 2/6. It can, 
however, be easily generalized for arbitrary » < 2/6, and has then the form 
si K _ x {(/r-1)R=Inft+ (oD yy aR 

K,. 2 (1/y—1)° 4 
where R = (6p/7)'” is the radius of the spherical obstacles. Woodside (1958), 
on the other hand, takes the cubical array of spheres and assumes that the 
isotherms are planes perpendicular to the direction of heat flow (equivalent 
to assuming an infinite lateral conductivity). He obtains the following ex- 


pression : 
i Hf is (<=!) fort) 
2a a—l 


4 1/2 
“— [ = ne ; 


There is one further model to consider, proposed by Russell (1935). He 
calculated the conductivity of cubes in a cubical array, assuming again that 
isotherms are planes perpendicular to the direction of heat flow. For this 
model he obtains 


(3) 


ee, 
Aly 
B | 

Il 


where 





: Kp 4p") 

“ Kn v(p’*—p)+(1—p"* +) * 

It is worth noting that for large » Russell’s equation reduces to 
F Moe Se 

(5) Km = i-s” 


and this is identical with the form obtained for large v if one treats the cubical 
array of cubes in the straight flow-line approximation. 

It may well be asked why one should consider any but Lord Rayleigh’s 
equation, which is the only exact one. The answer to that lies in the dual 
difficulty of the theoretical calculations mentioned before. Lord Rayleigh 
provides an exact mathematical calculation, but his model is so rigid and 
artificial, permitting the spheres to be located only in the center of their 
unit cubical spaces, that it radically departs from powders encountered in 
practice. The plane isotherm models, on the other hand, though their mathe- 
matical treatment is unrigorous, permit an arbitrarily irregular location of 
obstacles in their unit cubical spaces, and hence are more realistic. Further- 
more, Lord Rayleigh’s model (and the other models using spheres) limits the 
value of p to less than 2/6, while Russell’s model permits the whole range 
of p from zero to unity. 
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Figures 1 and 2 show some of the values of K/K,, obtained from the 
equations (1), (2), (3), and (4). The conductivities K(L.R), K(D and E), 


K(W) 


K(W) 





K(D+E) 


K/Km 


<---> 





; : [2° 36 a0 100 1000 
P Ko /Km 


Fic. 1. Variation of effective conductivity with obstacle volume fraction, for a fixed value 
of solid-to-gas conductivities, Ko/Km = 100. 

Fic. 2. Variation of effective conductivity with ratio of solid-to-gas conductivity for fixed 
obstacle volume fraction, p = 7/6. 


K(W), and K(R) refer to the conductivities calculated from equations (1) 
to (4) respectively. It is seen from these two figures that A(L.R) and K(R) 
do not differ very much from each other. This is true for all values of p and 
K./Km. On the other hand, K(W) and K(D and E) deviate considerably 
from K(L.R), particularly for large p and K,/K,. This deviation may con- 
demn these models from the mathematical point of view. From the practical 
point of view they may still be useful, however, as the final selection of the 
model is in this case a purely empirical matter. 

To determine which of these models, if any, best represents the actual 
powders, a set of experiments was undertaken with powders of various 
materials, various porosities, and various grain sizes. The effective thermal 
conductivity of these powders was measured in the range of 100° C to 1000° C, 
and the results compared with those predicted by the theoretical equations. 


EXPERIMENTAL PROCEDURE AND RESULTS 


The experimental arrangement is illustrated in Fig. 3. Essentially the 


apparatus consists of a long cylindrical furnace, inside which there is a 
centrally located heater 4 mm in diameter and 30 cm long. Two Pt/(Pt-10% 
Rh) thermocouples are held rigidly at radial distances of 5.6 mm and 10.1 mm, 
the tip of the outside thermocouple being about 1 cm above the tip of the inside 
one. The furnace was filled with the required powder, tapped, and filled 
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Fic. 3. The experimental arrangement for the measurement of the thermal conductivity 
of powders. 
again. This was continued until the powder no longer subsided on reasonable 
tapping. From the weight of the powder put in, and the known volume of 
the furnace, the porosity of the powder under experimental conditions was 
calculated. 

The following procedure was employed to measure the thermal conductivity. 
With no power through the central heater, the input into the furnace was 
adjusted until the outside thermocouple registered a required temperature 
To. At equilibrium, the radial gradient A7(Q) was measured. Then the power 
through the central heater was switched on, and again the furnace adjusted 
until the outside thermocouple read approximately 7». At equilibrium, the 
gradient A7(P), corresponding to the central heater power input of P watts 
was measured. The effective conductivity of the powder at the temperature 
Tot+[AT(P)]/2 was then calculated from the formula 


oh on aries) = In Ts ee. : Fr 
K( tet 2 ~ Qel AT(P)—AT(O) 
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where / is the length of the central heater in centimeters, rz and 7; are the 
radial distances of the outer and inner thermocouples respectively, and K(7) 
is the thermal conductivity of the powder at temperature 7 in watt/cm °C. 

The powders studied are listed in Table I, together with their pertinent 
properties. The powders marked ‘“‘D” were composed of dense crystalline 


TABLE I 


Properties of the powders 


Powder Volume fraction 


Nature of Particle diam., density, of particles, 
Powder particles cm g/cc r 
MgO D(10-20 Dense crystals 0.140+0.060 1.60 0.469 
MgO D(35-50) Dense crystals 0.040+0.010 1.70 0.475 
MgO D(100-200 Dense crystals 0.011+0.004 1.54 0.430 
Al.O; D(35-50) Dense crystals 0.040+0.010 1.82 0.475 
Al,O; B(35-50) Porous bubbles 0.040+0.010 1.15 0.29 
ZrO. B(35-50 Porous bubbles 0.040+0.010 2.02 0.33 
Al.O; R Graded size 1.00 0.25 


particles of porosity no larger than 4%, and of irregular shape. The powders 
marked ‘‘B”’ were in the form of bubbles, with the ceramic shells surrounding 
air. Their shape was nearly spherical, and surfaces smooth. The numbers in 
parentheses behind the powder material indicate the mesh of the sieves 
limiting the particle size. The Al,O; R powder was a reagent grade alumina 
with a wide range of particle size (Table II). Its low density indicates that 


TABLE II 
Particle distribution of Al,O; R 


Powder retained on 


Mesh sieve, % 
50 0 
100 13.1 
200 58.2 
325 19.9 
Through 325 8.8 
100.0 


the particles were highly porous. The experimentally obtained conductivities 
for these powders are shown in Figs. 4 to 8. 

The absolute accuracy of the measurements is not very high. The maximum 
deviation from the average observed in the measured conductivity for four 
separate fillings of the furnace with one and the same powder was +12%. 
The inaccuracy due to thermocouple positioning and possibly lack of con- 
centricity of central heater and furnace tube is estimated at +10%. Thus 
one has to put an uncertainty of +20% on the absolute value of the measure- 
ments. The relative accuracy, that is the relative values of thermal con- 
ductivity of the different materials is much better, of the order of 410%, 
This is also the accuracy one would expect for the relative values of the 


conductivity of one and the same powder at various temperatures. 
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ANALYSIS AND DISCUSSION 


The analysis of the dense particle powders being more straightforward, it 
is undertaken here first. In Figs. 4 to 6 are shown the results for the four 
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Fic. 4. Thermal conductivity of MgO D(1l0-20) 
@@@ [experimental points for central heater power input of 4 watts cm. A&A Experi 
mental points for central heater power input of 2 watts/cm, 

Fic. 5. Thermal conductivity of MgO D(35-50) em and ALLO; D(35-50). XXX Expert- 
mental points for AlOs. @@@ Experimental points for MgO 

Fic. 6.) Thermal conductivity of MgO D( 100-200). @@@ Experimental points 
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dense particle powders, together with the calculated effective conductivities 
K(L.R), K(R), and K(W). In these calculations the values for dense ceramics 
were taken from Kingery et al. (1954) and for air from Russell (1935). The 
calculated conductivities do not correspond to the experimental ones, either 
in absolute value or in their temperature variation. All the experimental 
conductivities increase more rapidly with temperature than any of the models 
predict. This points to the omission of some temperature-dependent factor 
in the theoretical analysis. Such a factor is, of course, radiation. 

The inclusion of radiation into any of the previously mentioned models is 
impossible without some structural changes. For, as these models stand 
now, continuous paths are open to the radiant energy. To obtain finite radiation 
path particles must be randomly positioned in their respective unit cells—this 
is possible only in the models of Woodside and Russell. As Woodside’s con- 
ductivities are already larger than the experimental ones, Russell’s model 
alone is left for modification to include radiation. Using that model, an 
expression for the contribution of radiation is derived in the appendix. The 
final formula for the effective conductivity is then 


(6) Ky, = K(R)+4oT"et (1—p**+p*”) 


ry <<,” 
(R) a K, (2/3 2/3 
vip —p)+(l—p +p) 

is the conductivity calculated by Russell. 

e is the emissivity of the ceramic (given in High Temperature Technology 

(1956)), 

a is the linear dimension of the particles, 

p is the fractional volume of the particles, 

o is the Stefan—Boltzmann constant. 

The conductivity A; accounts well for the shape of the temperature curve, 
but not for its absolute value. One can get a very good agreement with the 
experimental data, however, if one uses the formula 


(7) Ky = 2K(R)+407%e ; (1—p**+p""). 


The reason why A(R) has to be doubled can be seen qualitatively as 
follows. In actual powders there are densely packed regions where particles 
touch each other, and empty air pockets. The densely packed regions con- 
tribute chiefly by conduction while the air pockets contribute by radiation. As 
shown in the appendix, the air pockets have roughly a length of a/po, where 
a is the linear dimension of the particles, and po the average particle volume 
fraction for the aggregate. But then the p value for the densely packed region 
is not Po but po'*, and the equivalent conductivity of the powder becomes 
K(R, po!) . po?*, where A(R, fo'”’) is the Russell’s conductivity calculated for 
pb = po'*. For the powders used in these experiments 


K(R, p'). p?? = 2K(R, p), 


which gives a rough reason for the factor of two used in equation (7). 
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One can now determine whether equation (7) applies to the other powders 
tested. The blown powders can be considered in the following way. The 
radiation contribution will depend primarily on the size of the particles and 
their packing, and not on the true porosity. Thus the radiation contribution 
has to be evaluated on the basis of » = p1, where p; corresponds to the 
packing of the particles, and equals approximately 0.5. The conductivity, on 
the other hand, has to be evaluated for the true porosity, p ~ 0.3. Thus 


Ky = 2K(R, p)+4oT es (1 eT 
1 


where p, ~ 0.5, 
p ~ 0.3. 


Figure 7 shows that K,; accounts quite well for the experimental results. It 
can be pointed out here that for alumina and magnesia, the conductivities of 
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Fic. 7. Thermal conductivity for AlxO; B(35-50) and ZrO, B(35-50). @@@ Experi- 
mental points for AlxO;. XXX Experimental points for ZrOx. 
Fic. 8. Thermal conductivity for AlO; R. @@@ Experimental points. 


the solids are so high that equation (5), independent of Ko, can be used to an 
accuracy of 2 to 4%. This is not the case for zirconia, whose conductivity is 
relatively small. Thus, for a given particle size, the conductivities calculated 
for alumina and magnesia are the same, but different from that for zirconia. 

There remain the results for the graded alumina. Because of its very small 
average particle size one expects to have no significant contribution from 
radiation, and equation (7) reduces to 


(8) Ky, = 2K(R). 


a) 
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However, this does not explain the experimental results obtained (Fig. 8). In 
fact, the temperature variation of K(R), or of any other of the previously 
considered conductivities, is such that it cannot account for the graded 
powder. It may well be that such powder will require the construction of 
completely new models, but there is not as yet sufficient experimental informa- 
tion available to confirm that. Furthermore, some of the experimental work 
is rather confusing. For instance, Deissler and Eian published results on graded 
magnesia in argon, air, and helium for a p value of 0.58 in 1952 and 0.64 in 
1953. Their 1952 results for air and argon but not helium can be described 
quite well by equation (8), while their 1953 results are completely different 
from what is predicted by that equation. Particularly confusing is the curvature 
they obtained in 1953 for the magnesia in helium. 

Another interesting thing may here be mentioned. Figure 4 shows the 
results for the MgO D(10-20). The experimental points lie on two curves 
which meet at the lower and diverge at the higher temperatures. These two 
sets of points correspond to two different power inputs for the central heater, 
the higher conductivity being obtained from the large power input. This 
effect may be due to the transparency of magnesia to infrared, an effect which 
has been previously observed (see discussion after Russell 1935). Such trans- 
parency would increase the effective conductivity of the powder, the increase 
being a function of both temperature and temperature gradient. However, 
this contribution would only be significant in coarse powders, where there 
are relatively few reflecting grain surfaces. This is precisely what has been 


observed. 
SUMMARY 


The effective thermal conductivity of several powders of reasonably uniform 
particle size, and one powder of graded particle size, has been measured in the 
temperature range 100° C to 1000° C. The results for the powders with uniform 
particle size can be very well accounted for if the effective conductivity 
is assumed to be composed of two parts, the two-phase conductivity as calcu- 
lated by Russell (1935), arbitrarily doubled, and an equivalent radiation 
conductivity. For a large ratio of obstacle to medium conductivity, the formula 
for the effective conductivity can be put in the form 

a —2Km 4 407 a (1 pty gM), 
=F p 

The conductivity calculated by this formula does not agree with the measure- 
ments on graded powders. It is likely that for such powders a new model will 
have to be devised for the calculation of effective conductivity. 
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APPENDIX 


CALCULATION OF THE CONTRIBUTION OF RADIATION TO THE 
EFFECTIVE CONDUCTIVITY OF POWDERS 


The basic model is one given by Russell (1935). The obstacles are assumed 
to be cubes of length, a, located in cubical volumes of the medium of length, 
d. It is then assumed that the obstacles are randomly scattered in the cubical 
volumes, but in such a manner that the faces of the two cubes remain parallel. 
Such a scattering does not affect the thermal conductivity of the two-phase 
medium. Let the probability that radiation will pass a cube without hitting 
an obstacle be r, then 

a 2) 
¢ 
n= Ip: 
If ZL is the average length that radiation travels if it has not hit an obstacle 
in the first cube, then 


L = d+rd+r*d+... 
d d’ 


l-—r a” 


If p is the particle volume fraction, then 


i 
e* d° 
and 
a 
L= - 


The rate of heat transfer by that radiation is 
Q = 407%(d?—a")AT 
) ATd® 


px 
= 47 wi a L 


where T is the absolute temperature, ¢ is the emissivity of the obstacles, and 
AT/L is the gradient in the two-phase medium, assumed small compared 
to T. Therefore the contribution of that radiation to the effective conductivity 


1S 


tort §) = Aa te ; (1—p?"). 
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There will also be a contribution by the radiation that does hit an obstacle 
in the first cube of the medium. This contribution is difficult to calculate, 
but one can estimate it by assuming that the radiation that does hit an 
obstacle in the first cube travels on the average a distance of (d—a). Thus 
the rate of heat transfer by this process is 

a adK,Ky AT 
aK,+LK, d 
where 
K, = 4oT*eL 
and its contribution to the effective conductivity is 


= eee a 


d(aK,+LkK,) "1 K,+K,/p 


and as usually K,/p > K,, the contribution becomes 





pK, = p47 = 


The total contribution by radiation is then 


73 a 2/3 
4oT*e — (1—p?*?+-p*"*) 
p 
and the total effective conductivity 


3 a (1—p?*4 p*/*) 


K = K(R)+4oT'e > 


where K(R) is the two-phase conductivity calculated by Russell. 








THE PHOTODISINTEGRATION OF TANTALUM! 


R. W. Parsons? AND L. Katz 


ABSTRACT 


The cross section for the absorption of photons by Ta!®! has been deduced 
from measurements of the neutrons emitted when tantalum is irradiated in the 
X-ray beam of a 22-Mev betatron. The cross section has well-defined peaks 
at 12.75+.25 and 15.25+.25 Mev, and from this the intrinsic electric quadru- 
pole moment is deduced to be 5.7+1.3 barns. 


INTRODUCTION 

It has been suggested by Okimoto (1958) and Danos (1958) that the half- 
width of the photonuclear photon absorption cross section should be related 
to the nuclear eccentricity and hence to the intrinsic electric quadrupole 
moment. For sufficiently large eccentricities, the giant resonance is expected 
to show two distinct maxima at energies E, and EF, corresponding, on the 
hydrodynamical model of Goldhaber and Teller (1948) as developed by Danos, 
to oscillations of the proton and neutron fluids either parallel to the axis of 
rotational symmetry or normal to it. Danos predicts that, for a spheroidal 
nucleus with half-axes a and 8, 


(1) E,/E, = 0.911 a/b+0.089 


where a and bare related to the atomic number Z and the intrinsic quadrupole 
moment Q by the equation 


9 < - 
Q= -Z(a —b). 


If the separation E,—E, between the resonance energies is smaller than the 
half-widths of the component resonance curves, the total photon absorption 
cross section may show a single broad peak instead of two distinct ones. 
Nathans and Halpern (1954) found that the giant resonance half-widths 
for photon absorption by closed shell nuclei, for which Q = 0, are abnormally 
small. Katz and Chidley (1958) confirmed this and, in a study of 17 elements, 
found that the half-widths agreed reasonably well with those predicted by 
Okimoto (1958). Their measurements were not sufficiently precise to show 
evidence of double peaks, but such evidence has been found by Fuller and 
Weiss (1958), who studied tantalum and terbium, and Spicer et a/. (1958), 
who studied tantalum. Both these groups measured the neutron yield when 
a specimen is irradiated with X rays and deduced the photon absorption 
cross section from a knowledge of the X-ray bremsstrahlung curves and the 
neutron multiplicity. Both found two peaks in the giant resonance of tantalum, 
and from their separation the intrinsic quadrupole moment was estimated to 
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be 5.7+0.3 barns (Fuller and Weiss 1958) and 6.1+0.6 barns (Spicer 1958). 
Fuller and Hayward (1958) have recently reported measurements of the 
elastic scattering cross section for photons by tantalum, but although the 
results are suggestive the statistics are not good enough to provide convincing 
evidence for two peaks. 

The above results, however, are not confirmed by the recent work of Carver 
and Turchinetz (1959). They found the cross section o; for the reaction 
Tal'(y, 2)Ta'®® by irradiating tantalum with X rays and measuring the 
residual activity of the 8.15-hour isomeric state of Ta!®®. The ratio of the 
Ta'*!(y, 22)Ta!” cross section, g2, to o; had previously been deduced from 
residual activity and total neutron measurements by Carver, Edge, and 
Lokan (1957), and so the photon absorption cross section o;-++o2 was found. 
This showed a giant resonance of half-width 6.3 Mev, in reasonable agree- 
ment with the calculations of Okimoto (1958), but there was no evidence for 
two peaks. 

The present work gives an account of further measurements which provide 
additional evidence for the existence of two peaks. 


EXPERIMENTAL PROCEDURE 

A stack of tantalum foils, of total thickness 4.98 g cm~?, was irradiated 
in the X-ray beam from the University of Saskatchewan betatron, and the 
photoneutrons were detected by BF; counters; details of the method have 
been described previously (Montalbetti, Katz, and Goldemberg 1953). The 
room in which the detection system is housed was maintained between 19° 
and 20°C throughout the experiments, and the detection efficiency was 
monitored regularly with a Ra~—Be neutron source; the variation during a 
period of several weeks was + 3%. Except near threshold, the counting rate 
was kept near 60 per second, corresponding to a maximum rate at the instant 
of opening the gate of 2500 per second, and observed counts were corrected 
for both counting losses and small variations in the neutron detection efficiency. 

Measurements were made at intervals of 0.2 Mev between 7.6 and 22 Mev, 
the energy scale being monitored regularly by irradiating oxygen samples 
at 17.26 Mev and measuring the induced 2.1-minute 8+ activity. The dose 
was measured with an aluminum wall, transmission ionization chamber, and 
this was subsequently calibrated against a Victoreen thimble in an 8-cm 
lucite cube. The yield curve was normalized by comparing the yields from 
copper and tantalum at 22 Mev and assuming the former to be 2.71 X10° 
neutrons/mole 100 r. 


RESULTS AND DISCUSSION 


In Fig. 1, the photoneutron yield per mole per 100 r, A, is plotted against 
the maximum photon energy E, in the X-ray spectrum; except near the 
vy, n threshold, between 40,000 and 100,000 counts were taken for each point. 
The (y, ”) threshold was found to be 7.63+0.05 Mev, in good agreement 
with earlier measurements (Chidley, Katz, and Kowalski 1958). 

The yield curve was smoothed before being analyzed. It was fitted roughly 
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Fic. 1. Photoneutron yield from tantalum. The full curve is the sine approximation. 


by the function A’ = a{1+sin (8 E,,—6)}, the constants a, 8, and 6 being 
determined by trial and error. The deviations A— A’ of the experimental 
points from the sine curve are shown in Fig. 2; a smooth curve was drawn 
through the points by inspection. It was then assumed that the true yield 
curve is the sum of the sine curve 4’ and the smooth deviation curve. Excluding 
the points near threshold, the r.m.s. percentage deviation of the experimental 
points from the smoothed yield curve was 0.49%, and a rough check showed 
that the distribution was consistent with a normal error curve; this suggests 
that no significant features were removed from the yield curve by the smooth- 
ing process. 
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Fic. 2. Deviations of the experimental points from the sine curve. Where the standard 
deviation is not shown it is smaller than the radius of the circles. 


The yield curve was analyzed by the inverted matrix method of Penfold 
and Leiss (1958) using a bin width of 0.5 Mev; the resulting values of 01+2c2 
are shown in Fig. 3. These values were then corrected for neutron multiplicity 
using (a) the data of Carver and Turchinetz (1959) and (6) calculations based 
on the statistical model described by Blatt and Weisskopf (1952), to give 
the photon absorption cross section o;+¢2. Multiplicities derived from the 
data of Whalin and Hanson (1953) are between those used above. The cross 
section, which is very similar to that found by Fuller and Weiss (1958), has 
two distinct peaks at 12.75+0.25 and 15.25+0.25 Mev, and so, assuming 
that equation (1) is valid, the intrinsic quadrupole moment is 5.7+1.3 barns. 
This is in good agreement with the results of other observers which have 
already been quoted. 

The sudden increase that occurs above 21 Mev is almost certainly spurious; 
the calculated cross section is very sensitive to changes of slope of the yield 
curve in this region and if the observed yields for the last three points are 
reduced by approximately 3% the abnormality is removed. The cross section 
deduced by Fuller and Weiss also shows large fluctuations in this energy 
region. 

Several attempts were made to fit theoretical curves to the experimentally 
measured cross section; either a Lorentz curve, as used by Fuller and Weiss, 
or an optical dispersion curve can be made to fit the data reasonably well. 
The curve shown in Fig. 3 is of the latter type, viz. 
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Fic. 3. The cross section for neutron production ¢;+2¢2 and photon absorption o:+¢:. 
The full curve is compounded of two resonance curves whose parameters are given in the 
text. O o1+202; @ o1+02 (2); + o1+02 (0). 


es qo} eer _ 

on 14+ 1(E—E,)/7T,/2)" * 14+((E-E,)/T. 0’ 
where o; = 317 mb, o2 = 444 mb, E; = 12.5 Mev, Es = 15.5 Mev, Tr; = 2.3 
Mev, and I, = 3.6 Mev. Using these values of the parameters the over-all 
agreement with the experimental points is reasonable, although the separation 
between the peaks of the composite curve is rather too large. The ratio of the 
areas under the peaks of the component resonance curves is o;I°;/o2I2 = 0.46, 
the predicted value being 0.5. Not such a good fit, however, was obtained 
with the standard Breit—Wigner formula in which o; and o2 above are replaced 
by o,(£,/E)? and o2(E2/E)*. 


ACKNOWLEDGMENT 


One of the authors (R. W. P.) is indebted to the National Research Council 
for a Postdoctoral Fellowship. 


REFERENCES 
Biatt, J. M. and WeEtsskorr, V. F. 1952. Theoretical nuclear physics (John Wiley & 
Sons, Inc., New York), p. 379. 
Carver, J. H., EpGe, R. D., and Loxan, K. H. 1957. Proc. Phys. Soc. A, 70, 415. 
Carver, J. H. and TuRcHINETZ, W. 1959. Proc. Phys. Soc. 73, 69. 








814 CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 


CuibLey, B. G., Katz, L., and Kowatsk1, S. 1958. Can. J. Phys. 36, 407. 

Danos, M. 1958. Nuclear Phys. 5, 23. 

FuLLer, E. G. and Weltss, M. S. 1958. National Bureau of Standards Report 5687. 

FuLuLer, E. G. and Haywarp, E. 1958. Phys. Rev. Letters, 1, 465. 

GOLDHABER, M. and TELLER, E. 1948. Phys. Rev. 74, 1046. 

Katz, L. and CuipLey, G. B. 1958. Nuclear reactions at low and medium energies (Aca- 
demy of Sciences of the U.S.S.R.), p. 371. 

PENFOLD, A. S. and LeEIss, J. E. 1958. Analysis of photo cross-sections (Physics Research 
Laboratory, University of Illinois, Chicago). 

MonTAaLBeEttl, R., Katz, L., and GOLDEMBERG, J. 1953. Phys. Rev. 91, 659. 

NaATHANS, R. and HALPERN, J. 1954. Phys. Rev. 93, 437. 

OximoTo, K. 1958. Phys. Rev. 110, 148. 

Spicer, B. M., Tures, H. H., BaGuin, J. E., and ALtuM,F.R. 1958. Australian J. Phys. 
11, 298. 

Wuatin, E. A. and Hanson, A. O. 1953. Phys. Rev. 89, 324. 








THE OCCURRENCE OF (n,3n) REACTIONS IN A 
NUCLEAR REACTOR! 


T. A. Eastwoop AnD J. C. Roy 


ABSTRACT 


The reactions Bi?°9(n, 3n)Bi?°?7 and Au!97(m, 37)Au!® have been observed to 
occur in NRX, a thermal neutron reactor. Competing reactions leading from 
Bi? or Au!®? to the same products have been shown to be negligible. Within 
the uncertainties of the methods used the yields of Bi?°? and Au!® are about 
those expected on the basis of neutron evaporation from the compound nucleus 
formed by the amalgamation of the incoming neutron and the target nucleus. 


INTRODUCTION 


The complete exploration of the neutron spectrum in a nuclear reactor 
using activation detectors requires reactions with thresholds in a region above 
those in current use. The (7, 32) reactions have suitable thresholds (14-25 
Mev) for this purpose but their occurrence in a nuclear reactor has not been 
reported. There is, of course, no reason why such reactions should not occur 
if the excitation of the nucleus is sufficiently high. Marquez (1952) states 
that Cu® is produced during the irradiation of natural copper with neutrons 
of mean energy of about 370 Mev and this is probably the product of a 
(n, 3n) reaction. However, this energy region is much above that encountered 
in nuclear reactors. Watt (1952) has observed neutrons with energies up to 
17.2 Mev in the spectrum of neutrons from U** fission and it is thought that 
the neutron energy spectrum extends to 20 Mev (Hughes 1953). 

The detection of the (7, 3”) reaction is difficult because the product is an 
isotope of the target element and chemical methods cannot be employed to 
separate it from the products of the more probable (7, y) and (”, 2”) reactions. 
Nevertheless, the decay characteristics of the neutron deficient nuclei are now 
sufficiently well known to permit the choice of a few elements in which the 
(n, 3n) reaction might be observed. A consideration of the properties of the 
nuclides indicates that the most likely elements with which to detect the 
(n, 3n) reaction are those in which the product of the (, 32) reaction is 
radioactive and has a reasonably long half-life. The half-life of the product 
should differ by several orders of magnitude from the half-lives of the products 
of the (”, y) or (m, 2”) reactions. Then a long irradiation followed by suitable 
decay period facilitates the detection of the (”, 3”) reaction product. This 
approach requires rigorous chemical purification following the irradiation 
because even trace impurities can yield large amounts of radioactive products 
on long irradiations. 

The nuclear properties of bismuth and gold isotopes shown in Table I are 
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TABLE I 

Nuclear properties of bismuth and gold isotopes involved in neutron-induced reactions 
_  Bjeor —_ Bi208 Bj209 Bjz10 
28 years 3X 104 years Stable 2.6X10° years 
K and L electron K electron capture 100% abundance a—4.94 

capture 
y—0.57, 0.89, 1.06, y—2.615 

1.43, 1.77 

Au!% Au!% Au!97 Au!98 Au! 

180 days 5.6 days Stable 2.70 days 3.15 days 
KandLelectron Electroncapture 100% abundance B-—0.96, 0.29, B-—0.30, 0.25, 

capture 1.37 0.46 
+—0.031, 0.098, 8-—0.27 y—0.4118, 0.68, | y—0.05, 0.158, 

0.13 y—0.331, 0.354 1.09 0.208 





Note: Short-lived isomeric states are omitted. Energies are in Mev. 


such that these elements are suitable targets. The (7, 3”) reaction thresholds 
are 14.3 Mev for Bi?°*(n, 382) Bi? and 13.2 Mev for Au!®"(n, 3”) Au'® (Wapstra 
1958). In addition they are obtainable in pure form and their chemistry is 
relatively well understood. 

EXPERIMENTAL METHODS 
Irradiations 

With the exception of one bismuth irradiation all irradiations were made in 
empty fuel rod positions in the NRX reactor. This isa D,O-moderated, natural- 
uranium-fueled reactor, the physical characteristics of which have been 
described by Hurst (1955). The sample position was in the center of a hexa- 
gonal array of fuel rods and was separated from the fuel rods by 17.3 cm of 
D.O. The thermal neutron flux in the positions used was about 6X10" 
n/cm?-sec. 

Four bismuth irradiations were made. In the first a 160-mg sample of 
bismuth metal was irradiated to an integrated thermal neutron flux of 2.7 
X 102! n/cm*. It was allowed to decay for about 2 years before it was used 
for the measurements reported here. For identification this is called sample A. 
Two other samples of bismuth (called B and C) weighing about 3g each 
were irradiated near each other in the reactor. The integrated thermal neutron 
flux was 2.8X10?°n/cm?. In order to separate thermal neutrons from fast 
neutron effects, one sample (C) was placed inside a cadmium box with 0.064-cm 
thick walls. Measurements with these bismuth samples began 5 months after 
the irradiations. The fourth irradiation was the exception mentioned pre- 
viously. It was made in the center of a hollow natural-uranium cylinder, 3.5 
cm in external diameter and with walls 0.023 cm thick. This ‘‘neutron trans- 
former’’ was placed in an empty fuel rod position in the NRX reactor. The 
sample (D) weighed about 800 mg and was irradiated for 33 days. A period 
of 10 months elapsed before the measurements reported here were made. 
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One irradiation of 250 mg of metallic gold was made. The integrated neutron 
flux was 1.7 X 107° n/cm*. Measurements with this sample were begun 5 months 


after the irradiation. 


Chemical Purification 

(a) Bismuth.—As pointed out earlier, extensive purification of the target is 
required to detect a small amount of Bi?*’ produced. It was found that the 
main radiochemical contaminants were Po, Zn, Hg, Sn, Sb, and Cd. The 
chemical treatment used to remove these contaminants from sample A has 
been described previously (Roy et a/. 1958). This method was not suitable 
for the purification of samples B, C, and D owing to the large amount of 
bismuth involved. The samples were dissolved and the polonium was removed 
as described for sample A. The following preliminary steps were then made 
before the purification treatment described for sample A was applied. Mercury 
carrier as Hg++ was added to the solution and HgeCl, plus some Hg was 
precipitated by the addition of SnCly. The precipitate was removed by centri- 
fugation and discarded. After the addition of cadmium carrier to the solution, 
Bi and Sn in metallic form were precipitated by the addition of finely powdered 
Zn in dilute HCl; Cd remained dissolved under these conditions. The precipi- 
tated metal was treated with hot concentrated HCl, which dissolved the tin, 
and the precipitate was removed by centrifugation. The Bi was dissolved 
upon treatment with HNO;. This procedure was repeated several times. 
Finally, the Bi was further purified by the chemical treatment described for 
sample A. 

(6) Gold.—The irradiated gold sample was dissolved in aqua regia and, 
because a measurement of the amount of 48-day Hg? produced in the 
irradiation was required for reasons presented in the Discussion, a known 
amount of mercury carrier was added. The solution was adjusted to be 3 J 
in HCl and the gold was extracted with ethyl acetate. The organic phase 
was washed several times with 3 N HCl containing platinum carrier. The 
ethyl acetate was removed by evaporation and the residue was dissolved in 
2 NHNOs. Metallic gold was precipitated by the addition of stannous chloride. 
The gold was redissolved in aqua regia and the extraction was repeated. Finally 
the metallic gold was mounted on an aluminum tray for examination of its 


y-ray spectrum. 


Counting Techniques 

The decay of both Bi?®? and Au'® gives rise to well-known y-rays. There- 
fore, the y-ray spectra of the purified Bi and Au were examined with a 
scintillation y-ray spectrometer which consisted of a lin. X 1} in. Nal(TI) 
crystal coupled to an RCA 5819 photomultiplier tube, conventional electronics, 
and a 24-channel pulse height analyzer. The energy scale of the spectrometer 
was determined using y-rays from the decay of Am*4!, Na”, Cs87, Co®, and 
Na*™. The efficiency was determined using sources of Na”, Na’, Au!%8, and 
Cs'87 prepared from solutions standardized by 4x 8-counting techniques. 
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RESULTS 
(a) Bismuth 


The y-ray spectra of the four samples of irradiated, purified bismuth in the 
region of 0.3 to 1.3 Mev are shown in Fig. 1. The two well-known gamma 
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Fic. 1. Portions of the gamma-ray spectra of four samples of irradiated purified bismuth. 
Sample A was irradiated in a predominantly thermal neutron spectrum, integrated thermal 


neutron flux 2.7 X10! n/cm?. 
Sample B was irradiated in a predominantly thermal neutron spectrum, integrated thermal 


neutron flux 2.8 X 10? n/cm?. 
Sample C was irradiated with epi-cadmium neutrons near sample B. 


Sample D was irradiated for 1 month in a “neutron transformer’. 


rays of Bi?®? at 0.569 and 1.063 Mev (Alburger and Sunyar 1955) are clearly 
resolved in all four spectra. The other y-rays of Bi?’ are in low abundance 
and were not detected. It was therefore concluded that Bi?®’ was present in 
the sources. 

The complete y-ray spectra of each of the four samples also show peaks 
at 0.270, 2.1, and 2.6 Mev. These are not associated with the decay of Bi?®’. 
The 2.1- and 2.6-Mev peaks are due to Bi?® as described previously (Roy et 
al. 1958). The y-ray spectrum of sample A in the energy range 1.8 to 3.5 Mev 
was given there. The 270-kev peak is certainly due to bismuth since its peak 
height relative to the peaks of Bi?*? and Bi?®’ does not depend on the degree 
of chemical purification. It seems to be associated with the product of an 
(n, y) reaction since it is much more predominant in samples A and B which 
were irradiated in a thermal neutron flux than in samples C and D which 
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were irradiated with epi-thermal or fast neutrons. This point is illustrated 
by the spectra in Fig. 2. It is suggested that the 270-kev y-ray arises from the 
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Fic. 2. Gamma-ray spectra of irradiated, purified bismuth samples. 

The two samples were irradiated near each other in the reactor. Sample C (open circles) 
was enclosed in cadmium and sample B (full circles) was uncovered. The low abundance of 
the 270-kev y-ray in the spectrum of the covered bismuth indicates that the 270-kev y-ray 
is associated with a thermal neutron induced activity. 


decay of 2.6 10°-year Bi?!®, a hypothesis which is supported by preliminary 
a~y coincidence experiments.* Further work on the source of this y-ray is in 
progress and the results will be reported separately. 

The yield of Bi?’ was estimated from the observed counting rate of the 
1.063-Mev y-ray peak, the efficiency of the y-ray spectrometer (1.6%), the 
known decay scheme of Bi*®’, and the weight of the bismuth sample used for 
counting. The half-life of Bi?°? was taken to be 28 years (Sosniak and Bell 
1959). The results are given in Table II. For reasons to be discussed later it 


TABLE II 


Results of bismuth irradiations 


Observed Bi?*7 Effective (, 3n) 


Integrated Weight of 


Irradiation thermal neutron purified disintegration reaction 
position in flux, sample, rate, cross section, 
Sample NRX n/cm? g dis/sec cm 
A Empty fuel rod 2.7X107! 0.054 73 2.21073 
B- Empty fuel rod 2.8X 107° 0.308 47 2.4X1073! 
C Empty fuel rod 2.8X 102° 0.654 94 2.2X10"3! 
(Cd-covered sample) 
D ‘Neutron trans- 1.9 1020* 0.150 50 7.8X1073 


former” 


*Estimated integrated thermal neutron flux in this position in the absence of the ‘“‘neutron transformer’’. 


*NOTE ADDED IN PROOF.— This hypothesis is substantiated by the work of S. V. Golenetskii, 
L. I. Rusinov, and Iu. I. Filimonov, Zhur. Eksptl. i Teoret. Fiz. 35, 1313 (1958); JETP, 8, 
917 (1959). 
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is convenient to deduce an effective cross section for the Bi?°(n, 3m) Bi?®7 
reaction which, when multiplied by the thermal neutron flux, gives the observed 
reaction rate. The effective cross sections so obtained are also given in Table II. 


(b) Gold 
The y-ray spectrum of the purified gold sample is shown in Fig. 3. The 
y-ray spectrum was measured 5 months after the end of the irradiation so 
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Fic. 3. X- and y-ray spectrum of irradiated, purified gold. 


that the Au!®’, Au'®®, and Au'®6 activity had essentially all decayed. Cressman 
and Wilkinson (1958) have found that Au'® decays entirely by electron cap- 
ture and that the electromagnetic spectrum consists of a strong K X-ray 
peak at 67 kev, less intense peaks at 31 and 98 kev, and a very low intensity 
peak at 129 kev. All but the lowest intensity peak are present in the spectrum 
in Fig. 3. No y-rays above 98 kev were observed in the y-ray spectrum and 
no electrons with energy sufficient to penetrate 3 mg/cm? of aluminum were 
detected from the sample of gold, thus attesting to the radiochemical purity 
of the sample. The activity, observed for a period of 200 days, decayed with 
a half-life of about 190 days, in agreement with the average of published 
half-lives of 183 days for Au!® (Strominger 1958). It was therefore concluded 
that Au'® was present in the irradiated gold samples. 

It is difficult to estimate the yield of Au'!® from this irradiation because 
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the amount of internal conversion of the y-rays is not accurately known. 
However, approximately 1000 disintegrations per second of Au'® were ob- 
served in 230 mg of gold corresponding to an effective cross section of 
3<10-* em?. 


DISCUSSION 


The experimental results presented in the preceding section show that Bi?°7 
and Au!® are present in samples of bismuth and gold irradiated in the NRX 
reactor. The isotopes might be produced by a number of mechanisms but it 
will be shown that the most probable is the (m, 3) reaction. 


(a) Bismuth 

Let us first consider the possibility of two successive (n, 2”) reactions. The 
disintegration rate of the product of two successive (”, 2”) reactions is given 
by a relationship of the type: 


dNp N uncle oT. Bo? aes Nas 
(1) dt 0.698(T4—T a) [Tae *4)—Ta(1—e*?)] 
where the subscripts S, A, and B refer to the target nuclide, and the product 
of the first- and second-order neutron capture reactions respectively; Ng is 
the number of atoms of target nuclide; o is the effective cross section of the 
individual (n, 27) reactions; ¢@ is the neutron flux; 7 is the half-life; \ is the 
disintegration constant; and ¢ is the duration of the irradiation. 

All factors in equation | are known except the cross-section values. Ideally 
these would be measured but there are technical difficulties. Therefore we 
have made an estimate of the appropriate cross sections based on the cross 
sections observed for other (m, 2”) reactions under similar irradiation con- 
ditions. The results of these auxiliary measurements are given in Table III. 





TABLE III 
Relation between the effective > (n, 2n) cI cross section and d Bott 

Target Effective cross section,* 

nuclide | millibarns a / A?! Fett, Mev 
Th? 2. 7 +0.7 7. 2107? 7.540.2+ 
Pb? 0.094+0.02 2.8X10-3 10.2+0.5 
Tr 0.13 £999 3.6107 10.840.6 
[227 0.052+0.01 2.1X10-* 11.0+0.7 
As* 0.009 +0 .002 5.110 12.5+0.2 
ys 0.005+0 .002 2. 510-4 14.340.4 
Fe™ A. O+0 2) X10" ‘tf 7.2X10~ 6 16.7+0. 3 








*The ‘athens cross section is used only because it is convenient. It has no idepttiaees except that it gives the 
observed reaction rate in the particular neutron spectrum used in this investigation when multiplied by the thermal 
neutron flux. 

tThe errors given in this column are due to the uncertainties in the masses of the nuclei involved (Duckworth 
1958) plus an uncertainty of about +0.2 Mev in the estimation of Eegf. 

{Betts and Dahlinger (1957). 


It should be emphasized that the tabulations are effective cross sections. The 
only significance of the effective cross section is that when multiplied by the 
thermal neutron flux it gives the observed reaction rate. Effective cross 
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sections are used as a matter of convenience only. Clearly they apply to this 
particular neutron spectrum only and cannot be used directly to give the rate 
of reaction in a different neutron spectrum. 

In order to form the best estimate of the bismuth cross sections required 
it is desirable to interrelate the cross sections observed for the other (”, 2n) 
reaction on some theoretical basis. At bombarding energies below about 
20 Mev, according to current theories, fast-neutron-induced reactions in 
which neutrons are emitted proceed predominantly by the amalgamation of 
the incident neutron with the target nucleus to form a compound nucleus 
with all the available energy as excitation energy. This excitation energy is 
then carried off by one or more nucleons (‘‘evaporation”’). For heavy elements 
charged particle emission is hindered by the coulomb barrier and neutron 
emission predominates. The general feature of (m, 27) reactions in this energy 
range are reasonably well understood (see for example Blatt and Weisskopf 
1952). The cross section o is related to the excitation energy E by the relation 


(2) o = o,(1—(1+y)e] 
where 
_ £-Fr 
i 7 ’ 


Ey is the threshold energy, 7° is the nuclear temperature, and o, is the cross 
section for the formation of the compound nucleus. 

Equation 2 would give the desired relation for monoenergetic neutrons, but, 
of course, the reactions occurred in a spectrum of neutron energies from the 
reactor. It is therefore convenient to define an effective neutron energy, Ee, 
following Hughes (1953). The effective energy is chosen in such a way that 
the correct reaction rate is obtained on the assumption that neutrons below 
Ee produce no reaction while all those above it are equally effective in 
producing the reaction. Then 


(3) o = o(1—(l+y) e J ic vee |/| [owceae | 
where 
iz v(Eyae | / | fovea | 


is the fraction of neutrons with energy greater than Eg. 

For the present purpose this relation can be simplified by assuming (1) o, 
is equal to the geometric cross sections m 2.3 X 10~°* 4? cm’, (2) the factor 
1—(1+y)e~” is constant (in fact it varies from 0.30 to 0.44 for the reactions 
involved here), and (3) the fraction 


lf. veya | / | [vcnae| 


is proportional to e~** (a very good approximation for the fission neutron 
spectrum at energies above 6 Mev). Thus equation 3 reduces to 
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(4) o = kA ¢ Mat 
or 
(5) qin =k Pa 


The validity of equation 5 has been tested by plotting log ¢/ A?” against 
Een in Fig. 4. The experimental results are shown by the points with their 
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Fic. 4. The relationship between the observed effective (”, 2m) cross sections and the 
effective neutron energy. The effective cross section has no significance except that it gives 
the observed reaction rate in the particular neutron spectrum used in this investigation when 
multiplied by the thermal neutron flux. The points are the experimental results with an estimate 
of the errors. The line has the slope expected from theory. 


associated errors. The line in Fig. 4 has the slope predicted by evaporation 
theory. Since the points lie on a line of approximately the predicted slope, 
estimates of the cross section of other (m, 2m”) reactions can be made with 
some confidence. On this basis the maximum cross sections for the bismuth 


and gold reactions were deduced and are given in Table IV. 

The effective cross sections in Table IV were used to predict the yield of 
Bi2°? from two successive (7, 27) reactions starting with Bi?®. It was found 
to be 2X10-? disintegrations per second for sample B. This is about 1/3000 of 
the observed yield and shows that successive (”, 27) reactions are negligible. 

Inspection of the properties of nuclides in the neighborhood of Bi?” shows that 
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TABLE IV 


Expected effective (, 2”) cross sections 





Maximum expected 





Target effective cross section, 
nuclide Eetr, Mev millibarns 

Bi209 8.8+40.2 1.3 

Bi28 8.0+0.2 2.9 

Au!97 9.5+1.1 0.89 

Au!%6 7.0+1.1 19.0 

Hg! 9.3+1.1 1.1 





the only other means, apart from the (”, 3”) reaction, by which Bi?’ could be 
formed during a neutron irradiation of bismuth is by the series of reactions 
Bi2°"(n, y) Bi2!® - B-decay _, po210(q, y) Po?! __B-decay . ayo __a-decay _, i207, 


All of these reactions are known to occur except the 6-decay of Po*"'. Two 
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Fic. 5. Decay schemes of nuclei in the Bi?°7 region (Strominger et al. 1958). 
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isomeric states of Po*!! have been observed and both decay by the emission 
of an a-particle. No 6-decay has been found. It is clear from Fig. 5 that B-decay 
from the 25-second state of Po*" to At?" is energetically possible although 
highly forbidden if the spin assignments are correct. Therefore this mechanism 
for Bi?®? formation cannot be neglected. The proposed chain of reactions 
involves two (n, y) reactions which would be expected to be more probable 
in a thermal neutron flux than in a fast neutron flux. The yield of Bi?®? would 
therefore be expected to be smaller in the cadmium-covered sample C than 
in the uncovered sample B. This is contrary to the experimental observations 
and thus this possibility is eliminated. 

It is therefore concluded that Bi?” found in neutron-irradiated bismuth 
was produced by the Bi?°*(n, 32)Bi? reaction. 


(b) Gold 

In order to estimate the importance of two successive (m, 2”) reactions in 
the production of Au'® from Au!%’ the cross sections of the reactions involved 
were estimated by the method outlined for bismuth and the results are given 
in Table IV. Using these cross sections the expected yield of Au'® by this 
route was found to be 2 disintegrations per second or 1/500 of the observed 
yield. Here again it was concluded that successive (m, 2”) reactions make an 
insignificant contribution to the observed yield of Au!®. 

A possible source of Au'® in the present experiments is from mercury 
impurity in the gold sample which was irradiated. The Hg!**(m, 2”) reaction 
leads to Hg!® which in turn decays to Au'® by electron capture. The amount 
of mercury in the gold sample was determined by activation analysis. A known 
amount of mercury carrier was added during the gold purification procedure 
outlined in a previous section. From the chemical yield and the disintegration 
rate of the recovered 48-day Hg?®, determined by counting the 0.28-Mev 
y-ray, the amount of mercury was found to be 7X10~-* g. If the Hg!*°(m, 2m) 
Hg'® cross section is taken to be 1.1 millibarns (Table IV) the expected 
disintegration rate of Au'® from this source is 2X10~* per second. Clearly 
this reaction is an unimportant source of Au'®. 

On the basis of these results it is concluded that the reaction Au!*?(m, 37) 
Au'® with a threshold energy of 13.2 Mev has occurred in reactor irradiations. 

The emission of three neutrons from the compound nucleus at low excitation 
energy is predominantly an evaporation process and it should be possible to 
predict the cross section on the same basis as the (m, 2m) cross sections were 
predicted. When this is done the observed cross section for Bi?°*(n, 32) Bi?®? is 
found to be one-quarter of that expected. This is probably a satisfactory 
agreement in view of the approximations involved. The agreement for the 
Au!97(n, 3n)Au'® reaction cross section is not as satisfactory; the observed 
cross section is one-sixth of that calculated. However, as pointed out earlier 
there is insufficient information about the decay of Au'® to permit an accurate 
measurement of the effective cross section. 

The occurrence of (n, 3”) reactions indicates that neutrons of high energy 
exist in the sample irradiation positions of a well-moderated thermal neutron 
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reactor such as NRX and it would be of interest to search for neutrons of 
still higher energies using suitable threshold detectors. One promising reaction 
is the C(m,2n)C" reaction which has a threshold energy of 20.3 Mev 
(Ee = 24.5 Mev). Extrapolation of the relationship in Fig. 4 indicates that 
the yield of C'! should be sufficient to permit its detection by the 510-kev 
annihilation radiation. The results of such an investigation are reported in 
the paper which follows. 
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EVIDENCE FOR NEUTRONS WITH AN ENERGY OF 
20.3 MEV IN THE NRX REACTOR! 


j}--€. Roy 


ABSTRACT 


Experimental evidence is presented on the occurrence of the C!2(n, 2n)C¥ 
reaction in the NRX reactor. The probability of occurrence of the competing 
processes C!?( y,2)C" and C(p, pn)C" is discussed. A comparison is made 
of the observ A reaction rate with the one predicted by the formula re presenting 
the energy distribution of the prompt neutrons from the thermal fission of 
U5, Since the threshold of the C"(m, 2”)C" reaction is 20.3 Mev, observation 
of the reaction in the NRX reactor means that there are neutrons with an 
energy of at least 20.3 Mev produced in fission. 


INTRODUCTION 


It seems to be accepted that the fission neutron spectrum extends to about 
20 Mev (Hughes 1957). However, experimentally the fission neutron spectrum 
has been determined only up to 17.2 Mev (Watt 1952). Watt measured the 
ranges of knock-on protons by a telescope of proportional counters. The 
preceding paper on the observation of (m, 3”) reaction in the NRX reactor 
and the comparative studies of (m, 2m) and (n, 3n) effective cross sections 
supports the belief of Hughes that the fission neutron spectrum extends to 
20 Mev or more. Eastwood and Roy (1959) also suggested that the studies 
of some high-energy threshold reactions of the types (, 2”) and (n, 3”) could 
be used to examine the fission neutron spectrum in the 20-Mev range. One 
such reaction is C?(m, 2)C"™ which has a Q value of 18.72+0.02 Mev (Wapstra 
1955) and a threshold of 20.3 Mev. Experimental evidence on the observation 
of the C”(n, 2”)C" reaction in the NRX reactor is reported here. 


EXPERIMENTAL 

Samples of 0.2 g of anthracene were irradiated in a vacant fuel rod position 
near the center of the NRX reactor for 30 minutes at a thermal neutron flux 
of about 6.0 X10" n/cm?-sec. Immediately after the irradiation the anthracene 
was transformed to carbon duals and water by combustion in a Parr oxygen 
bomb. After the combustion the gases were passed through concentrated 
sulphuric acid and ‘‘Anhydrone’’; finally the carbon dioxide was absorbed 

1 ‘‘Ascarite’”’. The ‘‘Ascarite’’ was pressed onto an aluminum tray (3 cm in 
diameter) to form a cake about } cm thick. This cake constituted the source 
which was counted as such. The time necessary to complete these operations 
was about 20 minutes. 

Carbon-11 is a 0.96-Mev positron emitter with a half-life of 20.5 minutes 
(Strominger et al. 1958). Its radioactivity therefore can be detected easily by 
measuring its annihilation radiation or its positrons. In this work, the gamma 
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spectrum was examined with a Nal (TI) (1}in.X1 in.) scintillation spectro- 
meter. The positrons were measured with a low background beta-anthracene 
counter. 

RESULTS AND DISCUSSION 


In these experiments the only gamma radiation found in the source had an 
energy of 0.511 Mev as determined by comparison with the annihilation 
radiation from positrons from Na”. The pulse height spectrum is shown on 
the right side of Fig. 1. The gamma ray decayed with a half-life of 21 minutes 
as seen in Fig. 1. It was therefore ascribed to the annihilation of the 0.96-Mev 
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Fic. 1. Line 1 is the decay curve of the positrons; line 2 is the decay curve of the gamma 


radiations and the points represent 10 times the observed counting rate. The gamma-ray 
spectrum of C" is shown on the right side of the figure. 
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positron of C". Similarly the beta radiations detected on the anthracene 
crystal decayed with a half-life of 20 minutes as shown in Fig. 1. They were 
attributed to the positrons from C'!. The radiochemical purity of the carbon 
dioxide obtained by the combustion of anthracene was excellent; no activity 
other than C" could be detected either with the gamma-scintillation spectro- 
meter or the beta detector. The measured half-lives of the annihilation and 
positron radiations were obtained by subtracting from the observed counting 
rate the background of the detectors. For the anthracene counter, the back- 
ground was 10.0+0.3 counts/minute; for the gamma-scintillation spectro- 
meter it was 15 counts/minute under the energy range covering the annihilation 
radiation. The results of the measurements are summarized in Table I. The 


TABLE | 


Measurements leading to the identification of C" 


Half-life, Method of 


Experiment Method of Radiations 

No. detection detected minutes irradiations 
1 Anthracene Bt 19.0 Nylon capsule 
2 Anthracene gt 20.7 Nylon capsule 
3 Nal (TI) Annihilation 21.0 Nylon capsule 
4 Nal (TI) Annihilation 21.5 Nylon capsule 
5 Nal (TI) Annihilation 21.0 Iron capsule 
6 Nal (TI) Annihilation 22.0 Iron capsule 





half-lives listed in the last column have been calculated from the experimental 
points by the method of least squares, each point being weighted according 
to the equations developed by Behrens (1951). In terms of Behrens’ definitions 
the case under study was treated as a single exponential decay plus a back- 
ground of known mean strength. In each determination the calculated half-life 
is in agreement with the accepted value of 20.5 minutes for that of C!. 

In principle the nuclear reactions capable of yielding the C" activity 
produced under the experimental conditions of the irradiations are the 
following: the C(y,2)C" reaction which has a threshold of 18.73 Mev 
(Bromley and Rutledge 1958), the C(p, d) or (p, pn)C"™ reactions caused 
by recoil protons and which have thresholds at 18.0 and 20.3 Mev respectively 
(Bromley and Rutledge 1958), and finally, the C(n, 2m)C" reaction which 
has been discussed in the introduction. The probability of occurrence of each 
process will be now examined. 

The maximum energy of the gamma rays detected in fission is about 
8 Mev (Fraser 1956). The gamma rays arising from (m, y) reactions in the 
structural materials of the reactor have energies of no more than about 
8 Mev while those arising from the decay of fission products have energies 
ranging from a few kev to a few Mev. Thus the gamma rays present in the 
reactor cannot account for the (y, ”) process since its threshold is 18.73 Mev. 

Cowan and O’Brien (1955) conclude that the (p, pn) reaction does not 
yield a significant amount of C' when anthracene is used as a fast-neutron 
detector based on the C!(n, 2n)C"™ reaction. However, their conclusion is 
open to criticism because it was reached on the assumption that the cross 
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section of the (p, pm) reaction is one-sixth that of the (m, 2”) reaction. In 
fact the measured cross sections of the two reactions (p, d) and (p, pn) (Aamodt 
et al. 1952) are about seven times greater than the (, 27) cross section (Brolley 
et al. 1952) from threshold up to 40 Mev. 

In order to determine the relative importance of the (p, d) or (p, pn) and 
(n, 2n) reactions in the experiments described here, the following method 
was used. In one group of irradiations (experiments 1 to 4, Table I) the 
anthracene samples were irradiated in a nylon capsule so that the concen- 
tration of protons per square centimeter surrounding the anthracene was 10 
times greater than in the anthracene sample itself; in another group (experi- 
ments 5 and 6, Table I) the anthracene samples were irradiated in iron 
capsules. In both groups of irradiations the reaction rates were essentially 
the same, thus showing that the C' was produced mainly by a (, 27) reaction. 
Moreover this experiment demonstrated that the (p, pm) or the (p, d) pro- 
cesses do not cause appreciable error in the measurement of fast-neutron 
fluxes using the C#(n, 2”)C"™ reaction in anthracene. 

The measured reaction rate, R, for the reaction C?(n, 2”)C'! may be 
expressed as the product of an effective cross section, ¢, and the conventional 
neutron flux, ¢, as measured with a cobalt monitor 


(1) R = a¢. 


R was found to be equal to (6+2)X10~! disintegrations per second per 
atom of carbon. This value divided by the neutron flux, 6X10! n/cm?-sec, 
gives an effective cross section of 10-* cm’. 

It would be interesting to compare this observed reaction rate with the 
one predicted by the usual formula representing the energy distribution of 
the prompt neutrons from the thermal fission of U?* (Watt 1952) 


(2) f(E) = sinhv/2E exp (—E) 


where EF is the neutron energy in Mev. The integrated neutron spectrum, 
F(E), that is the fraction of neutrons having an energy greater than E is 
given by 


(3) F(E) = | fre | J. | [He yar |. 


In an empty fuel rod position of the NRX reactor the ratio fast flux/thermal 
flux is about 0.07. This value has been obtained from the ratio of the Ni®’ 
(n, p)Co®® cross sections measured respectively in an empty fuel rod and in 
a uranium receptacle slug (Hughes 1953; Mellish et a/. 1958). The fast- 
neutron flux is therefore 4X10" n/cm?-sec. Thus the product F(£)X4X10"” 
n/cm?-sec gives the neutron flux, ¢g, having an energy greater than E. The 
cross section for the reaction C?(m, 2”)C" is known from threshold up to 
340 Mev (Brolley et al. 1952; Hughes 1958). Now a reaction rate, Ri, can 
be evaluated by graphical integration from the values of ¢g and o as follows: 


20 20 
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where ¢ is the mean cross section over the energy interval AE. 

The summation gives a value for R; equal to 2.5 10~*! disintegrations per 
second per atom of carbon which is about two times lower than R, the observed 
reaction rate. Considering the uncertainties associated with the values chosen 
to calculate R; and the error in the determination of R, the agreement between 
the observed and calculated reaction rates is good. It means that the relation, 
given in equation 2, can still represent at least qualitatively the energy 
distribution of the prompt neutrons from the thermal fission of U®® in the 
range 17.3-25 Mev. 

Since the threshold of the C"(m, 2”)C™ reaction is 20.3 Mev, observation 
of the reaction in the NRX reactor means that there are neutrons with an 
energy of at least 20.3 Mev produced in the fission of fissile material present 
in the NRX reactor. Although the fuel in the NRX reactor is U*®, a significant 
fraction of fissions occur in the plutonium in irradiated fuel rods and in the 
Pu-Al booster rods; also about 6% of the fissions are due to fast-neutron 
fissions of U**. At this point it is interesting to consider what is the upper 
limit of the fission neutron spectrum. At the present time this question cannot 
be answered without further information on the fission process. However, one 
probably could still extend the limit to a higher energy than 20.3 Mev by 
the studies of fast-neutron reaction of the type (”, 3”) having thresholds in 
the 20-25 Mev range (Eastwood and Roy 1959). 
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Kr 86 AND Hg 198 WAVELENGTH STANDARDS! 


K. M. BarirD AND D. S. SMITH 


ABSTRACT 


New values are given for the vacuum wavelengths of six lines of Kr 86 and 
five lines of Hg 198 in the visible spectrum. The wavelengths were referred to 
the Kr 86 primary standard line, 6057.80211, by a series of very careful measure- 
ments over a period of several years. The accuracy with which wavelengths can 
be compared by present methods is discussed at some length. 


INTRODUCTION 


In order to present new values of wavelength standards in proper perspective 
it is important to consider rather carefully the accuracy to which such standards 
can be determined at present. Recent work in connection with the redefinition 
of the international metre in terms of optical wavelengths has thrown con- 
siderable light on this question and indicates that the accuracy is much less 
than has been generally realized. A few years ago the impression was widely 
held that an accuracy of .0001 A in the visible region, i.e. about two parts 
in 10%, was not out of the ordinary for wavelength values (Candler 1951; 
Tolansky 1955). In reality errors 10 times this amount existed in published 
values of secondary wavelength standards and even in recent determinations 
of the wavelengths of new isotope sources, which were expected to be much 
more precise than the earlier standards, it has been difficult to find agreement 
to +.0001 A. This is illustrated in Table I, which lists values for the orange 


TABLE I 


Recent published values of mace standards 


No. Laboratory Kr 86 ‘te 198 References 

1 N.B.S. (U.S.A.) 5462 .2707* Meggers 1947 

2 N.P.L. (England) 2707" Barrell 1947 

3. B.I.P.M. (France) 2107" Perard 1950 

4 N.B.S. (U.S.A.) 2707* Meggers 1950 

5 N.P.L. (England) 2702* Barrell 1951 

6 N.B.S. (U S.A. ) . 2702, Kessler 1955 

7 N.P.L. (England) 6057 . 8014 2702 Comité intern. pds. mes. 1956 
8 B.I.P.M. (France) 8021 2707* | 

9 P.T.B. (Germany) 8023 2707*t | 

10 N.P.L. (England) 8021 2707 | 

11) = C.ILI.W.M. (Japan) 8021 2708* Comité intern. pds. mes. 1958a 
12 N.R.C. (Canada) 8021 2704 [ 

13 ON. R C. (Canada) .2706t 

14. N.S.L. (Australia) 2707* 

15 B.I.P.M.( France) .2708t Terrien 1958e¢ 


*Uncorrected for effect of argon pressure (+ .0001 to 40,0003 A), 
tReferred directly to Kr 86 = 6057.8021, others referred to Cd = 6440.2491, 
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line of Kr 86 and the green line of Hg 198 which have been obtained at various 
standards laboratories during the last 10 years. All the values of this table 
have either been measured in vacuum or reduced to vacuum by the Edlen 
formula; in the former case where the Cd red line was used as a reference its 
value in vacuum was obtained by the same formula. The Kr 86 wavelengths 
were all measured for a hot cathode lamp of the form described by Engelhard 
(1951) in all cases at the triple point of Nz except No. 11, which was at the 
boiling point of No. Otherwise the lamps were used in various ways—viewed 
from the side, cathode end, or anode end, and at currents, whose exact value 
was not always specified, up to 6 ma/mm?. 

The Hg 198 wavelengths were measured for electrodeless lamps of the type 
due to Meggers (1950a) having argon as the carrier gas at an original pressure 
(at the time of filling) of from } to 5mm Torr. The lamps were cooled by 
water at from about 0° C to 15° C without always accurate specification of 
the water temperature. 

The table illustrates well the present uncertainty in practice in the deter- 
mination to the fourth decimal in angstroms of even very good wavelength 


standards. 


REASONS FOR SPREAD IN MEASURED VALUES OF WAVELENGTHS 
Although the quite large uncertainty in the fourth decimal of angstroms 
may be disappointing, it can quite reasonably be accounted for as indicated 


in Table II. 


TABLE II 


Various uncertainties contributing to scatter in wavelength comparisons 


Errors of observation 


Uncertainty due to (systematic and random- 
conditions of pointing error plus 
Source excitation instrumental error) 
Cd $1.5 (+ <1.0) +15 x10 A 
Hg 198 +3 (+0.2) +1.1 x10 A 
0 
Kr 86 +15 (+0. 2) +0.8 x1o-# A 
0.2 


The magnitude of the uncertainties given in Table Il which are ascribed 
to conditions of excitation are now known from actual measurements of the 
effects of current, temperature, carrier gas pressure, etc., on the emitted 
wavelengths. These uncertainties, which apply to many of the determinations 
of Table I, could be reduced by suitable control of the conditions of excitation 
to the values shown in parentheses. 

One important cause of shift of the Hg 198 green line results from foreign 
gas broadening which has been found for argon to be .00006 A to the red 
.0003 A for 5 mm pressure (Baird and Smith 


per millimetre pressure, i.e. 
1957). After some use it is likely that the pressure in the lamp becomes reduced 
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through cleanup but, on the other hand, electrodeless lamps are sometimes 
reactivated by heating and this procedure may well raise the internal pressure 
to above the original pressure. Another important cause of shift results from 
self-absorption if the lamp is not sufficiently cooled; the self-absorption is 
about 40% at 18°C and will affect the fourth decimal place in angstroms 
(Terrien 1958a). The effect of the nature and intensity of excitation on the 
Hg 198 lines has not been measured exactly; however, tests at this laboratory 
have not indicated any significant effect provided the power is not so great 
that the lamp becomes unduly heated. For most of the values for the mercury 
line given in Table I one cannot say exactly what total correction should be 
applied although the fact that the likely shifts are known to be toward the 
red would make one tend to favor the shorter values. It can be said, however, 
that if the Hg 198 lamp is excited at low power (so as not to have high- 
temperature gradients through the walls of the lanip) with cooling water at 
less than 10° C and having the argon carrier gas pressure less than 0.3 mm 
then the wavelength of the green line can be assumed with a reasonable 
degree of certainty to be reproduced to about +.00002 A. The best form 
of lamp is that with a large ballast volume operated at about 300 megacycles. 
Earlier observations that the line width was reduced by the use of very high 
frequency excitation have not been substantiated and are likely due to the 
fact that the large power used caused a doublet effect due to reversal which 
caused an apparent sharpening of fringes at certain path differences. 

The shifts which affect the Kr 86 orange line result from interatomic Stark 
effect and from Doppler effect due to ion drift in the cases where the lamp 
is viewed end on. The magnitudes of these effects were measured by Engelhard 
(1958) and by the present authors (see below). The uncertainties arising from 
these causes in the determinations listed in Table I would account for a total 
scatter of about .00015 A as indicated in Table I]. On the other hand, if 
the Késters-Engelhard Kr 86 lamp is operated at the triple point of No, at 
a current density of less than 4.0 ma/mm?’, and has the appropriate corrections 
applied for current and direction of viewing (where this applies), then the 
wavelengths of the orange line will also be reproduced to about + .00002 A. 

Although the large scatter of earlier measurements of wavelengths has 
sometimes been largely attributed to the Cd red line, which was used 
for reference, independent measurements by Baird (1958) and by Masui, 
Terrien, and Hamon (1958) showed no reason to expect a variation of much 
over + .0001 A as long as the conditions of excitation were reasonably close 
to any of the various specifications which have been officially adopted; this 
tends to be confirmed by the measurements of the Kr 86 orange line relative 
to the Cd line which are shown in Table I. More rigorous specifications 
could reasonably be applied which would reduce the uncertainty due to con- 
ditions of excitation to less than +.0001 A. Thus the Cd red line compares 
unfavorably with the other two lines discussed above, as they have been 
used in the past, only insofar as its width is two to two and a half times as 
great and consequently individual pointings are less precise by about this 
factor. The Michelson lamp is very inconvenient and no longer useful but 
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the electrodeless Cd 114 lamp is simple to make and to operate and would 
make an excellent source for secondary standards. 

The estimated error limits due to observation cannot be stated other than 
by an intelligent guess. However, a quantitative consideration of the many 
possible sources of error shows that the estimates indicated are certainly not 
too large and they are in fact consistent with the accuracy generally quoted 
by the authors of the values of Table I. 

As indicated in the table these errors can be ascribed to statistical fluctua- 
tions in the average of many pointings, systematic errors of pointing and 
errors due to the interferometer and atmosphere. They have been discussed 
in detail elsewhere (Terrien 19580; Baird 1958a) and will not be repeated 
here, beyond saying, as an example, that those who have had a great deal 
of experience in pointing on lines by the use of the most precise conventional 
methods and equipment have found that the loss of accuracy from this source 
alone is at least 1% of the line width (i.e. about .0001 A), however repro- 
ducible the pointing may seem under some given conditions. New methods 
of pointing on fringes are now being developed which appear to make possible 
a great improvement in precision. On the other hand there is no real assurance 
that a corresponding elimination of systematic errors of pointing has been 
attained as yet. 

It is quite clear that any reasonable convention for combining the errors 
of Table I] will easily account for the scatter in values, shown in Table I, for 
the intercomparison of any two of the wavelength standards. Indeed the good 
agreement of the last six values for the Kr 86 orange line appears even better 
than one might expect. 

In summing up the foregoing, one conclusion is that an uncertainty of 
greater than +.0001 A must be ascribed even to the best determinations of 
wavelength values up to the present. It may be argued that some of the 
previous determinations are better than others but there is really no positive 
basis on which to assess them. For example No. 6 of Table I is called a pre- 
liminary value although it is known that it was made in vacuum, that the 
lamp condition was well known and controlled, and that the result was 
quoted to + .000008 A; the main source of uncertainty is probably in the 
Cd red line to which it was referred but this is very unlikely to have contributed 
an error greater than .0002 A as mentioned above. On the other hand some 
of the later values in the table which are presumed to be definitive are known 
to have been made by the use of old lamps without application of accurate 
corrections; or they were made by new untried methods; or the lamp was at 
too high a temperature; and so on. 

A second conclusion is that as far as the sources themselves are concerned 
it ought to be possible to obtain approximate agreement in the fifth place 
of decimals in the case of lines emitted by Hg 198 and Kr 86, provided con- 
ditions of excitation are properly controlled. (Although only the green and 
orange lines have been considered, the conclusions as to reproducibility apply 
to many other lines from these sources.) Errors of observation can also be 
reduced by care even without a change from conventional methods, e.g. by 
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measurement always in vacuum or by reference to vacuum and by adequate 
statistical evidence in the observations. 

The wavelength determinations described below were made with the aim 
of obtaining the best possible values by the use of methods which are quite 
standard at present. 

The isotope lamps were made in the laboratory so that their characteristics 
were accurately known and perturbing effects due to the excitation conditions 
were taken into account. Measurements were made by the use of three 
independent methods over a long period of time, and by different observers 
with completely re-set-up equipment, in order to ensure good statistical 
reduction of systematic errors that may have been associated with any one 
method, set of interferometer plates, etc. 

It may be remarked that new techniques have been under development 
during the past few years in this laboratory and elsewhere which promise to 
yield an improvement in accuracy by a factor of 10 to 100. The attainable 
precision for the new methods, which is limited-by quantum effects, has been 
discussed theoretically by Hanes (1959). 

However, before wavelengths can be reported to such accuracy a great deal 
of work remains to be done, first to give assurance that systematic errors of 
observation are eliminated to an extent consistent with the improved precision 
of pointing and secondly to determine and control to the same degree of 
certainty the small shifts of wavelength due to variations in the conditions 


of excitation. 


EXPERIMENTAL PROCEDURE 
The vacuum wavelengths given in Table III for Kr 86 and Hg 198 were 
measured by reference to the Kr 86 line corresponding to the Paschen desig- 


TABLE III 
Vacuum wavelengths referred to Kr 86 6057.8021 A 

6458 071942 5792. 268343 
5651.1285 1 5771.1982 3 
4503.6159 1 5462.2705 1 
1464.9414 1 4359 .5621 2 
4455.1664 2 4047.7144 2 
4377 3500 1 


nated transition 2p, )-5d5; the wavelength of this line was taken as 6057.80211 
«*10~'" metres following the recommendation of the Comité International des 
Poids et Mesures (1958) and Commission 14 of the International Astronomical 


Union (1958). 

The measurements were made principally by the use of a Fabry—Perot 
interferometer which was mounted in a vacuum chamber whose temperature 
was controlled by a water bath. Photographs of the Haidinger fringes were 
made in the usual way by the use of a Hilger D260 prism spectrograph; they 
were made in symmetrical sets, such as Hg 198-Kr 86—Kr 86-Hg 198, on any 
given plate. The control of temperature and the vacuum were sufficient to 
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assure the introduction of no errors not completely negligible in the fourth 
decimal of angstroms. The usual precautions were taken regarding the 
parallelism of the optical flats, proper flooding of the apertures, the use of a 
lens system with negligible distortion and chromatic aberration, etc. The 
etalon plates were coated with aluminum of such density that the transmission 
was about 15%. The etalon length for most exposures was 67. 5 mm. Photo- 
graphs were also made with 150-mm and 5-mm etalons respectively; the 
results with the former showed no evidence of asymmetry in the lines; the 
5-mm etalon was used to eliminate the effect of dispersion of phase change 
on reflection. 

Measurements of each of the photographic plates were made by at least 
two observers who used visual setting with an accurate Zeiss Abbé comparator; 
some measurements were also made by means of a photoelectric comparator. 

These measurements were made over a period of several years with different 
setups and different etalon plates. Most of the values used for Table III are 
the average of 50 observations on 18 plates. 

In order to check for systematic errors some of the lines, listed in Table IV, 
were measured visually by the use of two other interferometers which were 


TABLE IV 


Comparison of wavelengths obtained by two different 
methods 





d (A), 
whole part Fabry-Perot Michelson 
O7195 O718; 
¢ 1285 1284, 
5792 2684 2682 
5771 1982; 1980 
5462 2704; 27055 


basically of the Michelson type. One of these (Baird 19580) made use of the 
half-shade technique and by using only highly collimated light with parallel 
mirrors and the same optical surface at the two end points of the reference 
distance; it avoids many of the systematic errors which are possible in other 
methods. The other determinations were made by the observation of Fizeau 
fringes by a method which was basically that of Késters in which fringe 
matching is used and in which the refractive index of the air is automatically 
eliminated. Both of these methods had a precision about the same as that 
of the Fabry—Perot, viz., about +.0003 A. standard deviation in single 
determinations. 

These measurements were completely independent of those made by the 
Fabry—Perot and were made by a different observer, using a different mono- 


chromator, lamps, ete. 


SOURCES 


The Hg 198 lamps were of the Meggers ‘reservoir’ type containing Hg 198 
of 99° isotopic purity and argon carrier gas at a pressure of 0.25 mm. The 








838 CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 


discharge region of the lamp was cooled to about 4° C and viewed from the 
side through the water jacket. It was excited gently by an oscillator of the 
type described by Smith (1957) operating at roughly 250 Mc/sec. 

Under these conditions the correction for pressure shift of the Hg 198 lines 
was .00002 A to the red for the lines reported; that for self-absorption is 
known to be negligible in the fourth decimal in angstroms (Terrien 1958a) ; 
the wavelength shift which may be associated with the excitation current 
was considered negligible in the present measurements on the basis of previous 
observations of the effect of excitation of much greater intensity than that 
used in the wavelength determinations. The effect of isotopic impurity can 
also be shown to be negligible in the fourth decimal in angstroms (Terrien 
1959). 

The Kr 86 lamps were similar to the hot cathode type developed at the 
German P.T.B. and contained only Kr 86 of 99.9% purity. They were cooled 
by immersion in a bath of Ne» at its triple point and excited by a current 
density of 2 ma/sq. mm in the thin-walled capillary in which the discharge was 
observed. The discharge was viewed from the side in the case of the Fabry-— 
Perot measurements and from the end in the direction cathode to anode 
for the Michelson measurements. 

The magnitudes of the wavelength shifts for this source were found by 
Engelhard (1958) to be as follows: the values for the Doppler shift when the 
discharge is viewed in the direction cathode to anode is independent of the 
current density and has the values: 


A) 6057 = 87X10-*A to the red 


and 
Ad 5651 = 74X10-°A to the red. 


The magnitude of the Stark shift, which is to the red, follows the law 
AX = const. (p.j)?* where p is the pressure of the Kr and j the current 
density. The pressure is calculated from the Meihuizen formula for the vapor 
tension of Kr at the temperature of the coldest part of the lamp. The value 
reported for the shift under the conditions of Engelhard’s determination 
correspond to values for the constant as follows: 


AX6057 = 240 (p.j)?* (7 in ma/sq. mm; p in mm Torr), 
AdXs5651 = 160 (p.7)**. 

These shifts were measured by the authors for the 6057 line by a pressure- 
scanning technique described previously (Baird and Smith 1957) and also by 
means of observations with the Michelson-type interferometers. The results 
are as follows: 

Doppler shift Ad 6057 = 48+9X10-6 A, 
Stark shift Ad 6057 = 200+30 (p.7)?”. 


The values obtained by the two methods agreed to within the experimental 
error indicated above. The Stark shift for the lamp as operated for the wave- 











| 
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length determinations amounts to 32X10-*A according to our results as 
compared with 38x 10-® A as found by Engelhard. The differences between 
our values and those of Engelhard are not very significant for the measure- 
ments reported here; however, it is hoped to determine these corrections 
considerably more accurately in the near future by the methods referred to 
earlier. This measurement is made difficult by the fact that the intensity 
and sharpness of both the 6057 and 5651 lines decrease at higher p.j and the 
intensity decreases at lower p.j thereby limiting the range over which AX 
may be measured precisely. It is useful to note that if the lamp is observed 
in the direction anode to cathode the Doppler shift and Stark shift tend to 
compensate each other, the compensation being exact at about 4.0 ma/sq. mm. 


RESULTS AND ACCURACY 

The best vacuum wavelength values derived from all the observations are 
given in Table III for unperturbed conditions. The values include the observa- 
tions made by means of the Michelson interferometers in the cases where 
such observations were made. The results obtained by the two methods 
(Fabry-Perot and Michelson) are shown separately in Table IV; in combining 
them they were weighted according to their respective statistical accuracies 
(which were approximately the same). The differences between the results 
obtained by the two different methods are not considered significant. 

The stated error limits in Table III are the calculated 99% confidence 
limits, i.e. about three times the calculated root-mean-square error for the 
average of all the observations on a given wavelength. 

Insofar as systematic errors are concerned the following is relevant: the 
6057-A reference line with the corrections which were applied can be con- 
sidered to have been accurate to +.00002 A; the same accuracy applies to 
the Hg 198 green line. Errors as great as .0001 A are known to be very 
unlikely in the production of the other lines reported (Baird and Smith 1957; 
Engelhard 1959). The diversity of measurements which were made indicates 
strongly that the results are free from significant systematic errors of observa- 
tion. In addition to the experiments described it may be added that the first 
two Kr 86 lines and first three Hg 198 lines in Table III as well as the 6057 
line of Kr 86 have been used in this laboratory in a large number of precise 
length measurements whose results give no evidence of systematic error in 
the values reported. 

In view of these considerations the values of Table III are presented with 
considerable confidence as being accurate to within the stated error limits. 
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MODIFICATION OF PbS NOISE SPECTRA BY RADIATION! 


R. L. WILLIAMs 


ABSTRACT 


The noise spectra of PbS photoconductive cells can be modified by optically 
excited carriers. Detailed wavelength studies at CO. temperatures have shown 
that g-r noise can be obtained using light of the 2.4-micron region, and the 
1/f form with radiation of the 1.0-micron region. The effect is explained by 
assuming that minority carrier lifetimes in Kodak Ektron detectors are very 
short. The 1/f noise is only observed when optical carriers are produced very 
near the cell surface. 


INTRODUCTION 
Noise spectra of Kodak chemically deposited PbS photoconductive cells 
have been measured with excess incident radiation. The form of the noise 
spectra has been found to depend on the wavelength of the radiation. 


EXPERIMENTAL TECHNIQUE 

Wavelength studies were carried out at CO, temperatures using an Axler 
linear wedge filter. Wavelength bands 0.05 micron wide were selected between 
1.0 and 2.4 microns. A typical set of measurements was made as follows. A 
series of noise spectra was taken with various levels of illumination falling 
on the cell, but using a fixed wavelength interval of the filter. A second set 
would then be taken in a similar fashion but using a different wavelength 
interval. The resistance ranges covered by the two sets of measurements were 
comparable. 

A schematic diagram of the cell cryostat is given in Fig. 1. The measurements 
at CO. temperatures were performed with the cell chamber immersed in a 
mixture of dry ice and acetone. For most of the measurements air was trapped 
inside the system to assist in maintaining the cell at the temperature of the 
bath. Noise measurements were made between 10 c.p.s. and 16 ke.p.s. with 
a General Radio Wave analyzer. The signal was tapped off one stage before 
the meter circuit and passed through a Ballantine r.m.s. voltmeter. The signal 
from the mean-square output terminals of the Ballantine meter was recorded 
for at least 30 seconds for each measurement. The values obtained from 
the recorder, being proportional to the mean square of the noise voltage of the 
photoconductive cell, were corrected for the zero current noise by a simple 
subtraction. 

The light beam used to produce resistance changes was modulated by a 
chopping disk with two narrow blades made of wire mesh. The response time 
of the cells, to this modulation, was recorded to determine the time constant 
corresponding to the various radiation conditions. 

1Manuscript received January 6, 1959. 

Contribution from the Canadian Armament Research and Development Establishment, 
Valcartier, Que. 
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Fic. 1. Cell cryostat used in the investigation of the noise-wavelength dependence. 
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Fic. 2. Noise frequency change-over from the f~! to the f~-? form (CO2 temperatures). 
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EXPERIMENTAL RESULTS 

Figure 2 is representative of the results that could be obtained with a cell 
operating at CO2 temperatures. The cell resistance R, in megohms, the cell 
voltage V, in volts, and the cell time constant, in milliseconds, are recorded 
on the graphs. Curves B and C were obtained by changing the level of illumina- 
tion to obtain arbitrary resistance values. The wavelength band in this case 
was between 1.75 and 1.80 microns. The curves shown are three of eight taken 
during this particular experiment. In general, the resistance was reduced by 
arbitrary amounts down to the minimum obtainable. When the same pro- 
cedure was followed, using the wavelength region 1.00 to 1.05, only the 1/f 
form was obtained. By contrast, radiation in the 2.35 to 2.40 region produced 
domination of generation-recombination (g-r) noise when the resistance 
changed by only a factor of 2. 

The results obtained with this cell can be expressed as follows. With suffi- 
cient radiation incident on the cell to maintain the resistance at approximately 
2 megohms, a 1/f noise spectrum can be obtained by using wavelengths in the 
1.0- to 1.5-micron region. Using wavelengths longer than 1.7 microns, the 
g-r form was always obtained. The wavelength at which the change from 
dominance by 1/f to g-r noise occurred was a function of the resistance of 
the cell. When the cell resistance was about 9 megohms, 1/f noise was domi- 
nant when the cell was irradiated with radiation in the 1.8-micron region, but 
a good g-r noise spectrum was obtained when the 2.4-micron region was 
used. The resistance of the cell with no external light source was 31 megohms. 

The particular wavelength dependence mentioned above is typical of the 
cells studied, but other cells had in detail a different wavelength noise charac- 
teristic. At CO: temperatures some cells did not show such a wavelength de- 
pendence; however, such cells were completely dominated by 1/f noise at all 
temperatures and under all conditions investigated. By contrast, one cell 
required wavelengths of near 1 micron before 1/f noise appeared. A cell which 
had a room temperature spectrum which was partially of the g—r form usually 
showed the wavelength noise spectrum dependence. 

Figure 3 illustrates the noise spectral change that could be observed at 
room temperature. With no excess radiation this cell exhibited no 1/f noise 
at 10 c.p.s., curve B. 

At room temperature the intensity of radiation that could be obtained 
with the light source was too small to allow filtering of the radiation. A fila- 
ment lamp operating at full voltage was used with apertures controlling the 
intensity. The operating temperature of the filament being near 3000° K 
ensures a predominance of short wavelength radiation (less than 1.5 microns). 
The spectrum change observed is shown by curve C. This type of spectrum 
change was usually complete when the cell resistance had changed by a 


factor of 3. 
Fourteen cells were tested in this manner at room temperature. All cells 


under the radiation conditions discussed above were completely dominated 
by 1/f noise; three of these cells were dominated by 1/f noise without 


radiation. 
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Fic. 3. Noise spectral changes affected at room temperature by increasing amounts of 
radiation from a tungsten lamp. 


A general feature of the effect was that cells which had a very small 1/f 
component at room temperature required, at CO: temperatures, the shortest 
wavelengths to produce the 1/f noise, and vice versa. 

Both a globular and a filament lamp were used for the CO: experiments to 
check that the radiation source was not producing the 1/f spectrum. Identical 
results were obtained with the two sources. The fact that g—-r noise was ob- 
served using long wavelength light, from the same source whose short wave- 
lengths produced 1/f noise, indicates that the spectrum could not arise from 
the radiation source. The cell contacts were shown not to be responsible for 
any of the effects, as identical results could be obtained with the contact area 
masked with black electrical tape. 

The noise voltage squared dependence was checked and found to be /'*. 
Both the g-r and the 1/f noise components followed the same law. Discussing 
Fig. 2, it was pointed out that maintaining the cell resistance constant, which 
ensures the same voltage across, and current through, the cell, the spectrum 
could be changed. The effects cannot be associated with the particular current 
values used. 

The dominance of the 1/f noise occurs because of an increase in magnitude 
of the 1/f noise component. This is demonstrated in Fig. 3. That the g-r noise 
component is below the 1/f component can be seen by using the formula for 
g—r noise as given by Petritz (1956). The validity of the noise formula has been 
established by Lummis and Petritz (1957). Using Wood’s (1957) result that 
the mobility does not change under illumination, the ratio of the noise voltage 
squared for the condition of curve C, subscript 2, to that for the conditions 
of curve B, subscript 1, is 


Us a (Ye Bote) Ree 


V., Ro. +R Re, 71 ; 
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In this expression V? is the mean-square noise voltage, V, the bias voltage 
across the cell, Ry the load resistor, and 7 the time constant. Using the values 
of the parameters given on Fig. 3 and the fact that Ry was equal to 1 megohm, 
the reduction of the g—r noise, as predicted by the above equation, is a factor 
of 16. Noise of the g—r form reduced by this factor might just be detectable, 
but is not apparent in the 1/f noise of curve C of Fig. 3. 

A similar analysis, attempted on the noise spectra at CO, temperatures, 
was unsuccessful. Using the expression for g-r noise employed to obtain 
equation 1, the equation for the mean-squared noise voltage can be written 
in the form 


— 
* 
t 


— (BA) Beet 
™ Re+Ry/ 1+w's? * 


Here, e is the electronic charge, » the mobility, a a geometrical factor, and 
Af the bandwidth. This expression was fitted to g-r noise curves obtained 
from the same cell as used to obtain the data of Fig. 2, but using radiation 
in the 2.4-micron wavelength region. Values of wa were found to be a strong 
function of R, at these temperatures, and the analysis was not continued. 


INTERPRETATION 

Noise of the 1/f form is believed to arise from trapping of carriers in slow 
surface states (Kingston 1956), which exist in or on the outside of a surface 
oxide layer (Lasser et a/. 1956). Wang and Wallis (1957) have established 
that radiation can alter the surface potential of a semiconductor. Trapping 
effects in slow states, which are believed to give rise to 1/f noise (Kingston 
1956), would certainly be affected by any change in potential at the surface. 
That the magnitude of the 1/f noise changes with radiation intensity, there- 
fore, is not too surprising. The wavelength dependence of the noise must arise 
because of some additional feature, dependent upon the wavelength of the 
light employed. 

Petritz et al. (1956) have suggested that the high sensitivity of PbS films 
may arise from the rapid trapping of the minority carriers into states which 
do not readily trap majority carriers, in this case, holes. Short minority carrier 
lifetimes prevent any appreciable diffusion of electrons. Under conditions of 
space charge neutrality the distribution of optically excited free holes is similar 
to that of the optically generated (trapped or free) minority carriers. The 
minority carrier lifetime in unsensitive films is of the order of 10-7 second 
(Petritz et al.) and in the sensitization process the minority carrier lifetime 
is believed to decrease to some small value. Assuming a value of 10~!° second 
and using Lummis and Petritz’s (1957) result of 2 cm*/v-sec for the mobility 
of carriers, the room temperature diffusion length is .006 micron. As this is 
much shorter than the reciprocal of the absorption coefficient, the distribution 
of trapped and free carriers is approximately of the form exp —A,, where A 
is the absorption coefficient for radiation of wavelengths considered, and x 
is the distance from the surface. 
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The noise spectra changes were observed when the resistance of a cell had 
changed by a factor of 2 or 3. This indicates that the effects are important 
when the optically excited carriers are dominating the conduction process. 
Under the condition of short minority carrier lifetime mentioned above, the 
distribution of carrier approximates to Pyexp— A,, where Py is the density of 
carriers per square centimeter at the surface. Neglecting mobility changes, a 
constant resistance requires that the integral of Poexp—A, across the film 
be a constant. This implies that Po will be some function of A, i.e., if the film 
is thick compared to A~, Po is proportional to A. The surface potential must 
then be different for different wavelengths of light. Po is expected to be some 
function of A approximating to A~, as the films are of the order of 0.6 micron 
thick (Koch 1956), and the absorption coefficient is about 2.7104 cm~ at 
1.8 microns (Avery 1954). 

It is to be noted that the reciprocal of absorption coefficient of PbS for the 
radiation in the wavelength region used (A ~ 2.7X104 cm~') is at least a 
factor of 10 larger than the room temperature Debye length. This indicates 
that the effects are not due to the position of excitation in the region of changing 
surface potential. 

For an increase of the carrier density at the surface the probability of a 
transition to a slow surface state should increase. This is expected as the 
capture probability is normally proportional to the number of particles avail- 
able for capture. A high probability of capture means that a larger percentage 
of the total number of carriers is involved in surface trapping and a larger 1/f 
noise component is expected. Lead sulphide has a larger absorption coefficient, 
the shorter the wavelength of light. Then on the basis of the above arguments, 
the shorter the wavelength the larger the 1/f component. Different surface 
conditions of different cells would give rise to a change, in detail, of the wave- 
length noise spectrum dependence observed. The different wavelength be- 
havior, cell to cell, is expected on such grounds. For a particular cell, different 
resistances imply different carrier densities and different potential configura- 
tions at the film surface. Different wavelength noise spectrum relationships 
are then expected for the same cell but for different resistances. 

The photons of short wavelength, as well as penetrating shorter distances 
into the film, are more energetic. If the difference in energy is considered as 
significant, the effects could result from excitation of some surface state which 
affects the surface trapping. Bray and Cunningham (1958) have found that 
rather narrow wavelength regions excite surface states on germanium. The 
effects, observed at liquid air temperatures, have been identified with slow 
surface states. In the present experiments the wavelengths required to produce 
a change of noise spectra was dependent on the level of radiation. This makes 
an explanation based on excitation of surface states unattractive, but it is 
appreciated that the first explanation suggested involves in some cases specu- 
lative values of physical parameters. 
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PARAMAGNETIC RESONANCE LINE SHAPES! 
PREM SWARUP? 


ABSTRACT 


A study of the line shape of the transitions of various concentrations of Cr*++* 
ion in potassium cobalticyanide and potassium aluminum alum single crystals 
has been made at 9300 Mc/s. The line shape is found to change gradually from 
Lorentzian to Gaussian with chromium concentration. Lorentzian shape of 
gadolinium ion in lanthanum ethyl sulphate and exchange narrowing in case of 
free radical and copper sulphate are reported. The line width and relative 
intensity of the transitions in a highly dilute cyanide salt at 4.2° K are plotted 
and it is found that the intensity followed the maximum of the matrix element 
squared for various transitions. The value of g is found to be 1.992. The actual 
concentration of Cr*** in dilute salts is less than the concentration of the solution 
from which they are grown by more than a factor of 10 as the dilution is 
increased. 


I. INTRODUCTION 


The paramagnetic behavior of the triply ionized chromium in potassium 
chromicyanide has been studied as well as other ions of the transition groups 
in an attempt to determine the shape of the resonance lines as a function of the 
concentration of the paramagnetic ions and of the microwave power. It is well 
known that at microwave frequencies both storage and dissipation of the 
magnetic energy are involved when a wave is propagated through a para- 
magnetic sample immersed in a magnetic field. The magnetic properties of the 
medium are described in terms of a complex magnetic susceptibility: 

x= x'—jx" 
where the imaginary part x” is responsible for the absorption of the energy 
in the medium. The expressions for the two quantities x’ and x”’ are the same 
as those given by Kronig—Kramers in the case of electrical susceptibility. 
The r-f susceptibility, x’’, at a frequency » is given by: 


ad 


x’ (vr) — Txovog(v) 
where xo is the static susceptibility, vo is the resonance frequency, and g(v) 
is the shape function which is normalized according to the relation 
fog r)dr = 1. The Bloch susceptibility expressions which hold well in the 
paramagnetic resonance are given by 


ay ) as eee 27. 
Ee Sade Tie) 
( ) o 427 2(vo—v) 


, 7 
x’ (v) = rxwol 2 ———3-3 3 . 
( 1+42°T3(vo—v) 

‘Manuscript received April 6, 1959. 
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The Bloch equations in absence of any saturation give the line a Lorentzian 
shape as is well known in collision-broadened spectra. The broadening of the 
paramagnetic resonance line can be pictured as caused by neighboring spins 
limiting the lifetime of a given spin state by inducing transitions from it, 
which is mathematically parallel to the collision broadening of the spectral 
lines resulting in Lorentz shape function: 
2T>2 

@) £) = Trae THe)" 
where 7» is the ‘transverse’ or spin-spin relaxation time. 

On the basis of the dipolar broadening of the magnetic resonance lines, 
Van Vleck (1948) derives expression for the moments of various orders of 
the shape functions g(v). The (2”)th moment is defined as: 


(3) ((Av)™) 4) = Jee) (Av)™ dy. 


An experimental check of his theoretical calculations to show that the nuclear 
resonance lines have essentially dipolar broadening, and that the lines have 
Gaussian shape, was successfully done by Pake and Purcell (1948) in the case 
of CaF,. The ratio (Av*)V/*/(Av*)\?, calculated theoretically as 1.25, was 
close to the experimental value of 1.24+0.02. 

Kittel and Abrahams (1953) have extended the calculations of Van Vleck 
to the magnetically dilute samples. If f is the probability that a lattice site 
of a nonparamagnetic crystal is occupied by a magnetic system, then the 
fourth moment of the shape function (assuming a simple cubic lattice with 
Ho along (100) axis) is: 

(4) (Av*)ay = 3[(Av”),y]°{0.742 +f [0.098 — 0.021 (.S°+.S)~"]} 
= 3[(Av*)ay]” for the Gaussian shape. 


For S = 1/2, the shape is Gaussian for f > 0.1 
and for S = 3/2, the shape is Gaussian for f > 0.4. 


Broer (1943) first tried to evaluate the higher moments to determine the shape 
function. He used Heisenberg’s theory of ferromagnetism in taking the shape 
function proportional to the Gauss function. The susceptibility would then 
have the form: 

x’ (v) = 2rxwol 2 exp[—4473(vo—v)’]. 


Pure dipolar interaction between the spins, in the absence of exchange 
effects, crystalline fields, or quadrupole moment, gives resonance line a 
Gaussian shape. This is true both in the case of nuclear magnetic and electron 
paramagnetic resonance. If it is assumed that the only interaction between 
the different spins are dipolar and exchange coupling, the Hamiltonian 
functions can be expressed as: 


§ = Zeeman + exchange + dipolar energies 


= Hgp>, Sut 2B , Set 2D B,US;. $)-On Gn. SdGa. S01 
J g 
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The constants A, = —2/] , and B, = £B/r x, J is the exchange integral, 
and 7 the interatomic spacing. S,; is the z component of the spin angular 
momentum vector of atom j. The presence of the exchange term (A) in the 
Hamiltonian increases the fourth moment relative to the square of second 
moment in the relation (4) resulting in the narrowing of the Gaussian line 
near the center and tapering gradually in the wings. The effect can be visualized 
by the fact that the exchange term A commutes with S and, therefore, does 
not contribute to the second moment. QDaipoiar and S do not commute and 
Dexchange does not commute with Daiporar, hence Dexchange Contributes to the 
fourth moment. (Reference is made to the expressions of various moments 
given by Van Vleck (1957).) It has been observed in cases of undiluted copper 
salts and manganese salts (MacLean and Kor 1956) that the calculated line 
widths are much greater than those observed. This exchange contribution 
to the fourth moment is important only in the undiluted paramagnetic salts 
since the exchange integral decreases rapidly with the distance between the 
paramagnetic ions. At very high magnetic dilution, on the basis of dipolar 
interaction, a Lorentzian line shape is expected. The half width is expressed as 
A = 5.3 f(g8u?/ha*) where the symbols have their usual meanings. The line 
width is linearly proportional to the concentration and the intensity is inde- 
pendent of the concentration. 


II. THE LORENTZ AND GAUSSIAN SHAPES 
The Lorentz and Gaussian expressions are: 
siecle 
1+42°T3(v—vo)*’ 
(6) x" (v) = Poexp{ —44T3(v—v0)"}. 


The spin-spin relaxation time 7» is related to the width of the line between 
the half-maximum points by the relations: Avy = 1/77, for Lorentzian and 


1 


(5) x") = 


Il 


Av = Vin2/x-1 72 for Gaussian lines. Substituting for 7, in the relations 5 
and 6, we get: 


(7) eo po/[r+a(") k 
(8) x (y= Ps exp| —4 In afte) | : 


The first derivatives of these are: 


dx''(v) _ | ex") 32V3 ea! osey 
(9) dv dy a) a Av ws Av 

dx"'(v) _ | exe) ated =) es fs a{z=z) "| 
(10) oe ayge ABR, 3.869 Ay ) oxP 4 In 2 res 


having maximum values when [(vy—v9)/Av] = 1/(2/3) in case for Lorentzian 


shape and when [(y—v)/Av] = 1/(2V2 In 2) in the Gaussian case. The 
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normalized graphs of the two derivatives, represented by 9 and 10, intersect at 
[(v—vo)/Av] = 1 when their magnitudes drop to quarter of their maximum 
values. Graph | shows the shapes of the absorption curves and their derivatives 
as represented by equations 7, 8 and 9, 10. 











0 oan al ai a s 


° 0.5 ' l 
V-% . 
av 


Fic. 1. Curves 1 and 2: Lorentzian and Gaussian line shapes of identical line widths. 
Curve 3: One side of the Lorentzian line derivative. Curve 4: One side of the Gaussian line 
derivative. The plot is between the normalized amplitude and the normalized frequency 
deviation from the resonance frequency. 


III. EXPERIMENTAL TECHNIQUES 

A conventional microwave paramagnetic-resonance spectrometer utilizing a 
magic T impedance bridge was used. A silver-plated cavity operating in 
TEj02 mode was coupled to the end of a two-foot-long, stainless-steel, X-band 
waveguide. The whole assembly passed through a polystyrene head such that 
the cavity was centered between the poles of an electromagnet inside a double 
Dewar assembly. The cavity could be cooled to temperatures down to that of 
liquid helium. A 1000 c/s modulation on the magnetic field was kept below a 
peak to peak value of 0.5 oersted. The static magnetic field was swept at the 
rate of 2 oersteds/minute and was controlled by the driving motor of a 
Rectiriter recorder. The operating frequency was 9300 Mc/s. The sample 
volume was kept roughly to 0.15 cc. In the angular rotation experiments 
the crystal was hung at the end of a glass rod passing through the waveguide 
and the coupling iris. A graduated disk mounted on top of the rod indicated 
the orientation of the crystal. The microwave power entering the cavity was 
monitored. The Rectiriter recorded the differential absorption curve. A 
Sperry 2K39 Klystron giving out about 500 milliwatts of power was used in the 


saturation experiments. 
The paramagnetic ion studied in detail is Cr+**+ in KsCr(CN).¢ diluted by 
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the isomorphous K;Co(CN).. It is a monoclinic crystal having its crystalline 
axis 6 perpendicular to the axes a and c, while K;Co(CN)¢ is an orthorhombic 
crystal with a # 6 ¥c and the three axes are perpendicular to each other. 
Although they differ slightly in their crystalline forms, they crystallize 
isomorphously. There are two magnetic complexes per unit cell which are 
oriented differently but their energy levels are identical when the magnetic 
field lies in either the (ac) or the (bc) planes of the crystal. The X-band 
transitions involving JJ = +1 and M = +2 (observed owing to a mixing 
of quantum states) at various angles of rotation are computed from the energy 
level diagrams. The four energy levels are derived from the spin Hamiltonian: 


© = g8 (H,S,+H,S,+H,S,)+D(S?—15/4)+E(S?—S}) 


with S = 3/2, D = 2.49kMc/s, E = 0.33 kMc/s, and g = 1.992. The energy 
levels can alternatively be determined by the methods of group theory and 
time inversion considerations as proposed by Koster and Statz (1959). 

Potassium chromicyanide is diluted with the isomorphous potassium 
cobalticyanide in various proportions. The two salts were dissolved in distilled 
water in various proportions by weight and the crystals were grown by slow 
evaporation. The nominal concentration of Cr in the crystals as given by the 
concentration of the solution was found to be greatly different from the 
actual concentration. The concentration of Cr in the diluted samples was 
determined at Atomic Energy of Canada Ltd., Chalk River, by irradiation of 
the sample and counting the radioactivity. Table I lists the nominal and 
actual concentrations of the samples used in the experiment. It also lists the 
relative concentrations of the samples as determined by comparison of the 
paramagnetic intensity of absorption (total area of the derivative) of the 
(—1/2 > +1/2) transitions in K3;CrCO(CN). with a known sample of di- 
phenyl picryl hydrazyl. Care was taken in arranging the same orientation of 
the crystals. 


TABLE I 
Actual concentrations 
Nominal concentration, Chalk River results, Intensity 
Salt % Cr % comparison 
K;Cr(CN)6 20 2.02 10.00 
10 0.65 5.58 
7 3.40 
5 0.43 82 
1 0.02 045 
0.5 Ol 
Cr*+*+* in KAI(SO,4)2.12H2O 4 16 


IV. RESULTS AND DISCUSSION 
The derivatives of the absorption obtained experimentally are plotted on the 
graphs with reference to the center of the derivative. The ratio {[dx’’(v)]/dv} / 
{[dx’’(v)|/dv}max was plotted along the axes of ordinates. The distance in 
oersteds on the recorder chart divided by the distance at which the amplitude 
drops to one quarter of its maximum value (at (v—vo)/Av = 1) representing 
(H—H,)/AH was plotted along the abscissae. 





| 
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Potassium Chromicyanide 

Samples having various chromium concentrations have been examined. 
Graph 2 shows the experimental curves for the nominal concentrations of 
0.5%, 1%, 5%, 7%, 10%, and 20% of chromium. It was observed that 
nominal 0.5% and 1% salts gave truly Lorentzian lines which broadened 
gradually to Gaussian at high concentrations. Table II tabulates the results. 


TABLE II 


Line shapes of Cr+** transitions at 4.2° K 








0.4% 0.2% 


Actual dilution: 2% 0.6% 
0 0 


Transition Nominal dilution: 20% 10% 7% 5% 1% 0.5% 
—1/2— +1/2 G GL GL GL E L 
+1/2 > +3/2 GL GL GL E E 
—3/2— -1/2 GL GL GL L L 
—1/2 — +3/2 FP E, 
—3/2— +1/2 L L 


Nore: G: Gaussian, L: Lorentzian, and GL: between Gaussian and Lorentzian. 


Potassium Chromium Alums 

The shapes of the single resonance line corresponding to the (—1/2 — +1/2) 
transition in diluted potash chrome alums were studied. Undiluted chrome 
alums have been reported earlier to have Gaussian line shapes (Davis and 
Strandberg 1957). The crystals studied here had nominal Cr*+** concentrations 
of 4%, 10%, 20%, 50%, and 80%. Analysis of the 4% crystal gave the actual 
concentration of 0.16%, a factor of 25. Other crystals were not analyzed. 


Lor fs, 1.0r 





Pmax 


0.5; 


Pan 


0.5 








F Fic. 2. Shapes of the differential absorption curves for different concentrations of Cr+++ 
in potassium cobalticyanide. The curves are represented by: nominal concentration 0.5% 
and 1% ——; 5% and 7% —-—-; 10% — —; and 20% ——-. 

Fic. 3. Shapes of the differential absorption curves for different concentrations of Cr+++ 
in potassium aluminum sulphate. The curves are represented by: nominal concentration 10% 
by solid line, 20% by dot-dash line, and 80% by dotted line. The curve for 50% is closer to 
80% and is not shown. 
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Graph 3 shows the shapes of the derivatives. Nominal 80% salt had a true 
Gaussian shape, 50% (not shown) was found to be close to 80% curve, 20% 
lies between the two shapes, and 10% is closer to the Lorentzian shape. 


Gadolinium Ethyl Sulphate 

A dilute sample of Gd+++ in lanthanum ethyl sulphate (about 1%) had 
shown Lorentzian line shape for all of its seven AM = +1 transitions and six 
AM = +2 transitions at liquid helium temperature. Graph 4 shows the 
Lorentzian line shape and the experimental points. The line width of (—1/2 
— +1/2) transition was 20.8 oersteds. 
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Fic. 4. Shape of the differential absorption curve for 1% gadolinium in lanthanum ethyl 
sulphate. The experimental points are shown by small circles on the true Lorentzian curve. 


Exchange Narrowing 

Two cases of exchange narrowing were studied. A powdered sample of 
free radical (diphenyl picryl hydrazyl) showed a line shape between Gaussian 
and Lorentzian at room temperature. The shape was not much changed at 
4,.2° K even though the line width had increased to nearly two and one-half 
times its room temperature value of 2 oersteds. A compact crystal of copper 
sulphate showed a narrow resonance line resembling Lorentzian. 


Haif Widths of Lines 

The half widths AH (between half-maximum points) of the various transi- 
tions of the potassium chromicyanide crystal were measured at 4.2° K. The 
line widths remained practically constant for a particular transition on 
rotation of the crystal in the (bc) plane. The widths were measured with 
presumably no power saturation and no correction was made for the longi- 
tudinal relaxation time. The earlier reported r.m.s. widths of the undiluted 
crystal by Baker, Bleaney, and Bowers (1956) were 750 oersteds along a axis, 











SWARUP: RESONANCE LINE SHAPES 855 


240 oersteds along b axis, and 700 oersteds along c axis. The magnetic field 
spectra of the resonance lines at different dilutions were plotted and the value 
of g tensor was determined by comparing the position of the resonance line 
with that of the resonance line of diphenyl picryl hydrazyl (g = 2.0036). The 
value was found to be 1.992+.002 in all the cases as compared with the 
value of 1.998+.008 in the undiluted crystal. The value of g was not found to 
vary throughout the temperature range (4.2° to 80° K). The relative peak 
intensities and the line widths of the five transitions observed in the nominal 
0.5% sample at liquid helium temperature and at a crystal orientation of 90° 
to the direction of the magnetic field in the (dc) plane are shown in Fig. 5. 
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Fic. 5. Line widths and peak intensities, relative to (—1/2 — +1/2) transition, of the 
five transitions in 0.5% Cr*+** in potassium cobalticyanide crystal are shown in (a) and (6) 
respectively. The transitions are represented in order by (—1/2 — +1/2), (—3/2 — —1/2), 
(+1/2 — +3/2), (—1/2 > +3/2), and (—3/2 — +1/2). The dotted blocks in (6) represent 
the square of the matrix element relative to maximum for (—1/2 — +1/2). 


Since the intensity measurements were made in the case of all transitions 
under similar conditions of the microwave power and temperature, the relative 
peak intensities could be compared to the maximum of the square of the 
perturbation matrix element of the magnetic moment of Cr+++ corresponding 
to the various transitions (Chang and Siegman 1958). The calculated matrix 
element square is shown by dotted blocks. The deviation of the experimentally 
observed intensities with the calculated one is marked in the (3/2 > —1/2) 
and (—3/2 — +1/2) transitions. The line width showed no variation on 
rotation in the (bc) plane. Great care had to be taken in setting the orientation 
of the crystal because a slight deviation from the (bc) plane separated the two 
spectra due to the two magnetic complexes per unit cell, and the line width 


might appear to have broadened. 
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Line Shape at the Onset of Saturation 

The line shape of the nominal 0.5% Cr crystal (—1/2 — +3/2) transition 
had been studied at different power levels at the onset of saturation at 77° K. 
The shape remained Lorentzian although the width increased from its low 
power value. Holcomb (1958) had shown that in case of nuclear resonance, 
‘saturational’ narrowing was observed and the shape changed from Gaussian 
to Lorentzian. Further work in this direction will be reported later. 


V. CONCLUSIONS 


The concentration of the paramagnetic impurity in the isomorphous 
diamagnetic crystal as compared with the concentration in the mother solution 
did not appear to follow a definite rule. Giordmaine, Alsop, Nash, and Townes 
(1958) have reported that in gadolinium magnesium nitrate actual con- 
centrations were found to be less than the nominal values by somewhat more 
than a factor of 10. They found no such large differences in other materials 
examined such as copper ammonium sulphate. The present measurements in 
case of potassium chromicyanide indicated that at higher nominal concen- 
trations, the factor (ratio of nominal to actual concentration) was 10 which 
increased to 50 at very low concentration. A sample containing a nominal 
concentration of 0.5% Cr had given a line width for the (—1/2 — +1/2) 
transition of 9 oersteds as compared with Townes’ (1958) value of 12 oersteds. 
The same crystal had been successfully used as a Maser crystal in two setups 
in this laboratory operating at 10 cm and having same operating characteristics 
as reported elsewhere. 

The shapes of the resonance lines of chromium salts were found to depend on 
the distribution of the paramagnetic ions in the isomorphous diamagnetic 
crystal lattice. The results could be explained on the basis of purely dipolar 
type of interaction between the spins. 
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A SIMPLE GAMMA-RAY INSENSITIVE FAST-NEUTRON 
COUNTER' 


G. M. Grirritus, P. P. Sincu, Y. I. Ssu, AND J. B. WARREN 


ABSTRACT 


A gamma-ray insensitive fast-neutron counter is described which consists 
of a number of thin rectangular sheets of lucite coated with zinc sulphide and 
sandwiched together to form a rectangular block which is mounted on a photo- 
multiplier. Pulse height spectra and absolute efficiency curves are presented 
for neutrons from 300 kev to 15 Mev. 


INTRODUCTION 

Chadwick in 1932 first identified fast neutrons by observation of the knock- 
on protons ejected from hydrogenous materials. Since that time many neutron 
counters have used the same neutron—proton interaction in either hydrogen- 
filled proportional counters (Coon and Nobles 1947), ionization chambers 
(Stafford 1948), or solid hydrogenous materials. Bell (1948) and Segel et 
al. (1954) employed organic phosphors successfully ; however, these suffer from 
the fact that they are gamma-ray sensitive if the volume used is sufficiently 
large to give good neutron-detection efficiency. A thin layer of zinc sulphide 
powder placed on the face of a photomultiplier has been used to detect proton 
recoils from a layer of plastic resulting in a neutron counter with a very low 
gamma-ray sensitivity; however, the neutron efficiency is also small. Hornyak 
(1952) made a great improvement in detection efficiency while maintaining 
low gamma-ray sensitivity by using a molded button of lucite with zinc sulphide 
powder dispersed in it. As zinc sulphide is opaque to its own radiation the 
useful size of this device is limited. Since then several attempts have been 
made to improve the neutron-detection efficiency without increasing the 
gamma-ray sensitivity (Emmerich 1954; Brown and Hooper 1958). Recently 
a method, based on the fact that the ratio of the intensity of the slow light 
decay component to that of the fast component is different for protons and 
electrons in organic scintillators, has been used to discriminate between gamma 
rays and neutrons (Litherland et a/. 1959; Owen 1958). 

Below we describe the properties of a counter used in our laboratory for 
several years, which has an efficiency comparable to others we have seen 
reported, along with an excellent gamma-ray rejection ratio and simplicity of 
construction. 

COUNTER CONSTRUCTION 

The counter consists of a number of thin rectangular lucite sheets coated 
with zinc sulphide over their large faces and then sandwiched together to form 
a rectangular block which is mounted on an RCA 6342 photomultiplier. Pro- 


1Manuscript received February 24, 1959. 
Contribution from the Department of Physics, University of British Columbia, Vancouver, 
B.C. 
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tons knocked out of the lucite produce light in the thin zinc sulphide layers 
and the lucite plates conduct the light to the photomultiplier. To make the 
block, the lucite sheets were dipped in ethyl formate and while still wet a 
weighed quantity of zinc sulphide powder was sprinkled on evenly through 
a fine mesh nylon cloth. Then the plates were stuck together under pressure 
until dry. The block was then polished on a face perpendicular to the plane 
of the sheets and this face was mounted on the photomultiplier with Dow— 
Corning silicone oil. The other faces of the block were smoked with burning 
magnesium ribbon. The construction of a typical phosphor block is shown in 
Fig. 1. 


SE SEE, 
 oneeeenetietestapeenmeneetememmmemmmeeetmeene eee 
oenenneneniaieesihdateiteteeaintemmmmemenemmemeee ee 


Fic. 1. The lucite — zinc sulphide sandwich block. 


Both Patterson D and RCA 33-Z-20A zinc sulphide phosphors were tested. 
With a radium—beryllium neutron source the RCA phosphor gave about 30% 
more counts above a selected bias than the Patterson D phosphor and so was 
used for all further work. Tests to determine the optimum thickness of the 
zinc sulphide powder indicated that for thicknesses greater than 12 mg/cm? 
the gamma-ray sensitivity becomes appreciable and thicknesses less than 
10 mg/cm? gave a satisfactory neutron response. Several lucite sheet thick- 
nesses were tried and it was found that for radium—beryllium neutrons the 
efficiency decreased for thicknesses greater than 0.125 inch. The plate thick- 
ness determines a limit to the height of the block that can be usefully em- 
ployed, since for a given height the light from the top of the block will not 
reach the photomultiplier as efficiently with thin plates as it would with 
thicker ones due to absorption in the zinc sulphide. Thus with increasing 
height a decreasing proportion of the light will get into the photomultiplier 
and after a certain height, which depends on the plate thickness, any additional 
height contributes very little light. However, greater plate thickness, although 
it means that the detector may be made longer, does not necessarily mean 
greater efficiency since for thicker plates more of the protons lose all their 
energy in the lucite. Protons stopped in lucite alone produce very small 
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pulses which are lost in the noise. A check on the transmission of light down 
the block was made by preparing a plastic sandwich with 0.125 inch thick 
plates 5 inches long mounted on a photomultiplier. A radium—beryllium source 
was placed 25 cm away to one side and the counting rate was observed, as a 
function of the length, which was altered by cutting sections from the end of 
the sandwich block. Little change in counting rate was observed until the 
length was reduced to 3 inches and the counting rate became roughly propor- 
tional to length for lengths less than 1.5 inches. 

A rectangular block consisting of 1/8 inch square lucite rods packed to- 
gether with zinc sulphide powder on all internal faces was made and tested. 
Its efficiency per unit volume for neutron detection was little if any better 
than that of the flat sheet assemblies and as it was more difficult to construct 
no further ones of this type have been made. 

Tests were performed to find out whether the neutron-detection efficiency 
of the counter (using radium—beryllium neutrons) depended on the direction 
of neutron entrance (i.e. parallel or perpendicular to the lucite sheets). One 
might expect such an effect since the protons are scattered forward in the 
laboratory system, but within the errors of about 10% no effect was found 
either in the shape of the spectrum or in the number of counts obtained above 
a given bias. 

Two blocks with the following characteristics were tested in some detail: 


Block I Block II 








Lucite thickness 1/32 inch 1/16 inch 
Length lcm 2 cm 

Cross section 4.2X4.2 cm 4.44.4 cm 
No. of plates 50 28 

ZnS thickness 7 mg/cm? 7 mg/cm? 





The efficiency per unit volume for block I was 1.75 times that for block II. 
However, since its volume was less than half of that of block II it gave fewer 
counts for a given neutron flux. The results described below were obtained 
with block IT. 
COUNTER PERFORMANCE 

(a) Pulse Height Spectra 

The pulse height spectra and absolute efficiencies for the neutron counter 
were measured for neutrons of energies from 225 kev to 4.4 Mev using neutrons 
from the T(p, ~)He® and D(d, m)He’ reactions and for 17-Mev neutrons from 
the T(d, m)He* reaction. The differential pulse height spectra are shown in 
Figs. 2 and 3 along with spectra for photomultiplier noise and y-rays of 1 
Mev from Co® and 6.14 Mev from the F'*(p, a, y)O" reaction. The pulse 
height scale is the same for all curves but the vertical scale is arbitrary and 
the curves are not normalized to each other. 


(b) Absolute Efficiency 
The absolute efficiencies were determined from knowledge of the neutron 
fluxes from the targets used. The tritium content of the tritium—zirconium 











GRIFFITHS ET 


COUNTS 


COUNTS 


° 400 


200 


800 


AL.: GAMMA-RAY INSENSITIVE COUNTER 


940 KEV 


280 KEV 


300 400 500 600 
PULSE HEIGHT 


1200 1600 2000 2400 
PULSE HEIGHT 





700 





2800 


861 


Fic. 2. Pulse height spectra for low-energy neutrons and for gamma rays. The curves are 
not normalized with respect to each other. 
Fic. 3. Pulse height spectra for high-energy neutrons. The curves are not normalized with 


respect to each other. 
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target, kindly supplied by the Oak Ridge National Laboratory, was measured 
in terms of the absolute gamma-ray yield from the T(p, y)He* reaction at 
800 kev, using absolute differential cross-section data of Perry and Bame 
(1955) for this reaction and theoretical estimates of the scintillation counter 
efficiency. Low-energy neutrons from 280 kev to 940 kev were produced by the 
T(p, n)He® reaction and high-energy neutrons were produced by the T(d, 
n)He* reaction using the same target. Intermediate energy neutrons from 
1.9 Mev to 4.2 Mev were produced by D(d, m) He? reactions using a heavy ice 
target of known thickness. The D.O targets were made by condensing heavy 
water vapor from a constant volume dispenser onto a liquid air cooled gold 
plate. The dispenser—target system was calibrated by observing the shift in 
the 340-kev resonance of the F!9(p, a, y)O'® reaction after a given pressure of 
vapor was deposited on a thin F™ target. From the measured energy loss of 
protons in passing through the ice layer and from the stopping power data 
given by Whaling (1958) the number of target deuterons per square centi- 
meter was calculated. Neutron yields were obtained from cross-section data 
given by the review article of Fowler and Brolley (1956). 

The neutron-detection efficiency for the counter is dependent on the neutron 
energy as shown in Fig. 4. Curve (a) shows the angular distribution of 
D(d, n)He* neutrons for 1.0-Mev incident deuterons as obtained by Hunter 
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Fic. 4. The angular distribution of D(d, m)He® neutrons at 1.0 Mev. Curve (a): The dis- 
tribution of Hunter and Richards (1949) obtained with a long counter. Curve (5): The dis- 
tribution as seen by the lucite — zinc sulphide sandwich counter. Curve (c): The ratio of the 
curves (b) and (a), which are normalized to each other at zero degrees. 
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Fic. 5. Absolute efficiency vs. discrimination level for the low-energy neutrons. The ordinate 
for the gamma-ray curves should be multiplied by 1072. 


Fic. 6. 


Absolute efficiency vs. discrimination level for the high-energy neutrons. 
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and Richards (1949) using a long counter and making small corrections for the 
change in efficiency with neutron energy as given by Hanson and McKibben 
(1947). Curve (6) shows the curve obtained with the present counter normalized 
to curve (a) at zero degrees. The ratio between the curves is shown by the 
upper line and indicates the decrease in efficiency for decreasing neutron 
energies. The absolute intrinsic efficiencies defined as the number of counts 
observed divived by the number of neutrons incident on the front face of the 
counter for various neutron energies as a function of discriminator bias level 
(integral bias curves) are shown in Figs. 5 and 6. The efficiency falls off 
rapidly for neutrons with energies less than 300 kev and is very small for 
thermal neutrons. It should be noted that the efficiency figures given here are 
smaller by a factor A, equal to the area of the front face of the counter, than 
those quoted by some authors who give the number of counts per unit neutron 
flux at the position of the counter. Figure 5 also shows on an expanded vertical 
scale the absolute efficiencies for Co® and 6-Mev gamma rays. Note that for 
the gamma-ray curves the ordinates should be multiplied by 10-*. The curves 
indicate that a bias may be chosen without serious loss of neutron counts to 
give gamma-ray efficiencies of 10-* to 10-"° of the neutron efficiency at that 















bias. 

As a crude first approximation one would expect the neutron-detection 
efficiency to be proportional to the neutron—proton scattering cross section, 
o, and the mean range, Rk, of the scattered protons in the lucite, since for 
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Fic. 7. Absolute efficiency vs. neutron energy for various biases. The thick curve shows a 
rough theoretical estimate of efficiency with energy on an arbitrary vertical scale. 
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all energies studied here the proton ranges are smaller than the thickness 
of the lucite plates. The validity of this assumption is indicated in Fig. 7, 
which shows experimental absolute efficiencies as a function of neutron energy 
for several discriminator levels and also shows a curve of the product ¢-R asa 
function of neutron energy with arbitrary vertical normalization. The agree- 
ment in shape of the curves is quite good; the theoretical curve is rather 
flatter than the experimental ones suggesting that for lower-energy neutrons 
a larger proportion of the pulses are lost in the noise. 

This counter has been used to detect the small yield of secondary neutrons 
produced by proton bombardment of heavy ice targets below the D(p, 2) 
threshold in the presence of a much larger yield of 6-Mev y-rays, and as a 
general monitor for fast neutrons. 
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THE NEUTRON YIELD FROM HEAVY ICE TARGETS 
BOMBARDED WITH PROTONS BELOW THE 
D(p, n)2p THRESHOLD! 


P. P. Stncu, G. M. Grirritus, Y. I. Ssu, AND J. B. WARREN 


ABSTRACT 


In connection with some experiments on the gamma-ray yield from the 
D(p, vy) He?’ reaction using heavy ice targets a considerable yield of neutrons was 
found even for proton bombarding energies well below the D(p, 2)2p threshold 
of 3.3 Mev. The yield, excitation function, and angular distribution of this 
neutron yield have been investigated both experimentally and by means of 
theoretical calculations. These studies confirm the suggestion that the neutrons 
are produced by a secondary reaction in which deuterons, scattered in the 
target by incident protons, collide with further target deuterons to produce 
D(d, n)He? reactions. 


INTRODUCTION 


The production of neutrons when heavy ice targets are bombarded with 
protons below the D(p, 2)2p threshold has been reported by Jennings et al. 
(1950), who suggested that the neutrons arise from deuterons, scattered by 
incident protons, colliding with other target deuterium atoms and producing 
D(d, n)He’ reactions. In connection with some measurements on the gamma- 
ray yield from the D(p, y)He* reaction this effect was found to contribute a 
considerable background (Griffiths and Warren 1955). Consequently we have 
made a detailed experimental and theoretical study of the effect. 


EXPERIMENTAL 

The gamma-ray insensitive neutron counter used in these experiments has 
been described previously (Griffiths et a/. 1959). It consists of a block of lucite 
sheets with thin layers of zinc sulphide powder sandwiched between them, 
mounted on an RCA 6342 photomultiplier. A bias setting which made the 
6-Mev gamma-ray sensitivity less than 10~° of the neutron sensitivity was 
used. The absolute efficiency for neutron detection varied with neutron 
energy from 0.15% for 2-Mev neutrons to 0.3% for 4-Mev neutrons at this 
bias. The pulse height spectrum was analyzed with a 100-channel Computing 
Devices of Canada kicksorter. 

The target chamber used is shown in Fig. 1. A gold-plated liquid air cooled 
copper plate formed the target backing. This could be turned to face an inlet 
hole through which heavy water vapor from the dispenser was admitted to 
the target chamber. For thin targets, the constant volume dispenser was 
evacuated, taps T; and T, were closed, and tap T; was opened to admit 
the required pressure of water vapor as measured by the oil manometer. Then 
tap T; was closed and T; was opened to admit the vapor slowly via the 

1Manuscript received March 12, 1959. 
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Fic. 1. Target chamber and heavy water vapor dispenser. The graph on the right shows 
the dispenser calibration in terms of target thickness versus dispenser pressure. 


glass-wool diffuser to the target chamber. The target thickness as a function 
of dispenser pressure was calibrated by freezing the vapor on top of thin 
fluorine targets for various initial pressures in the dispenser and measuring 
the shift in the F!*(p, a, y)O'® 340-kev resonance due to the energy loss of 
the protons in passing through the ice layers. From a knowledge of the 
stopping power of the protons in ice (Whaling 1958) the number of deuterium 
atoms per square centimeter on the target could then be calculated as a 
function of the manometer pressure (Fig. 1). When a thick target was required 
so that the incident protons would be stopped completely in the ice layer, 
taps T, and T; were opened for 2 minutes and then a measurement of the 
yield was made. The taps were opened for a further 2 minutes and a second 
yield measurement was made. When no increase in yield occurred for additional 
ice on the target, the target was judged to be ‘‘thick’’. 


RESULTS 

1. Thick Target Yield and Excitation Function 

A thick heavy ice target was bombarded with protons of energies from 
350 kev to 1.5 Mev and the neutron yield at zero degrees to the proton 
beam per millicoulomb of beam per steradian was obtained. Since the detection 
efficiency is a function of the neutron energy and since at any angle the 
neutron energies incident on the counter varied from approximately 2 Mev 
to 4 Mev, the mean neutron energy at 0° was estimated for each bombarding 
energy by assuming that on the average the deuterons were scattered at 45° 
to the proton beam direction and then by calculating the neutron energy in 
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the direction of the counter from the mechanics of the D(d, )He?’ reaction. 
The experimental excitation function, shown in Fig. 2, has been corrected 
for the change in efficiency due to the change in mean neutron energy. The 
theoretical neutron yield as a function of proton energy calculated on the 
basis of the secondary process is also shown in Fig. 2. 
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Fic. 2. The theoretical and experimental thick target excitation functions. 
2. Thick Target Angular Distribution 


With the counter subtending a solid angle of 0.15 steradian at the target 
the thick target neutron vield was measured for a proton energy of 1.5 Mev 
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Fic. 3. The thick target angular distribution for 1.5-Mev incident protons. Curve 
(a), experimental results—uncorrected; curve (b), theoretical result; curve (c), experimental 
results—corrected for change in neutron-detection efficiency with mean neutron energy. 
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at seven different angles. The neutron yield observed at each is shown by the 
lower curve in Fig. 3. Since, for each counter position, neutrons from deuterons, 
scattered with various energies and travelling in various directions in the 
target, may reach the counter a rough estimate of the mean neutron energy 
and thence the mean detection efficiency for each counter position was made. 
The correction for efficiency change results in the upper curve of Fig. 3. The 
theoretical curve is also shown in the figure. The increase in neutron yield 
towards backward angles after a minimum at 90° results from the correspond- 
ing increase in the D(d, m)He* neutron yield at backward angles. In this case 
the effect is less pronounced than for a D(d, 2)He*® experiment since it is 
averaged over the scattered deuteron directions and energies in the target. 


8. Pulse Height Spectrum 
The pulse height spectrum from the fast-neutron counter produced by 
bombarding a thick heavy ice target with 1.5-Mev protons is shown in Fig. 4. 
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Fic. 4. The differential pulse height spectra from the neutron detector for 4-Mev and 
2-Mev monoenergetic neutrons and for the neutrons from a heavy ice target bombarded with 
protons of 1.5-Mev energy. 


For comparison the spectra produced by monoenergetic neutrons of 2-Mev 
and 4-Mev energy from the D(d, 2) He’ reaction are also shown. It is apparent 
that the neutrons from the ice target are comparable in energy to those from 
the D(d, n)He? reaction. 


4. Thin Target Yield and Angular Distribution 

A thin target containing 24X10'* deuterium atoms per square centimeter, 
approximately 70 kev thick for 1.5-Mev protons, was bombarded by 1.5-Mev 
protons and the neutron yield was measured as a function of angle. The 
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results are shown in Fig. 5. Since a large number of scattered deuterons 
escape from the back of the ice, it is not easy to make an estimate of the 
mean neutron-detection efficiency because it is difficult to estimate the neutron 
energy spectrum as a function of counter angle. Therefore no correction for 
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Fic. 5. The thin target yield and angular distribution function for 1.5-Mev protons. Curve 
(a), experimental results; curve (}), theoretical results. 


the change of efficiency with neutron energy was applied to the curve. The 
theoretical curve is also shown. The two curves show good agreement in 
shape and the factor of about two between the absolute values may be partly 
accounted for by the fact that in the theory for the thin target it was assumed, 
for the sake of simplicity, that all the deuterons scattered in the ice layer 
started out at the front face of the layer which gives an overestimate of 


the neutron yield. 


5. Effect of DsO Concentration 

If the neutron production from heavy ice targets is due to primary processes 
one would expect the neutron yield to depend linearly on the deuterium 
concentration in the target, whereas if the yield is due to secondary pro- 
cesses, as suggested here, one would expect the neutron yield to depend on 
the square of the deuterium concentration. This was checked by measuring 
the zero-degree neutron yield for 1.5-Mev protons using thick targets with 
various concentrations of heavy water in ordinary water. The results are 
shown in Fig. 6. The 6-Mev gamma-ray yield from the D(p, y)He® reaction 
is directly proportional to the deuterium concentration as expected since it 
arises from a direct interaction of the incident protons with the target nuclei. 
However, it is clear that the neutron yield increases as the square of the 
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Fic. 6. The effect of D2O concentration on the gamma-ray and neutron yields. The gamma - 
ray yield is plotted on a linear scale and the neutron yield is on a square-root scale. 


concentration. This leaves no doubt about the secondary nature of the 


neutron production. 


THEORETICAL ANALYSIS 

Calculations have been made of the yield, excitation function, and angular 
distribution of neutrons from thick and thin targets of heavy ice bombarded 
with protons on the assumption that they are produced by scattered deuterons 
colliding with other deuterons in the target and producing D(d, 2) He’ reactions. 

In the case of a thick target the calculations involved the following steps. 
The target was divided into six layers perpendicular to the proton beam 
direction as shown in Fig. 7(a), so that each layer has approximately one 
sixth of the incident proton energy of 1.5 Mev dissipated in it and in each 
layer the protons were considered to have a constant energy equal to the 
mean proton energy in the layer. The data of Sherr et a/. (1947) and Taschek 
(1942) on proton—deuteron scattering were used to calculate the energy and 
the angular distribution of the deuterons from each layer. The scattered 
deuterons from each layer were divided into overlapping cones, Fig. 7(0), 
symmetrical about the proton beam direction, with angular limits for each 
cone such that the difference in the deuteron energies at the outside and the 
inside of the cone was equal for all cones. It was assumed that all deuterons 
entering a cone have an energy equal to the mean energy of the deuterons 
entering the cone. Further, assuming Rutherford scattering and using the 
analytical expression for the energy and angular dependence of the proton— 
deuteron scattering cross section the number of deuterons scattered into each 
cone was obtained by analytical integration between the angular limits of the 
cone. To allow for the fact that the scattering does not follow the Rutherford 
formula for all energies and angles due to some nuclear interaction, the results 
from the Rutherford formula were multiplied by the ratio of the measured 
cross section to the Rutherford cross section. This ratio was interpolated for 
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Fic. 7. The division of the thick target into layers, cones, tubes, and sections as used in the 
theoretical calculations. 


the required energies from the results of Sherr et a/. (1947) and Taschek (1942). 
Each cone was subdivided into 12 “‘tubes’’, Fig. 7(6), each tube containing 
1/12 of the deuterons in the cone. All deuterons in each tube were assumed to 
travel down the center of the tube. For each tube the angle between this 
direction and the line joining the target to the detector was obtained geo- 
metrically. These angles were computed for the counter at 0°, 45°, 90°, and 
135° with respect to the incident beam direction. Finally each tube was divided 
into ‘‘sections’’ along its length such that the energy loss of the deuterons in 
traversing each section was approximately equal to the difference between 
the deuteron energies at the outside and the inside of the particular cone. In 
all there were 672 sections. The stopping power data for heavy ice of Wenzel 
and Whaling (1952) was used to compute the thickness and hence the number 
of target deuterons in each section of the target. The yield and angular 
distribution of D(d,)He*® neutrons has been expressed by Hunter and 
Richards (1949) in terms of Legendre polynomials which were used to provide 
the relevant neutron yields from each section in the direction of the counter 
for four counter positions at 0°, 45°, 90°, and 135°. The total neutron yield 
from the thick target in the counter direction can be obtained by summing 
over all sections of each tube, over all tubes of each cone, over all cones 
of each layer, and over all layers. 


“a 
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Table I gives typical neutron yield data at 0° for each of the six layers 
with the energy of the protons in each layer. 





TABLE I 
Proton energy, 0° neutron yield per steradian 
Layer kev per millicoulomb 
A 1500 to 1322 3.25105 
B 1322 to 900 7.02 X 105 
c 900 to 563 1.55 X 108 
D 563 to 338 1.08 X 104 
E 338 to 113 1.08 X10 
F 113 to 0 Neglected 





From this table the thick target yield of neutrons as a function of energy 
can easily be obtained by adding the zero-degree contribution from the 
bottom of the table to the layer which has an incident proton energy equal 
to the particular energy at which the total yield is required. The theoretical 
results are plotted in Figs. 2, 3, and 5 along with the experimental results. In 
principle the calculations for the thin target yield and angular distribution 
function were the same as for the thick target calculations except that only 
one target layer was considered; however, since a large fraction of the higher- 
energy scattered deuterons escape from the thin target layer the yield from 
the low-energy scattered deuterons was considered in greater detail than for 
the thick target calculation. 

The over-all agreement between theory and experiment leaves no doubt 
as to the secondary process responsible for the neutron production. 
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MAGNETIC CONTROL ON THE VARIATIONS OF THE 
CRITICAL FREQUENCY OF THE F, LAYER OF THE 
IONOSPHERE! 


R. G. Rastocr 


ABSTRACT 


The paper discusses the comparative influence of the true magnetic and 
smooth geomagnetic latitudes on the diurnal and latitudinal variations of the 
critical frequency of the F2 layer (foF2) at low latitudes. The diurnal variations 
of foF2 are shown to differ considerably at stations having the same geomagnetic 
latitude, but the discrepancies disappear when the true magnetic latitude is 
taken into consideration. The latitudinal variation of noon values of foF:2 is 
also shown to present discrepancies for low latitude stations in the geomagnetic 
latitudes plot, but on true magnetic latitude plots the points fall regularly along 
a smooth curve. 


Anomalous features in the diurnal and seasonal variations of the F» layer 
critical frequency were evident from the earliest ionospheric observations at 
Washington (Kirby, Berkner, and Stuart 1934) and Slough (Appleton and 
Naismith 1935). Longitudinal differences in fof: were first noted by Ranzi 
(1939), who showed that noon fo/2 were higher at Rome (41.8° N.) than at 
Washington (38.7° N.). Similar differences in fo, were noted between Washing- 
ton (38.7° N.) and Tokyo (35.7° N.) by Maeda, Uyeda, and Shinkawa (1942), 
and it was realized that the Washington values fitted well with the values at 
the Asian stations if plotted according to geomagnetic rather than geographic 
latitude. It was found (National Bureau of Standards, U.S.A. 1945) that the 
stations at Guam, Kwajalein Atoll, Christmas Island, and Huancayo, which 
were widely separated in geographic latitude and longitude, experienced similar 
diurnal variations of fo/». This resulted from the fact that these stations were 
at nearly equal geomagnetic latitude. Appleton (1946) was first to find con- 
siderable anomalies in the /2 layer at the equator, and expressed the opinion 
that better conformity could be achieved if the stations were arranged in 
accordance with geomagnetic latitude. A belt of low values of fo/2 circling 
the earth and centered roughly on the magnetic equator was also discovered. 
The geomagnetic distortion of fo/2 was further studied by Liang (1947), 
Bailey (1948), Lung (1949), Maeda, K. (1954), and Maeda, H. (1955) and 
they all confirmed the geomagnetic latitude control of fo/, distributions. 

The importance of true magnetic inclination in F, variations has been 
stressed by Rawer (1951). An effort was made by Burkard (1954) to examine 
the relative importance of geomagnetic latitude and magnetic inclination in 
the distribution of noon fo/2, but no significant difference could be found. 
Recently Rastogi (1959), while studying the influence of solar activity on the 
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geomagnetic control of the critical frequencies for the F, and F; layers, showed 
that the latitudinal distribution of fof: indicated less scatter when plotted 
against magnetic rather than geomagnetic latitudes. However, during a dis- 
cussion of the geomagnetic influences on the F2 layer at the recent U.R.S.1. 
Conference at Boulder, Martyn (1959) has expressed the view that it is better 
to plot fos: with respect to the smoothed geomagnetic equator rather than 
the actual magnetic equator. In this note it is intended to show that the 
latitudinal distribution of noon fo/2 in low and middle latitudes is controlled 
by the true magnetic latitude and not by the idealized geomagnetic latitude. 

Magnetic latitude » as referred to in this paper was derived from the 
magnetic inclination J of the station according to the relation 


tan wp = 1/2 tan I, 


which is analogous to the relation between geomagnetic latitude and geo- 
magnetic inclination. 

In Fig. 1 are compared the seasonal variations of fo¥: during the low sun- 
spot years 1952 to 1954 at pairs of ionospheric stations having approximately 
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Fic. 1. Seasonal variations of the critical frequency of the F; layer during the years of 
low sunspot number (1952-1954) at pairs of ionospheric stations having same geomagnetic 
but different magnetic latitudes. 


the same geomagnetic latitudes but different magnetic latitudes. The pairs of 
stations chosen and the corresponding latitudes are given below. 
In Fig. 2 are shown the diurnal variations of fof: (averaged for the year 
1954) at the pairs of stations Leopoldville-Madras and Bombay-Ibadan. 
Comparing observations at Bombay and Ibadan (both having geomagnetic 
latitude, = 10° N.), one finds that Bombay experienced much higher fo/2 
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Fic. 2. Yearly average diurnal variations of the critical frequency of the F; layer at pairs 
of stations Leopoldville-Madras and Bombay-Ibadan. Note that the ordinate scales for the 
two pairs are shifted with each other. 


than Ibadan. From a study of the ionosphere in the Asian Zone (Rastogi, un- 
published) it was found that Bombay was very near the zone of foF: peak 
during these years, whereas the various ionospheric phenomena observed at 
Ibadan, specially the association of E, with magnetic disturbances (Skinner 
and Wright 1957), justify listing Ibadan as an equatorial ionospheric station. 
Referring to Fig. 2 one finds that the diurnal variation of fo: at Ibadan had 
two maxima in the morning and evening hours, characteristic of any equatorial 
station. The diurnal variation of fo¥:; at Bombay and Ibadan can be explained 
only in terms of the magnetic latitudes of these stations. The magnetic latitude 
of Bombay (13° N.) places it in the group of middle latitude stations and that 
of Ibadan (1° S.) justifies it as an equatorial station. 
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Comparing observations at Dakar (# = 21°N.) and Yamagawa (® = 
20° N.) one finds that foF; at Dakar was much greater than at Yamagawa. 
The average values of fo/: at Dakar were comparable with those at other 
stations, e.g. Ahmedabad, Formosa, Calcutta, or Bombay, occupying the 
peak region for fo: Yamagawa was well outside the anomalous belt at all 
times (Rastogi, unpublished). These features can be reconciled with the fact 
that the magnetic latitude of Yamagawa is 26° N. while that of Dakar is only 
14° N. 

In the top curve of Fig. 1 one notices that fo: values at Leopoldville were 
as high as at the middle latitude stations Bombay and Dakar. Further, the 
values of fo/2 at Leopoldville and Madras were much different, even though 
the geomagnetic latitudes of the two stations are the same. Referring to Fig. 
2 one finds that diurnal variation of fo, at Madras had two maxima similar 
to that of Ibadan. The diurnal variation of foF: at Leopoldville showed only 
a single maximum, though the geomagnetic latitude of Leopoldville is only 
3° S. The non-equatorial type of the diurnal variation and the high values of 
foF: at Leopoldville can be interpreted if one notices that the magnetic lati- 
tude of the station is 19° S., and then one expects fof» at Leopoldville to behave 
as at other middle latitude stations. 

The latitudinal distribution of noon fy: during the low sunspot year of 
1953-54 is now plotted against both geomagnetic and magnetic latitude for 
comparison. In Fig. 3 separate plots are given for three periods in the year: 
the northern solsticial months, May, June, July, and August 1953; the southern 
solsticial months November, December 1953 and January, February 1954; and 
the equinoctial months September, October 1953 and March, April 1954. The 
points corresponding to stations in three different zones, East, West and 
Intermediate, are plotted in different symbols. 

A cursory glance of the figure shows that the points corresponding to low 
and middle latitude stations are much scattered in the geomagnetic latitude 
plots whereas these fall fairly along a smooth curve on magnetic latitude 
plots. The points most noticeably out of the smooth curve on geomagnetic 
plot correspond to Leopoldville, Lwiro, and Nairobi, Dakar, and Ibadan. The 
points for Leopoldville, Lwiro, and Nairobi represent the highest frequencies 
observed in the southern hemisphere but are very near to the geomagnetic 
equator. On the magnetic latitude plot these points are close to the region of 
peak fo». Similarly the point for Ibadan represents one of the lowest frequen- 
cies but on the geomagnetic latitude plot it lies at the region of northern 
foF, peak in any of the seasons. On the magnetic latitude plot the point corre- 
sponding to Ibadan occupies the region of minimum fo/2 and so low values 
of foF2 at Ibadan are explained. Further, the point for Dakar lies far removed 
from the smooth curve on geomagnetic latitude plot. The points for stations 
Dakar (@ = 20°), Panama (@ = 21°), Mauii (@ = 21°), and Okinawa (@ = 
20°) having about the same geomagnetic latitudes indicate the values of 
foF: differing by about 3 Mc/s. On the magnetic latitude plot all these points 
fall close to the smooth curve and the above discrepancies are removed. 

These facts indicate that the latitude distribution of the noon values of 
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foF2 within the anomalous equatorial belt may be controlled by the real 
magnetic latitude and not by the idealized geomagnetic latitude. 


The latitudinal distribution of fo/2 at high northern latitudes does not show 
any comparative preference between geomagnetic and magnetic latitudes; 
the points are equally scattered on the two plots. Further, the points corre- 
sponding to the West Zone (North America) and Intermediate Zone (Europe) 
lie close to the smooth line in the summer months, but in the winter months 
these lie along two different smooth lines. In other words, stations in North 
America and Northern Europe having the same geomagnetic (or magnetic) 
latitude record the same fo/2 in summer, but during winter months fo/» is 
greater at the American than at the European stations. It indicates that the 
factors predominantly controlling the F: ionizations at low latitudes may have 
comparatively less control at higher latitudes. 

It may be noted that the points corresponding to southernmost stations of 
the West Zone lie fairly close to the smooth line in the geomagnetic latitude 
plot but are far from the line in the magnitude latitude plot. This may indicate 
a comparatively greater control of the geomagnetic than the magnetic latitude 
in the South American Zone. The abnormal features of the /, layer at the 
South American stations have been discussed by Rastogi in a separate paper. 


The low values of fo/,: at these stations are suggested due to the abnormal 
distribution of magnetic dip lines in that region. 
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THE C-X SYSTEM OF N} 


P. K. CARROLL? 


ABSTRACT 
The C *2%-X *2¢ system of N+ has been investigated under high resolution 
in the region 1650-2000 A. Fifteen bands of the system have been analyzed and 
accurate rotational constants and other data for the C?Zé state have been 
derived. Anomalies in the rotational and vibrational structure are discussed. 
Some evidence in favor of a predissociation between the levels v = 2 and.v = 3 
of the C ?2% state is presented. 


I. INTRODUCTION 


The C — X bands of Nf form a fairly extensive system which stretches from 
about 1370 A to 2070 A. The system was discovered by Hopfield (1930) and 
later extended by Watson and Koontz (1934), who made a vibrational analysis 
and identified the lower state as the ground state of Nf. Setlow (1948) revised 
the vibrational scheme of Watson and Koontz and the correctness of the new 
analysis was verified by Miescher and Baer (1953) by means of the isotope 
effect. Setlow also made rotational analyses of a number of the bands and 
reported a marked anomaly in the B,-v curve for the upper state, C *2f, at 
low v values. Wilkinson (1956), using much higher resolution, studied a 
number of bands arising from the levels v’ = 3 and 4. 

The system presents a number of features of interest. Apart from the 
anomalously large B values reported by Setlow for the lower vibrational levels 
in the upper state, there appears to be an irregularity in the vibrational 
structure as well. Furthermore, the intensity of the system is particularly 
sensitive to excitation conditions. In discharges through helium which contains 
a trace of nitrogen, the bands for v' > 3 appear very readily with great 
intensity. In discharges through pure nitrogen, however, the bands are 
difficult to excite and lower vibrational levels in the upper state are favored. 
It was originally believed that the enhancement of levels with v’ > 3 in the 
helium discharge was due to collisions of the second kind between Het and 
nitrogen molecules in the ground state. Douglas (1952), however, suggested 
that an inverse predissociation could be responsible for the observed intensity 
distribution. 

The work described below was undertaken with the following objects in 
view: 

(1) To obtain accurate data on the rotational and vibrational structure of 
the C*Z* state. In Setlow’s work the resolution used was insufficient to 
separate the P and R branches and this leads to uncertainty and inaccuracy 


9 


in the analyses while Wilkinson has studied the levels v = 3 and 4 only. In 
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the present work attention was particularly focused on the bands with v’ < 3 
as these are difficult to excite and have been very inadequately studied. 
Extensive data for the lower state are available in the studies of Douglas (1952) 
and previous workers on the near ultraviolet bands of Nf. 

(2) To prove or disprove the presence of irregularities in the vibrational and 
rotational structure, with particular reference to the bands of low v’ values. 

(3) To elucidate, if possible, the process by which selective excitation occurs 
in helium discharges. In particular a predissociation limit between the levels 
v = 2 and v = 3 of the C*Z* state was sought as this would confirm the 


- 


Douglas theory of inverse predissociation at v’ = 3. 


Il. EXPERIMENTAL 

Three different sources were used in the experiments. The first source 
employed was a hollow cathode discharge tube with a cylindrical anode and a 
slot-shaped water-cooled aluminum cathode. When helium mixed with a trace 
of nitrogen was passed through this source, the bands with v’ > 3 were excited 
with great intensity. Discharges were also run in this tube using pure nitrogen 
at low pressure. It was hoped that in this case the excitation of the C *2+ 
state would be due to electronic collisions so that a more normal intensity 
distribution than that obtained in the helium discharge would result. Bands 
of the system with v’ = 0, 1, and 2 were in fact observed with this source, but 
the Lyman-Birge—Hopfield bands of Ne also appeared with considerable 
intensity and partially masked the Nf bands. In an attempt to get higher 
excitation a hollow cathode of the type described by Naudé and Hugo (1953) 
was employed. In this tube a thin-walled iron cathode is mounted in the 
center of a large glass bulb and by varying the current through the tube the 
temperature of the cathode can be controlled. It was found that as the tem- 
perature of the cathode rose, the pressure necessary to maintain the discharge 
decreased, while the degree of excitation of the nitrogen increased. The con- 
dition most favorable for the excitation of Nf spectra was reached at the stage 
when the cathode was just beginning to glow. At lower temperatures, bands of 
neutral nitrogen were enhanced whereas, if the cathode became red hot and if 
the pressure decreased, the discharge became white and omitted the spectrum 
of atomic nitrogen. The third source investigated was a mildly condensed 
discharge through pure nitrogen at low pressure. This source gave spectra 
rather similar to that given by the hot cathode discharge described above, i.e., 
the C-—»>X bands appeared in moderate strength while the spectrum of 
neutral nitrogen molecule was largely absent. 

The spectra were investigated in the third and fourth orders of the 3-meter 
vacuum spectrograph in the National Research Council of Canada. The optical 
arrangement described by Herzberg (1956) was used to eliminate the lower 
unwanted orders. Wavelengths were measured in the usual way against 
second-order iron arc lines. Comparisons made between different plates 
suggest that within a given set of measurements the wave numbers of lines 
are correct to about 0.08 cm~!. Between one set of measurements and another, 
however, approximately constant differences of as much as 0.20 cm~ have 
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been observed so that in absolute magnitude the wave numbers of lines may 
be in error by this amount. On the best plates, lines of equal intensity separated 
by more than 0.30 cm~! were clearly resolved, indicating a resolving power 
of about 170,000. Lines as close together as about 0.18 cm, although not 
fully resolved, could be identified as two distinct features and were measured 


as such. 


III. RESULTS 


In all, the fine structure of 15 bands was studied. These bands were, in order 
of decreasing wavelength: 3-10, 4-11, 5-12, 6-13, 0-6, 1-7, 2-8, 3-9, 4-10, 
1-6, 2-7, 3-5, and 0-2. The analyses of the bands are given in 


5-11, 
Table 


6-12, 


I. 


Table I includes data on four of the bands studied by Wilkinson (1956). 
Apart from the general desirability of having a set of homogeneous measure- 
ments on a band system, there are other reasons for including the new measure- 
ments of these bands. In the first place, over the region of the 3-10 and 4-11 
bands, the new measurements differ from Wilkinson’s by an amount varying 
from about 0.60 cm! at 50,500 cm! to 1.10 cm at 50,200 cm—'; the present 
measurements are believed to be correct to within the limits mentioned 
above. Furthermore, Wilkinson’s analysis of the 4-10 band seems to be 
incorrect. Finally the resolution in the present work appears to be about twice 
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that attained by Wilkinson. As a result of this, not only is there less blending 
of lines but there is also a partial resolution of the spin fine structure at higher 
N values of the 3-9 and 3-10 bands and in these bands, a small perturbation 
can be detected (see below). 

Transitions from the levels v’ > 3 were analyzed from plates taken with the 
helium—nitrogen discharge. A typical spectrogram taken with this source is 
reproduced in Fig. 1(a). On this and on similar spectrograms, transitions from 
the levels v’ = 2 and 1 are extremely weak. Curiously enough, however, 
transitions from v’ = 0 occur with moderate strength and it was in fact from 
plates taken with the helium-—nitrogen source that the 0-2 and 0-6 bands 
were measured. The effect is illustrated in Fig. 2, which shows the 0-2 band 
appearing with an intensity comparable to that of the adjacent 3-5 band 
while the 1-3 and 2-4 bands are completely absent. Presumably this intensity 
distribution is due to the fact that the 3-5 band, which has a small transition 
probability according to the Franck—Condon principle, is enhanced by the 
excitation conditions while the 0-2 band, although not favored by the excita- 
tion conditions, has a large transition probability. 

A comparison of Fig. | and Fig. 2 illustrates another interesting feature of the 
system. Whereas at longer wavelengths, for the Av = —8, —7, —6, and —5 
sequences, the vibrational structure of the system degrades in the wrong 
direction, i.e., in the sense opposite to the rotational structure, at shorter 
wavelengths, as for instance for the Av = —2 sequence shown in Fig. 2, the 
vibrational structure and rotational structure degrade in the same direction. 
At the intermediate stage, particularly in the Av = —3 sequence, the heads of 
the bands almost coincide resulting in a crowded and complex structure. 

Apart from the 0-6 and 0-2 bands the other transitions for v’ < 3 were 
analyzed from plates taken with the hollow cathode discharge through pure 
No». A typical spectrogram obtained from this source is reproduced in Fig. 1(d), 
which shows the 1-7 and 2-8 bands. 


IV. MOLECULAR CONSTANTS 

The band origins and rotational constants were derived by the usual 
graphical methods and the results are shown in Table II. The bands were 
chosen in such a way that each vibrational level in the upper state from v = 0 
tov = 6 was studied at least twice. In this way the correctness of the analyses 
could be checked and more accurate B values obtained. Table III gives the 
mean values of B for the upper state from v = 0 to v = 6. Here it should be 
pointed out that the accuracy and the extent of the data vary from band to 
band. In preparing Table III, therefore, extra weight was given to B values 
which were considered more reliable. Figure + shows a plot of B against v for 
the C*Z* state based on the data in Table III. It is seen that over the whole 
range of v investigated, the curve is non-linear with the departure from linearity 
being most pronounced at low v values. It is of interest to note that the 
dependence of B on v is quite different from that given by Setlow (1948), 
according to whom the B values increase rapidly with decreasing v reaching the 
value 1.62 cm~! for v = 0. 
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TABLE II 





3-10 ; 02 50374.15 
4-11 .4903 70% 50513. 
5-12 -479 5 6 50658. 
6-13 . 468 5 6 50806 . 
0-6 O11 806 51986. 
1-7 .5122 ‘ 52061. 
2-8 5043 . 765 52148. ¢ 
3-9 4972 7452 52246 . 2. 
4-10 489, : 52350. 
5-11 481, 70. 52459. 
6-12 4703 . 68: 52571. 
1-6 .508 5 805 54037 .: 
2-7 .5027 . 186 54090 .9, 
3-5 1,497) 825 60080. 2 
0-2 .509 . 88: 60225 .0o 


Average D’ 4.0X 10-6 
Average D” = 6.11076 


TABLE III 


Data for state C 2 


v 0 1 2 3 4 5 


1.489, 1.480, 
1946.1 1913.0 


B 1.5095 1.5096 1.5035 1.4970 
AG(v+3) 2051 . 1 2029.7 2004.3 1976.0 


° i 2 3 4 5 


Fic. 4. The variation of B with v in the C2 state. The shape of the curve for v < 2 is 
somewhat uncertain. Clearly the curve could be made to pass through the point By, in which 
case By would be considerably removed from its extrapolated position. Unfortunately B, is 
not known with sufficient accuracy to determine the course of the curve definitely. 
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Table III also gives the vibrational intervals in the C?Z state. These have 
been obtained from the band origins and are believed to be correct to about 
0.30 cm~!. As in the case of the B values, the dependence of AG(v+3) on v is 
non-linear over the whole range investigated. A graph of AG(v+ 4) against v 
gives a curve which is qualitatively very similar to the B, versus v curve. In 
this case the curvature is again most pronounced at low v values. Because 
of this, it is not possible to represent the vibrational levels in the C*Z state 
by the orthodox vibrational formula. Furthermore, it is not possible to 
determine an accurate value of v,. However, one can derive the approximate 
value 64,628 cm! which should not be in error by more than +2 cm—!. An 
accurate value of vo can of course be obtained by combining the present 
results with the known data for the ground state of Nf. The value obtained is 


voo = 64,542.0 cm—'. 


V. ROTATIONAL PERTURBATIONS 

Three small but definite perturbations are observed in the rotational 
structure of the bands. The most obvious of these occurs in the level v = 3 
of the C *X state. It has been mentioned that only at higher N values (> 24) 
is it possible to detect a partial resolution of the doublet structure of the P 
and R branches. Nevertheless in the 3-10 band the R(16) and P(18) lines 
obviously consist of two incompletely resolved components while R(17) and 
P(19) each consist of two distinct lines. The effect is also observed in the P 
branch of the 3-9 band; it is not detected in the R branch due to the blending 
of lines. It is present in the 3-5 band also but here the effect is less noticeable 
due to the decrease in resolution at the shorter wavelengths. 

It is obvious that the appearance of spin structure is due to an increase in the 
splitting of the levels N = 17 and N = 18 in v = 8 of the C*Z state. The 
average value of the observed splitting is 0.21 cm~ for the lines arising from 
N = 17 and 0.31 cm for the lines from N = 18. The higher-frequency 
members of the doublets fit smoothly into their appropriate P and R series; 
clearly therefore only one component of the spin doubling in the state C°*Z+ 
is affected. One would normally attribute such a perturbation to a ?II, state. 
Only one such state is known for Nf, namely the upper state, A *II,, of the 
Meinel infrared bands. This state, however, very probably converges to the 
limit N(4S)+Nt(P) at 70,358 cm and if this is so, it cannot be responsible 
for the perturbation in C?Z which is observed at 71,084 cm~'. Furthermore, 
according to Mulliken (private communication), no other II, state is expected 
theoretically in this energy region. However, the configuration 


KK (a,2s)?(ou2s) (4y2p) (0 ,2p) (4 ,2p), 


which produces the C *2+ state, also gives rise to a *Z> state and this may 
well be responsible for the observed perturbation. A consideration of the 
selection rules shows that in this case only one spin component of the C2? 
state will be affected. 

A second perturbation is observed in the level vy = 2 at N = 23. It is found 
that the P(24) and R(22) lines in the 2-7 and 2-8 bands are displaced to higher 
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energies by about 0.30 cm—!. The perturbation may also effect transitions from 
higher WN levels, but the effect is smaller and its detection less certain. In 
contrast to the perturbation in v = 3, no spin splitting is observed in the 
present case so that both spin components of the upper state must be equally 
displaced. A similar perturbation is found at N = 24 in the level v = 6. 
Again both spin components are affected and in this case the level is shifted to 
lower energies by about 0.40 cm~'!. These perturbations could be due to the 
states 427 or *A, (case b) which arise from the same electron configuration as 
C** and which should lie at about the same energy. It is also possible 
that the perturbation in v = 2 is due to some high levels in the B 22 state. 


VI. VIBRATIONAL STRUCTURE OF THE C?2z STATE 

The band measurements given by Watson and Koontz (1934) and Tanaka 
(1953) cannot all be represented by a vibrational formula of the usual kind. 
In particular if one uses a formula which represents the higher levels of the 
C *Y state satisfactorily, it leads to considerable error for the levels v = 0, 1, 
and 2. Setlow (1948) believed that rotational perturbations were produced in 
these lower levels by a 7II state and that this effect, combined with the limited 
accuracy of the available measurements, could account for the observed 
anomaly. As stated in Section IV, however, the data based on band origins 
leave no doubt that the AG(v+3) versus v curve is abnormally shaped. The 
curve shows that the vibrational intervals at lower v values are smaller than 
predicted by an extrapolation of the intervals at higher v values. If we consider 
this effect as being due to a perturbation, it would correspond to an upward 
shift of the lower levels from their original unperturbed positions. The shift 
falls off rapidly, but smoothly, with increasing v; there is no indication of a 
discontinuity at v = 3 as mentioned by some of the earlier workers. 

The variation of B with v, described above and shown in Fig. 4, can be 
correlated directly with the anomalous vibrational structure. It is particularly 
interesting to note that the B versus v and AG(v+3) versus v curves follow 
the same pattern, such that where the vibrational intervals are abnormally 
small the B values are also abnormally small. This is a feature of vibrational 
perturbations and occurs for example in the upper state of the AgH bands 
(see Bengtsson and Olsson 1931; Geré and Schmid 1943) and in the B ?3 state 
of Nf (Douglas 1952). 

Douglas (1952) has interpreted the unusual shape of the AG curve for the 
Bz state in terms of an interaction between the potential curves for the 
B *Z and C 2 states. He suggested that the interaction was such that a 
repulsion takes place between the two curves leading to a decrease in the 
equilibrium internuclear distance for the C *2 state. At the time this appeared 
to be supported by the large Bo value given by Setlow (1948). We know now 
that in fact the By value is abnormally small. But this does not invalidate 
Douglas’s suggestion. It may well be that the repulsion between B?Z and 
C* is such that upward displacement of the lower part of the CZ state is 
much more pronounced than any lateral displacement of the potential mini- 
mum. This could lead to anomalies in the vibrational and rotational structure 
of the kind described above. 
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VII. PREDISSOCIATION 


When the system is observed in helium—nitrogen discharges, the level v’ = 3 
is excited with great intensity. Higher vibrational levels are also excited, but 
to an extent that falls off fairly rapidly with increasing v’. Bands with v’ < 3 
are absent or weak. It was suggested by Watson and Koontz (1934) that this 
curious intensity distribution was due to resonance excitation of Nf in collisions 
between ionized helium atoms and nitrogen molecules in the ground state. 
There are, however, objections to this theory. In the first place Tanaka (1953) 
has pointed out that in this process maximum excitation should occur at v’ = 4 
and not at v’ = 3 as observed. In the second place it is difficult to see how levels 
as high as v’ = 9 and 10 can be excited as they lie considerably above the 
ionization potential of helium. 

Douglas (1952) has proposed that the intensity distribution in the helium 
discharge should be attributed to an inverse predissociation from a state 
arising at the dissociation limit N(4S)-+Nt+t(P). In this process, energy con- 
siderations show that the level vy = 3 should be preferentially excited. The extra 
energy required to excite the higher vibrational levels could be supplied by 
the relative kinetic energy of the “4S atom and *P ion as they approach along 
the potential curve of the repulsive state. The process of inverse predissociation 
would be expected to be important in discharges such as the helium—nitrogen 
discharge where it is very likely that a large number of nitrogen atoms and 
atomic ions is present. On the other hand, as pointed out by Douglas, if the C 
state is excited directly by electron collisions and if the density of atoms and 
atomic ions is small, one would expect a decrease in the intensity of the bands 
due to predissociation to be observed at v’ = 3. It was hoped that in present 
work the spectra taken with the hollow cathode discharge through pure 
nitrogen would provide some definite evidence on this point. In particular a 
breaking off in the rotational structure of bands from the level v’ = 2 was 
sought as this would settle the question decisively. However, the branches 
in the 2-7 and 2-8 bands appeared to behave in a normal way without any 
sudden breaking off or drop in intensity. This of course does not mean that 
predissociation is absent because it could well occur above the point at which 
the branches fade in the normal way. The highest rotational level observed 
was N = 31, which is about 517 cm~ below N = 0 of v’ = 3 and it is quite 
possible for predissociation to occur above this point. In fact, from the appear- 
ance of the vibrational structure of the bands, it is highly likely that this is the 
case. Thus Fig. 3 shows, side by side, spectra of the Av = —5 sequence taken 
with the helium—nitrogen discharge and the cold, hollow cathode discharge 
through pure nitrogen. In the latter case it is clearly seen that bands with 
v’ = 2, 1, and O are well developed while only a faint trace of bands with 
v’ > 2 can be detected. It is very difficult to interpret this behavior in any 
way other than in terms of a predissociation. 

It is evident that the predissociation must be a weak one. The fact that the 
bands with v’ > 3 appear strongly in the helium discharge suggests that the 
lifetime with respect to predissociation must be at least of the same order as 
the lifetime with respect to radiation. Further evidence that the predissociation 
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is weak is provided by the fact that in direct electronic excitation the bands 
with v’ > 83 still appear, although very feebly. 

The configuration N(*S)+N+(P) gives rise to 12 states of species =f, 
z, 8, and II, each of which can occur as doublets, quartets, and sextets. 
Of these *2* and *=* are undoubtedly the upper and lower states respectively 
of the first negative bands, while the “II, state is very probably the upper state 
of the Meinel bands. Bearing in mind that g states cannot predissociate u 
states and that predissociation of a doublet state by a sextet state would be 
much more strongly forbidden than by a quartet state, we are led to the 
conclusion that the predissociating state is probably 42+ or ‘II,. According 
to Mulliken (private communication) the *2* state should be stable and lie 
1 to 2 volts below the C *2 state. We therefore conclude that the state respon- 
sible for the predissociation is very probably “II,. Thus it is seen theoretically 
that a predissociation of state C?Z must be a forbidden one—which is at 
least consistent with the experimental observation that it is weak. 


I am very grateful to Dr. G. Herzberg and Dr. A. E. Douglas for their valuable 
advice and suggestions in the course of this problem, and to Professor R. S. 
Mulliken for providing important information on the electronic structure of 
Nf. I am also indebted to Professor T. E. Nevin for his interest in the work 
and for reading the manuscript and to Sister Cabrini for her assistance in the 
measurement and analysis of several of the bands. Finally I should like to 
thank Mr. J. Shoosmith for his expert help in photographing the spectra. 
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ON MAINTAINING AN OBJECT AT A UNIFORM TEMPERATURE 
C. M. CARMICHAEL AND A. D. MISENER 


INTRODUCTION 

When it is desirable to maintain an object at a constant temperature and 
without temperature gradients in the object itself, one common practice is to 
suspend it in a uniformly heated enclosure. If high precision is required, par- 
ticularly with respect to the absence of temperature gradients in the object, 
considerable attention must be paid to the uniformity of the temperature of 
the enclosure as well as its constancy. We have undertaken to test experi- 
mentally what external conditions are necessary to maintain a temperature 
uniform to within 0.003° C over metal objects of reasonable size. 

Interest in this problem was aroused by the work of Thompson (1956) in 
developing a new pendulum apparatus for determining the relative accelera- 
tion of gravity. By placing his penduli in a thick-walled, insulated, and heated 
metal case he could maintain their average temperature constant to the re- 
quired +0.01° C, but there was a gradient of 0.05° C along their length and 
their average temperature was about 0.08° C below the average of the case 
as read by 20 thermistors placed symmetrically over its inner surface. 

Two types of experiment were performed. In the first type the objective was 
to produce as uniform a temperature in the suspended bodies as possible using 
ordinary laboratory techniques. In the second type, constant temperature 
gradients were produced on the enclosure and the resultant gradients on the 
suspended objects were measured. Since no enclosure can have an absolutely 
uniform temperature all over its surface, it is hoped that these latter results 
will be of some value in estimating the gradients likely to exist on enclosed 
specimens where direct measurement is impracticable. 


APPARATUS AND PROCEDURE 

The chamber used consisted of an 18-in. length of 6-in. diameter steel 
tubing of }-in. wall thickness onto which were bolted flat steel end plates using 
0 rings for vacuum seals. The outer surface was wrapped with 28-gauge copper 
wire spaced so that adjacent wires were about 3 inch apart and divided into 
six separate heating coils each covering about one-sixth of the surface. Insu- 
lation for the chamber consisted of a layer of fiber glass 1 in. thick surrounded 
by 6 inches of foamed polystyrene which was strong enough to act as the 
supporting base. 

Heat was supplied to the system by passing direct current through the 
heating coils under the control of a relay operated by an a-c. bridge circuit. 
One arm of the bridge consisted of a thermistor, fastened to the outside of the 


Can. J. Phys. Vol. 37 (1959) 





890 








NOTES 891 


steel chamber, which threw the bridge off balance as its temperature changed 
from any predetermined value in the range 20° C to 100° C. The phasing of 
the circuit was such that a change in temperature of the control thermistor 
of plus or minus 0.001° C would turn the heating current off and on respectively. 
Because of the time lag in the heat flow through the chamber walls the tem- 
perature of the inner surface showed a fluctuation of +0.005° C periodic 
with the heating cycle. 

It was found that while any point on the inner surface remained constant 
to the nearest 0.01° C it was not possible to prevent differences in temperature 
of about 0.1° C over the surface. This was not uniform enough, so a second 
container was placed inside the chamber to act as the uniform temperature 
enclosure. Of the several enclosures tried the most successful was a 12-in. long 
copper cylinder, 4 in. in diameter, and having a §-in. wall thickness. This was 
supported on cork feet in the center of the steel chamber. Inside the enclosure 
objects with various shapes and conductivities were suspended on nylon 
threads. Best results were obtained with an object consisting of a cylindrical 
bar of aluminum 9 in. long and 13 in. in diameter (Fig. 1). 


STEEL CHAMBER —————, 
COPPER ENCLOSURE — / 









ALUMINUM 
OBJECT , 


PLAN THROUGH -A- A- 


¢ 7 - THERMISTOR 
-'™ - SCALE 


-—— SECTION THROUGH 
CHAMBER 


Fic. 1. Diagram of chamber, enclosure, and object showing location of thermistors. 


Temperatures were measured using 20 Veco 32A38 bead thermistors, 14 on 
the enclosure and 6 on the object. These were calibrated against each other 
for several ranges of about 1° C using a Beckman thermometer which could 
be read to the nearest 0.001° C. For calibration the thermometer, with the 
thermistors taped around the bulb in a single layer, was placed in a small 
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mercury bath which was in turn immersed in a water bath. The resistances of 
the thermistors were measured with an H.W. Sullivan Ltd. Wheatstone bridge 
with which the resistance could be measured to the nearest 0.1 ohm. This 
corresponds to reading the temperature, for the thermistors used, to +0.001° C. 
In practice it was found that the time necessary to balance the bridge for this 
precision was too long so that measurements of temperature were not repro- 
ducible to better than about 0.002° C. Best thermal contact between ther- 
mistor and metal was obtained by drilling small holes in the wall of the en- 
closure and object and inserting the thermistors in these along with a drop 
of mercury. 

Measurements were made with a range of air pressures in the system. It 
was found that for ordinary pressures the effects of convection made it difficult 
to adjust the power supplied to the individual heating coils so as to produce the 
uniform temperature condition, particularly when using odd-shaped objects. 
Thus for most of the work the pressure was maintained at 0.01 to 0.02 mm of 
Hg. 

The general measurement procedure was to read the resistance of each of 
the 20 thermistors at hourly or half-hourly intervals each day during a period 
of 1 week or more for each series of measurements. After any adjustment of 
power to the heating coils or relocation of the thermistors the system was 
left for 1 day or longer before readings were taken. 


RESULTS 
In investigating the equilibrium temperature over the surface of an object 
in relation to the temperature over the surface of the enclosure in which it 
was suspended, two types of measurements were made. 


Uniform Temperature 

With careful adjustment of the power supplied to the individual heating 
coils the temperatures over the copper enclosure showed a standard deviation 
of about 0.003° C, which was about one and one-half times the reading error. 
Under these conditicns the temperatures over the suspended object showed a 
standard deviation of about 0.002° C or about the same as the reading error. 
The temperature of either enclosure or object did not, however, remain con- 
stant. There was a random drift of about 0.002° C per hour with the object 
showing a time lag of several hours. Thus over a period of days the control 
circuit used could not maintain the temperature of the object constant to 
better than 0.01° C. Because of this drift it was difficult to compare the equi- 
librium temperature of the object with that of the enclosure to the limit of 
temperature-measuring accuracy. The nearest to an equilibrium condition 
was reached when the random drift was small and roughly periodic about a 
constant temperature over a period of several hours so that the temperature 
of both the enclosure and the object showed a similar periodic change. Under 
these conditions the object reached the same temperature as the enclosure. 
For example, in the same region of the temperature scale as Thompson worked, 
during a period of 4 hours the average of 56 readings of the 14 thermistors on 
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the enclosure was 41.1328° C with a standard deviation of 0.0033° C and for 
the same period the average of the 24 readings of the 6 thermistors on the 
aluminum object was 41.1332° C with a standard deviation of 0.0025° C. 


Temperature Gradients 

The second set of measurements, where a temperature gradient was set up 
along the length of chamber, enclosure, and object, were more valuable in 
determining the usefulness of this type of apparatus for holding an object at 
a constant and uniform temperature. The power supplied to the heating 
coils on the chamber was adjusted so as to produce temperature gradients, 
from top to bottom of the chamber, of up to 2° C, both with the top hotter than 
the bottom and vice versa. There was no appreciable difference between the 
measurements for the top hotter than the bottom and those for the top colder 
than the bottom for pressures of the order of 0.01 mm Hg, indicating heat 
transfer was mainly by radiation and conduction. For the dimensions and 
materials used in the apparatus there was a roughly 10-fold reduction in 
temperature gradient from chamber to enclosure to object as indicated in 


Table I. 


TABLE I 


Temperature differences over outer chamber, inner enclosure, and suspended object 





Temperature (°C) at: Temperature differences 
- between top and bottom, 
* 








Location of thermistors Top Middle Bottom 
Outer steel chamber 39.8 41.0 41.9 2.1 
Inner copper enclosure 40.68 40.86 41.01 0.33 
Suspended aluminum cylinder 40.859 40.864 40.870 0.011 





Tuompson, L. 1956. Ph.D. Thesis, University of Western Ontario, London, Ontario. 


ReEcEIVED Marcu 31, 1959. 
DEPARTMENT OF Puysics, 
UNIVERSITY OF WESTERN ONTARIO, 
LONDON, ONTARIO. 


NEGATIVE AND OSCILLATORY LONGITUDINAL MAGNETORESISTANCE 
R. BARRIE 


The purpose of this note is to report the result of numerical calculations 
of the longitudinal magnetoresistance of semiconductors and semimetals, for 
the case in which the scattering of the electrons is due to the acoustic modes 
of the lattice vibrations and the magnetic field is so large that the quantiza- 
tion of the electron orbits is important. The conductivity, o, in the presence 
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of the magnetic field, H, is related to that in zero field, oo, by the expression 
(Barrie 1957; Argyres 1958) 


(1) c/ay = ~3{In(Lte)}" J es) 2 fen) dx 
H’/2 ox 
where 


—1 
g(x,H’) = Do en DH} be oy" | 


n 


and 


f(x—nw) = (7 +1). 

In this equation, x = «/kT, H’ = heH/kTm*c, no = £0/RT, and ny = ¢x/kT? 
f> and ¢q are the Fermi energies in the absence and presence, respectively, 
of the magnetic field; « = (w+ 3)(heH/m*c)+¢’ is the electron energy, with n 
the oscillator quantum number associated with the transverse motion of the 
electrons and ¢’ the translational energy in the direction of the fields. The 
sums over the non-negative integer ” are from zero to that maximum value 
of m which leaves {x—(n+})H’} positive. For the assumptions involved in 
the derivation of (1), the reader is referred to the papers quoted above. The 
function g(x, H’) suffers discontinuities, so no assumptions concerning the 
step-function nature of f(x—n,) were made. 

In addition to (1), a relation is needed between ny and no. If the assumption 
is made that the number of electrons in the band is not altered by the magnetic 


field, then 


(2) Fy (no) = (H'/2) D0 Fy{na— (n+3)H’} 


n=0 


where 


os 
F,(n) = J x*(e7 "+1) dx. 
0 

The general nature of the dependence of n/n and of the magnetoresistive 
ratio Ap/po (= (a0/0)—1) on H’ as np is increased is illustrated in Figs. 1, 
2. wand. 

If there are unfilled levels below the band being considered, the magnetic 
field can, through its alteration of the energy-level distribution, change the 
number of electrons in the band. Equation (2) would then have to be suitably 
modified. For semimetals, (2) will be valid, but for semiconductors it is 
not. The value 40 for no would represent a rather extreme case of a semimetal 
and the value —20, a rather extreme case of a semiconductor. If the alteration 
of the levels in the lower band is less marked than in the upper, the suitable 
modification of (2) will lead to an increase in Ap/po. This would, for example, 
be the case in an intrinsic semiconductor in which the hole effective mass is 
much greater than the electron effective mass. 

In comparisons of theory and experiment, it is often assumed that the 
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Fics. 1-3. Plot of ny/no and Ap/po as functions of H’. — — — ny/n0o; ———— Ap/po. 
Fic. 1. no = — 20. Fic. 2. no = 0.92. No calculations were made for H’ < 0.5. Fic. 3. 


no = 40. The isolated points are for H’ = 2, 4, and 5. No curve is drawn through these points 
as the oscillations became very rapid. 


oscillations can be approximated by an analytic function, e.g. Broom’s use 
of the theory of Davydov and Pomeranchuk for impurity scattering (Broom 
1958). This approximation is made in the theory when H’ ny (Argyres), 
a condition not always satisfied in the experiments on semimetals. For the 
problem discussed here, the oscillations are represented analytically by 
—cos {2(anq/H')—(2/4)}. The numerical analysis was undertaken to check 
this approximation. If the reciprocals of the reduced fields H’ at which Ap/po 
is a minimum are plotted against the appropriate series of integers, a straight 
line should be obtained, whose slope is 1/nq and whose intercept on the 
number axis is —1/8. When this graph is drawn for no = 40 (the minimum 
at H’ = 32 being that for the integer unity) a straight line is obtained of 
slope 1/40. The intercept on the number axis, however, is —1/4. The analytic 
approximation thus indicates the correct period of the oscillations, but the 


wrong phase. 

The most interesting result of the calculations is the appearance of a 
negative Ap/po. This is in general agreement with experiments on indium 
arsenide (Sladek 1958). In this material, Ap/po oscillates between positive 
and negative values. To obtain better agreement, the calculated mean value 
of Ap/po during one oscillation would have to be increased. This could perhaps 
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be achieved by one of two modifications required in the theory. First, the 
effect of energy-level broadening was neglected in the derivation of (1). This 
would remove the discontinuities in g(x, H’) and thereby tend to smear out 
the oscillations. It might also raise the mean value of Ap/po. Secondly, the 
scattering of the electrons in indium arsenide, at the temperatures involved 
in the experiment, is due to ionized impurities; the equation corresponding 
to (1) would have to be derived for this type of scattering. 
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SUGGESTED MODIFICATIONS TO A METHOD FOR THE 
DETERMINATION OF THE ABSOLUTE AMPERE* 


M. ROMANOWSKI AND R. BAILEY 


While the absolute unit of electrical resistance is determined by several 
different methods, the unit of current is determined exclusively by means of 
current balances (coaxial coils or torque dynamometer types). Silsbee (1949), 
after having examined the results obtained from various current balances, 
commented as follows: ‘There is an evident need for the development of some 
radically different method for measuring current in absolute amperes that 
will not be subject to the same type of systematic error.”’ 

Only one such method is known to exist. It is based on the fact that we can 
express the magnetic induction B, produced in a high-quality magnet, by 
means of two equations from which the induction vector is eliminated, thus 
obtaining a relation between current and quantities which can be expressed 
in terms of mechanical units. 

The first of these equations, F = BIJ, expresses the force F resulting from 
the action of the field of magnetic induction B on a conductor of length / 
carrying a current J (Fig. 1). The second equation, BK = MI, is deduced 
from the relationship of the e.m.f. induced in a search coil? rotating between 

*Issued as N.R.C. No. 5233. 


{The inductive constant K of this coil has to be determined beforehand, by an independent 
operation, in terms of the value of a computable self-inductor. 
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the poles of the magnet, and the e.m.f. induced in the secondary of an inductor 
of known value M. The primary of this inductor is energized with the same 
current J and this current is reversed synchronously with the 180° flip of the 
search coil. The two e.m.f.’s are compared with a ballistic galvanometer. 

This method has been used by Dupouy and Jouaust in 1935. In 1936, 
Briggs suggested that its precision might be increased by substituting a 
balance equipped with a freely suspended rectangular coil in place of the 
Cotton-type balance used by Dupouy and Jouaust. 

The difficulties associated with the reversals of current, in addition to those 
already existing in the remaining parts of the experiment, have prevented 





Fic. 1. Fic. 2. 


this method from becoming a serious rival to the current balance. In the 
course of a re-examination of the problem we came to the conclusion that 
the reversals of current and the difficulties associated with the calibration 
of the mutual inductor MV (in terms of a computed standard) can be avoided 
by using a homopolar generator instead of a search coil. In fact, this substi- 
tution transforms the method into a completely direct-current one. 

The values of the e.m.f. E (Fig. 2), generated between the rim and the 
center of a conductive disk D (thin metal platings on glass) rotating in a 
magnetic flux ®, is expressed by E = &n where n is the number of revolutions 
per second. This electromotive force can be balanced by the drop of potential 
FE’ across a resistor X (measured in terms of absolute ohm) through which 
flows the current J. When E = E’, the relation becomes @2 = BSn = RI, 
S being the area of the disk. 

With a good magnet it is relatively easy to obtain a field of a density of 
0.5 Weber/m?*. In a disk of 7-cm radius the flux @ will be equal to 7.510% 
Weber. Using = 80r.p.m. the induced e.m.f. / will be of the order of 
10 mv. This corresponds to a drop of potential E’ created by a current 
F= 0:1. 4 mi@2Q, 

The equation BSn = RI, combined with F = BIJ, gives J* in terms of 
S, # R,- andl. 

In an actual experiment the position of the brush at the axis of the disk 
cannot be known with sufficient accuracy; this contact can be replaced by 
another brush contact on a second disk of much smaller diameter. The resulting 
loss of flux will be small. 

An alternative technique would be to use a radius sector rotating in a 
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conducting ring. It may be necessary to try both experimentally before taking 
the final decision. 

By using two disks which rotate in opposite directions it is possible to 
create more degrees of freedom for analyzing secondary phenomena. The 
e.m.f. E may be opposed to the e.m.f. created by another homopolar generator 
of computable characteristics. This procedure, although probably very difficult, 
should not, a priori, be excluded. The technique of computed generators has 
been demonstrated from the work done on the Lorenz apparatus (Smith 
1914). 

The problems connected with creating a stable and reproducible field with 
the required homogeneity over the surface of the disk are difficult but not 
insurmountable with modern techniques based on the use of proton resonance. 

We may summarize by observing that the Briggs freely suspended coil is 
capable of yielding a high accuracy in the operation of the ‘“‘weighing”’ of the 
current. If, in addition, we substitute the homopolar generator for the flip-coil 
technique, we should obtain a higher accuracy for the second operation by 
eliminating the necessity for a well-known mutual inductance of high value 
and the problems related to the reversals of current. 

These two operations appear to constitute a method of sufficient accuracy 
to warrant an experimental trial. This would provide an independent check 
of the present classical balance technique for determining the absolute ampere. 
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